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ABSTRACT 


The  theory  of  operation  of  Laser  Doppler  Velocimeters  is 
discussed  for  both  reference  beam  and  differential  Doppler 
modes.  Current  entrainment  models  for  velocity 


distributions  in  vertical  axisymmetric  turbulent  buoyant 
jets  in  a  quiescent  ambient  are  reviewed.  A  specific  Laser 
Doppler  Velocimeter  is  used  to  measure  the  velocity 
distribution  found  in  a  vertical,  axisymmetric  turbulent  jet 


discharged  into  a  quiescent  ambient  and  the  results  are 


compared  to  theory. 
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I. 


INTRODUCTION 


The  study  of  the  characteristics  of  turbulent  jets  and 
plumes  has  become  an  item  of  widespread  interest  as  the 
concern  for  the  environment  has  grown.  An  understanding  of 
how  the  velocity  profile  varies  throughout  che  jet  is 
central  to  the  development  of  knowledge  of  the  behavior  of 
the  jet  with  respect  to  heat  and  mass  transfer,  both 
subjects  of  intense  interest  to  the  engineer  or  scientist 
who  must  predict  the  impact  of  a  discharge  on  the 
environment.  Further,  changes  in  temperature  or 
concentration  may  someday  have  significance  in  non-acoustic 
antisubmarine  warfare  applications. 

Researchers  have  used  a  variety  of  techniques  to  study 
turbulent  jets  and  plumes.  Qualitatively,  the 
characteristics  may  be  disclosed  through  photographic 
studies  involving  the  injection  of  tracer  dyes.  For  the 
purposes  of  determining  the  widths  of  the  plumes  and  their 
trajectories,  this  technique  is  useful,  but  precise  velocity 
information  or  numeric  data  regarding  the  temperature  or 
concentration  profiles  within  the  fluid  are  not  available 
from  this  method.  Velocity  information,  as  well  as 


temperature  can  be  determined  using  hot  wire  anemometry,  but 
this  technique  requires  that  a  probe  be  inserted  in  the  flow 


in  order  to  obtain  the  measurement.  The  presence  of  the 
probe  is  intrusive  in  nature  and  will  invariably  alter  the 
flow  characteristics  of  the  system.  Similar  objections 
arise  when  the  use  of  rakes  of  thermocouples,  conductivity 
probes,  or  any  similar  device  is  contemplated.  For  years, 
however,  no  alternative  existed  and  researchers  were 
required  to  make  do,  minimizing  the  effects  of  the 
measurement  devices  if  possible. 

In  1964  Cl]  [2],  it  was  noted  that  laser  light  in  a  flow 
in  which  particles  were  present  was  scattered,  and  that  the 
scattered  light  was  shifted  in  frequency  by  an  amount 
proportional  to  the  velocity  of  the  scattering  particles. 
Subsequent  work  developed  the  technique,  known  as  Laser 
Doppler  Velocimetry  or  Laser  Doppler  Anemometry,  with  the 
result  that  there  is  now  a  reliable,  non-invasive  method  of 
measuring  fluid  velocity  in  a  variety  of  environments. 

For  the  purposes  of  the  present  work,  a  commercial  Laser 
Doppler  Velocimeter  was  obtained.  The  device  was  interfaced 
with  a  Hewlett-Packard  9826  microcomputer  for  data 
acquisition,  a  tank  was  constructed  and  a  suitable  nozzle 
was  installed  for  the  introduction  of  a  vertical  submerged 
jet  of  water  into  an  ambient.  Here  the  intent  was  to  gather 
data  on  a  jet  discharging  to  the  quiescent  medium  of  the 
tank,  and  to  develop  the  measurement  procedures,  interfacing 
programs,  and  data  reduction  facilities  necessary  to  analyze 
the  behavior  of  the  jet.  The  results  were  compared  to  those 


.2 


obtained  by  earlier  researchers  under  similar  conditions 
with  the  objective  of  validating  the  technique  and  ensuring 
proper  functioning  of  the  system. 

Data  was  collected  at  a  variety  of  flow  rates  in 
quantities  sufficient  to  constitute  a  valid  statistical 

base.  Raw  data  gathered  by  the  microcomputer  in  conjunction 
with  the  LDV  electronics  was  transmitted  via  telephone  line 
to  a  Digital  Equipment  Corporation  VAX  11/780  computer, 
which  was  large  enough  to  perform  the  required  manipulations 
to  reduce  the  data  to  usable  form.  Results  were  expected  to 
be  consistent  with  those  obtained  by  other  researchers 
utilizing  other  methods. 

Specific  information  to  be  gathered  and  compared  was: 

1.  The  jet  centerline  velocity  as  a  function  of 

streamwise  coordinate  and  the  velocity  of  discharge  at 
the  nozzle  axis. 

2.  The  variation  of  jet  characteristic  width  with 

streamwise  coordinate. 

3.  The  rate  at  which  the  ambient  fluid  is  entrained  into 
the  jet,  as  characterized  by  the  entrainment  constant. 

Information  which  was  also  gathered  in  order  to 
demonstrate  the  capabilities  of  the  system  included: 

1.  The  mean  velocity  component  normal  to  the  streamwise 
coordinate. 

2.  The  turbulence  intensities  in  the  radial  and 

streamwise  directions. 
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3.  The  mean  Reynolds  stress  at  the  point  of  measurement. 

The  measurements  were  made  using  the  LDV  to  gather  data 
in  first  the  streamwise  direction,  then,  after  rotating  the 
plane  of  the  beams,  in  the  direction  normal  to  the 
streamwise  coordinate.  No  attempt  was  made  to  take 
instantaneous  readings  in  the  two  directions  simultaneously, 
as  this  model  LDV  possesses  no  capability  to  do  so. 

As  the  data  was  gathered,  the  velocity  information  was 
monitored  on  the  screen  of  the  microcomputer,  which  was 
programmed  to  provide  a  display  similar  to  that  of  a 
spectrum  analyzer.  As  each  sample  was  received,  the  data 
point  was  plotted,  accumulating  vertically  with  an  abscissa 
representative  of  the  velocity.  The  resulting  plot 
represented  the  number  of  "hits"  at  a  particular  velocity 
over  the  total  sampling  time.  In  general,  the  accumulation 
of  100  or  more  points  of  data  provided  a  distribution  which 
was  Gaussian  in  nature,  and  which  could  then  be  saved  on 
flexible  disk  for  further  processing. 

The  ability  to  accurately  determine  mean  velocity  at 
selectable  points  within  the  jet  allowed  confirmation  of  the 
shape  of  the  velocity  profiles  within  the  body  of  the  jet. 
This  is  of  particular  value  as  most  models  so  far  postulated 
assume  a  Gaussian  profile  radially  and  that  the  centerline 
velocity  decays  inversely  with  the  streamwise  coordinate 
dimension.  It  was  expected  that  these  assumptions  would  be 
verified  during  the  course  of  the  investigation. 
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II. 


LASER  DOPPLER  VELOCIMETRY 


Laser  Doppler  velocimetry  (also  known  as  Laser  Doppler 
Anemometry),  is  a  relatively  new  technique  for  the 
measurement  of  velocities  in  fluids.  It  is  dependent  on  the 
wave  mechanics  of  light  in  the  medium  of  interest,  but 
requires  no  physical  probe  to  be  inserted  into  the  flow,  and 
is,  therefore,  non-intrusive  in  nature.  The  method  has  been 
used  to  measure  mixing  in  chemical  reactors,  flow  through 
nuclear  reactor  fuel  bundles,  chemically  reacting  flows,  and 
flows  in  other  environments  where  the  geometry  or  the  nature 
of  the  fluid  to  be  studied  precludes  the  use  of  other 
techniques.  It  has  been  used  in  both  gases  and  liquids. 

Drain  [3]  presents  the  theory  and  techniques  of  Laser 
Velocimetry  on  an  introductory  level  and  covers  the 
currently  used  methods  and  processes  in  some  depth.  Only 
those  topics  necessary  for  background  or  understanding  of 
the  methods  used  in  the  present  work  will  be  discussed 
here. 

A.  THE  REFERENCE  BEAM  TECHNIQUE  AND  THE  DOPPLER  SHIFT 

The  Doppler  shift  is  a  well  known  phenomenon  which  has 
been  widely  utilized  in  acoustics  and  in  radar 
applications.  Because  light  is  a  form  of  electromagnetic 
radiation,  scattered  light  from  particles  in  a  fluid  flow  is 
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also  Doppler  shifted  by  an  amount  determined  by  the  relative 
velocities  of  the  source,  the  particle,  and  the  detector. 
For  a  stationary  source  and  detector,  the  amount  of  the 
Doppler  shift  observed  is  a  function  of  the  half  angle  of 
scattering  and  the  component  in  the  direction  of  scattering 
of  the  velocity  vector  of  the  particle  scattering  the 
light.  As  the  geometry  of  the  system  is  a  known  quantity, 
the  velocity  of  the  particle  can  be  obtained  directly  from 
the  Doppler  shift  frequency.  This  is  the  principle  of 
operation  of  the  reference  beam  mode  in  Laser  Anemometry. 

Because  of  the  speed  at  which  light  propagates,  the 
derivation  of  the  Doppler  shift  must  include  relativistic 
effects  in  order  to  be  rigorous.  Applying  the  equations  of 
special  relativity  and  using  the  final  approximation  that 
the  velocity  of  the  particle  is  very  small  compared  to  that 
of  light,  the  Doppler  frequency  shift  can  be  shown  as: 

fj  ■  ( 2v/x )cos#sin(f/2) 

where  v  is  the  particle  velocity,  A  is  the  wavelength  of  the 
light,  and  4  is  the  angle  between  the  direction  of 
scattering  and  the  particle  velocity  vector,  and  #  is  the 
angle  between  the  original  beam  path  and  the  line  of  sight 
from  the  particle  to  the  detector.  The  scattering  direction 
is  defined  in  Mie  scattering  theory  as  the  direction  of  the 
bisector  of  the  obtuse  angle  formed  by  the  lines  from  the 
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source  to  the  particle  and  the  particle  to  the  detector. 
Figure  1  shows  the  geometry. 

In  practice  the  Doppler  shift  frequency  is  very  low  in 

comparison  to  the  frequency  of  the  laser  light  (32  MHz 

Doppler  shift  at  10  m/sec  for  red  light  such  as  from  a  He-Ne 

14 

laser  at  a  frequency  of  4.7  X  10  Hz,  a  ratio  of  7  parts  m 

g 

10  ) .  Such  a  low  variation  in  frequency  is  extremely 
difficult  to  detect.  Furthermore,  the  output  frequency 
response  of  most  detectors  is  far  lower  than  optical 
frequencies.  Accordingly,  a  technique  must  be  used  whereby 
improved  frequency  response  is  obtained  and  the  Doppler 
shift  more  readily  detected.  Mixing  Doppler  shifted  light 
with  light  of  the  original  frequency  is  a  technique  similar 
to  that  used  in  superheterodyne  radio  and  radar  transmitters 
to  move  signals  to  more  favorable  parts  of  the 
electromagnetic  spectrum.  The  light  of  the  original 
frequency  is  the  reference  beam  from  which  the  technique 
takes  its  name. 

When  two  coherent  (i.e.,  of  a  fixed  phase  relationship) 
signals  are  mixed,  the  result  is  a  signal  with  constituent 
waves  of  frequencies  corresponding  to  the  sum  and  the 
difference  of  the  original  frequencies.  In  optical  detector 
systems,  the  sum  component  is  at  a  frequency  much  too  high 
for  the  bandpass  characteristics  of  the  system  to  process 
and  is  thus  filtered  out.  The  difference  signal,  however. 
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is  at  a  lower  frequency  which  is  easily  processed,  and  which 
corresponds  to  the  Doppler  shift  frequency. 

Physically,  the  reference  beam  is  generally  split  off  of 
the  illuminating  beam  before  it  passes  through  the  fluid. 
It  is  redirected  to  be  finally  focused  on  the  detector 
pinhole,  mixing  with  the  scattered  light  from  the  main 
beam.  It  is  sometimes  necessary  to  attenuate  the  reference 
beam  in  order  to  assure  similar  intensities  of  the  signals 
at  the  surface  of  the  detector. 

The  reference  beam  technique  has  been  used  with  success 
in  a  variety  of  applications.  However,  its  applicability  is 
limited  by  a  low  signal-to-noise  ratio  which  renders  it  less 
useful  than  the  differential  Doppler  technique.  The 
arrangement  of  a  typical  reference  beam  LDA  is  found  in 
Figure  2. 

B.  THE  DIFFERENTIAL  DOPPLER  TECHNIQUE  AND  THE  FRINGE  MODEL 

If  two  beams  of  coherent  light,  both  of  equal 
intensities  are  crossed  in  the  body  of  a  moving  fluid,  light 
scattered  from  the  intersection  volume  will  be  shifted  by  an 
amount  equal  to  the  difference  between  the  Doppler  shifts  of 
the  light  scattered  from  the  individual  beams.  This  shift 
can  be  analyzed  similarly  to  the  case  of  the  reference  beam 
method,  but  a  simpler  model  is  available  to  explain  the 
results  obtained.  This  is  the  fringe  model. 
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The  formation  of  interference  fringes  has  been  observed 
using  diffraction  gratings,  pinholes  of  specified 
separation,  and  many  other  simple  optical  devices  [4].  These 
fringes  are  the  result  of  constructive  and  destructive 
interference  between  the  wave  fronts  of  the  light  where  they 
are  coincident.  The  same  result  occurs  when  beams  from 
coherent  sources  are  crossed  in  space.  The  maxima  of  the 
light  intensities  resulting  from  the  mixing  are  the  bright 
planes  of  the  fringe  pattern  and  correspond  to  the  loci  of 
points  at  which  the  wavefronts  are  in  phase.  If  the  beams 
are  crossed  at  angle  0,  the  spacing  of  the  fringes  can  be 
shown  to  be  given  bys 


d  =  X/2sin(  6/2) 

The  fringes  are  parallel  to  the  bisector  of  the  intersection 
angle  of  the  beams.  If  a  particle  passes  through  the 
pattern  at  velocity  v,  the  scattered  light  will  be  modulated 
in  intensity  at  a  frequency  given  by; 

fd  =  (2v/A )cos*sin(« /2) 


This  is  the  differential  Doppler  frequency,  and  the  process 
is  the  basis  of  the  differential  Doppler  method.  Here  *  is 
the  angle  of  the  velocity  relative  to  the  normal  to  the 


fringe  planes.  Figure  3  depicts  the  geometry  of  the  fringe 
model. 

The  frequency  of  light  scattered  from  the  intersection 
volume  using  the  differential  Doppler  technique  is 
independent  of  angle.  This  fact  allows  collection  of  the 
signal  at  any  angle  which  is  advantageous  to  the  purpose  at 
hand.  If  the  detector  is  positioned  in  the  hemisphere 
beyond  the  intersection  volume  from  the  laser  or  beam 
source,  the  term  "forward-scatter”  is  used.  Conversely,  if 
the  detector  is  on  the  same  side  as  the  source,  the  term 
"backscatter "  is  used.  Collection  in  the  forward-scatter 
mode  provides  greater  signal  strength  while,  if  backscatter 
is  used  at  180  degrees,  the  receiving  optics  can  be  combined 
with  the  transmitting  optics.  Backscatter,  however,  is 
limited  by  a  much  lower  signal  strength. 

Light  intensity  across  the  beam  fronts  is  not  uniform, 
but  varies  in  a  Gaussian  manner  radially.  As  a  consequence, 
the  volume  of  intersection  of  the  beams  is  ellipsoidal  in 
shape  and  the  intensity  of  the  fringes  varies  across  the 
cross-section  of  the  intersection  volume.  This  results  in  a 
lower  frequency  modulation  of  the  differential  Doppler 
signal  which  is  called  the  pedestal.  The  pedestal  is 
unwanted  noise  and  is  usually  removed  by  filtering  before 
signal  processing  begins. 

A  diagram  of  a  typical  differential  Doppler  arrangement 
is  found  in  Figure  4. 
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C .  FREQUENCY  SHIFTING 

As  a  particle  passes  through  the  intersection  volume  in 
the  differential  Doppler  mode,  the  light  scattered  is 
modulated  at  the  differential  Doppler  frequency 
corresponding  to  the  component  of  velocity  of  the  particle 
normal  to  the  fringe  planes .  For  flows  in  which  the 
direction  of  travel  is  always  known,  sufficient  information 
can  be  obtained  to  characterize  the  flow.  In  highly 
turbulent  flows  where  flow  reversals  are  present  or  in 
similar  cases  such  as  mixing,  an  ambiguity  exists  as  to  the 
direction  of  flow.  Furthermore,  since  the  signal  is 
generated  by  the  movement  of  particles  with  respect  -to  the 
fringes,  a  zero  velocity  results  in  an  absence  of  signal. 
Both  of  these  problems  are  eliminated  through  the  use  of 
frequency  shifting. 

The  technique  of  frequency  shifting  consists  of  altering 
the  frequency  of  one  or  both  of  the  crossed  beams  by  a  fixed 
amount.  This  can  be  accomplished  by  the  use  of  mechanical 
or  electro-optical  devices,  the  most  common  method  being  the 
use  of  the  electro-acoustic  Bragg  Cell. 

When  frequency  shifting  is  used,  the  fringes  will  move 
with  a  velocity  determined  by  the  frequency  shift. 
Directional  ambiguity  and  loss  of  signal  at  zero  velocity 
are  now  eliminated,  and  the  true  velocity  of  the  particles 
may  be  determined  by  subtracting  the  frequency  shift  from 
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the  measured  differential  Doppler  frequency,  and  then 
proceeding  with  the  computations  as  though  no  frequency 
shift  had  been  applied. 

The  Bragg  Cell  consists  of  a  block  of  acousto-optic 
material  on  one  side  of  which  a  piezoelectric  transducer  is 
affixed  to  impart  the  acoustic  excitation  to  the  block.  The 
acoustic  waves  traveling  through  the  cell  diffract  the  laser 
light  causing  the  it  to  be  shifted  in  frequency  by  an  amount 
equal  to  the  frequency  of  the  acoustic  wave.  Successive 
acoustic  waves  cause  the  signals  to  reinforce  yielding  an 
improved  efficiency  in  transmission.  It  has  been  found  that 
a  shift  frequency  of  about  40  MHz  can  provide  80%  efficiency 
and  result  in  sufficient  separation  of  beam  orders  to  allow 
for  easy  blocking  and  selection.  In  many  currently  used 
systems,  the  plus  first  order  beam  is  used  at  a  40  MHz  shift 
because  the  wavelength  of  the  acoustic  waves  is  comparable 
to  the  diameter  of  the  beams  to  be  shifted. 

Either  one  or  two  Bragg  Cells  may  be  used  in  frequency 
shifting  applications.  Usually  if  one  cell  is  used,  a 
downmix  frequency  of  32  to  35  MHz  is  used  to  reduce  the 
observed  signal  in  frequency  to  a  range  where  the  frequency 
response  characteristics  of  the  installed  instrumentation 
are  better  adapted  to  handle  it.  In  the  two  Bragg  Cell 
technique  one  is  usually  operated  at  the  most  efficient 
frequency  for  transmission  (40  MHz  is  typical)  and  the  other 
operated  at  the  amount  of  shift  desired  above  the 
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reference.  The  net  shift  detected  is  the  difference  of  the 
two,  providing  selectable  frequency  shifting  on  the  order  of 
magnitude  of  most  flow  measurement  velocities  encountered, 
while  maintaining  optimum  beam  power. 

D.  SIGNAL  DETECTION  AND  PROCESSING 

The  signal  of  interest  in  laser  velocimetry  is  optical 
in  nature  and  is  hence  detected  by  means  of  optical 
instrumentation.  A  photomultiplier  tube  or  photodiode  is 
most  often  used  as  the  primary  detector.  Photomultipliers 
provide  inherent  amplification  features  and  have  good 
frequency  response  characteristics  up  to  10  to  100  MHz.  They 
are  most  responsive  in  the  blue  to  green  portion  of  the 
spectrum  but  can  be  obtained  with  good  response  in  the  red 
to  infra-red  regions.  Photodiodes  usually  have  no  inherent 
amplification  attributes,  but  can  be  avalanched  to  provide 
some  multiplication.  They  are  usually  most  sensitive  in  the 
red  to  infra-red,  although  adequate  sensitivity  may  be 
obtained  across  the  visible  spectrum. 

Signal  analysis  can  be  accomplished  through  a  variety  of 
means.  The  method  used  is  dependent  on  the  type  of 
information  required.  For  low  Doppler  frequencies,  spectrum 
analyzers  can  be  used  to  provide  frequency  information  from 
the  signals.  Counter-tracker  arrangements  and 
autocorrelators  can  be  used  to  separate  the  noise  from  the 
signal  and  to  yield  the  frequency  information.  Any  signal 
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processing  arrangement  must  include  provision  for  filters  to 
adjust  the  bandpass  characteristics  of  the  system  to 
maximize  the  signal-to-noise  ratio.  Low  frequency  filters 
(high  pass  filters)  are  used  to  remove  the  pedestal  while 
low  pass  filters  are  adjusted  to  eliminate  high  frequency 
noise  such  as  that  from  thermal  sources  and  photomultiplier 
tube  "shot  noise". 

For  the  present  work,  a  counter  type  signal  processor 
was  used  to  provide  frequency  information  for  each  burst  as 
sampled.  The  processor  used  a  high  speed  clock  (250  MHz)  in 
conjunction  with  Schmitt  triggers  to  measure  the  duration  of 
a  burst  of  a  fixed  number  of  cycles.  The  frequency  was  then 
obtained  from  the  number  of  cycles  measured  and  the  time 
required  for  the  burst.  Filtering  was  provided  in  the 
apparatus  before  the  counter  tracker. 

E.  VELOCITY  AND  FRINGE  BIASING 

The  results  obtained  by  laser  anemometry  can  be  biased 
in  two  principal  ways.  The  first,  velocity  biasing,  is  the 
biasing  of  the  mean  velocity  by  the  disproportionate 
inclusion  of  high  frequency  bursts  over  low.  This  occurs 
primarily  at  low  data  rates  when  the  registration  of  a  data 
point  results  from  the  passage  of  a  single  particle  through 
the  intersection  volume.  In  a  flow  with  uniform  particle 
density,  more  particles  moving  at  a  high  flow  rate  will  pass 
through  the  measurement  volume  than  will  particles  moving  at 


24 


low  velocities .  This  is  a  consequence  of  a  high  velocity 
particle  having  a  higher  probability  of  encountering  the 
intersection  volume.  If  a  sufficiently  high  data  rate  is 
present  that  the  data  acquisition  system  samples  at  uniform 
intervals  at  a  rate  much  slower  than  the  data  rate,  velocity 
biasing  is  eliminated  and  the  time  average  frequency  will 
yield  the  correct  mean  velocity.  If  such  a  high  data  rate 
cannot  be  maintained,  the  data  must  be  weighted  with  the 
duration  of  each  burst  sampled.  This  will  correct  for  the 
low  probability  of  sampling  a  low  velocity  particle. 

The  second  biasing,  fringe  biasing,  is  a  result  of  the 
parallel  nature  of  the  fringes.  If  a  particle  passes  on  a 
course  parallel  to  the  fringes,  no  signal  will  be 
generated.  Therefore,  even  though  a  significant  number  of 
particles  having  a  zero  component  of  velocity  in  the 
direction  being  measured  can  exist,  only  particles  having  a 
non-zero  component  will  be  "seen".  The  use  of  frequency 
shifting  causes  the  fringes  themselves  to  move  and 
eliminates  the  bias  by  generating  a  relative  motion  between 
the  particles  and  the  fringes. 

F.  SEEDING  OF  FLOWS 

Many  flows  of  interest  to  the  LDV  user  do  not  include 
enough  particulate  matter  to  scatter  the  light  and  provide  a 
signal.  In  those  cases,  an  appropriate  seeding  particle 
must  be  chosen.  There  are  two  basic  requirements  for 
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scattering  particles.  First,  they  must  follow  the  flow 
sufficiently  well  to  adequately  represent  the  dynamics  of 
the  system.  Second,  they  must  scatter  sufficient  light  to 
provide  a  good  signal.  Consideration  must  also  be  given  to 
special  conditions  such  as  high  temperature  or  chemical 
reactivity  of  the  flows  in  the  selection  of  the  seeding 
particles.  The  manufacturers  of  most  commercial  laser 
Doppler  velocimeters  provide  guidance  in  the  selection  of 
appropriate  particles  for  use  with  their  apparatus. 


III.  CIRCULAR  BUOYANT  JETS 


The  characteristics  and  behavior  of  circular  buoyant 
jets  have  been  the  subject  of  many  studies  in  the  past  four 
decades.  Most  of  these  studies  have  involved  not  only 
experimental  observations  but  also  mathematical  model 
development  in  an  effort  to  accurately  predict  plume  and  jet 
behavior  under  a  variety  of  conditions.  By  far  the  most 
widely  utilized  model  is  the  entrainment  model  which 
attempts  to  predict  jet  behavior  based  on  integral  equations 
of  continuity,  momentum,  and  energy.  Other  approaches 
utilized  an  analysis  based  on  Prandtl's  mixing  length 
hypothesis  or  on  turbulence  kinetic  energy.  Schlichting  [5] 
discusses  these  more  traditional  approaches. 

The  entrainment  model  has  been  used  by  numerous  authors 
in  their  studies,  and  general  agreement  on  the  forms  of  the 
equations  and  constants  involved  has  been  reached .  Gebhart 
et  al.  [6j  have  provided  a  comprehensive  discussion  of  the 
various  models  and  have  employed  a  number  of  them  for 
comparative  calculations.  The  models  predict  not  only 
trajectory  and  velocity  behavior  in  jets  but  also  deal  with 
the  temperature  and  concentration  fields  found  in  jets  in 
uniform  and  stratified  ambients,  either  quiescent  or 
flowing.  For  the  present  work  only  the  information  which 
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pertains  to  velocity  distribution  and  jet  growth  will  be 
considered. 


A.  GENERAL  PROPERTIES  OF  JETS 

The  physical  extent  of  a  buoyant  jet  can  be  divided  into 
three  zones.  The  first,  the  Zone  of  Flow  Establishment,  is 
that  region  in  which  the  flow  transforms  from  the 
characteristics  found  in  the  nozzle  to  those  of  a  free  jet. 
In  this  region  momentum  effects  predominate  and  the 
velocity,  temperature,  and  concentration  profiles  found  in 
the  flow  field  are  typical  of  the  nozzle  itself.  The  flow 
at  the  core  of  the  jet  is  unaffected  by  the  surrounding 
ambient,  and  viscous  shear  effects  act  at  the  periphery  of 
the  jet  to  initiate  mixing  with  the  surrounding  fluid.  The 
Zone  of  Flow  Establishment  ends  where  the  turbulent  mixing 
reaches  the  core  of  the  jet. 

The  second  region  is  the  Zone  of  Established  Flow.  Here 
buoyant  effects  and  momentum  effects  of  the  jet  are  on  the 
same  order  of  magnitude.  Velocity  profiles  (as  well  as 
those  of  temperature  and  concentration)  are  similar 
throughout  the  region.  Interaction  with  the  ambient  remains 
through  viscous  shear  at  the  boundaries  and  mixing  proceeds 
throughout  the  jet.  The  initial  discharge  conditions  of  the 
jet  play  a  progressively  smaller  role  as  the  jet  grows 
toward  plume-like  behavior. 
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In  the  final  region,  the  far  field,  momentum  effects  of 
the  jet  are  negligible  and  the  behavior  is  influenced  by 
buoyancy  effects.  The  jet  is  free  to  be  convected  passively 
by  the  ambient.  Diffusion  acts  as  the  jet  loses  definition 
and  becomes  indistinguishable  from  the  surroundings. 

The  terms  jet  and  plume  refer  to  submerged  discharges 
which  have  primarily  momentum  and  buoyant  behavior 
respectively.  A  pure  jet  can  be  considered  to  be  a 
discharge  of  fluid  of  the  same  density  as  the  ambient 
injected  at  some  velocity  relative  to  the  surrounding 
fluid.  A  jet's  behavior  is  determined  by  its  momentum  when 
discharged  and  it  reacts  with  the  ambient  only  through 
viscous  shear  at  the  boundary.  Conversely,  a  pure  plume  can 
be  considered  to  be  a  discharge  of  fluid  into  an  ambient  of 
significantly  differing  density  injected  at  a  low  enough 
velocity  that  momentum  effects  are  negligible.  A  pure  plume 
may  also  be  created  by  the  injection  of  energy  to  the  fluid 
such  as  through  immersion  of  a  hot  wire  in  the  ambient.  A 
plume  may  be  either  positively  buoyant  or  negatively 
buoyant.  Clearly,  real  world  discharges  are  rarely  either 
pure  jets  or  pure  plumes.  In  this  work  the  terms  jet  and 
plume  will  be  used  interchangeably  and  will  be  understood  to 
apply  to  a  discharge  of  a  fluid  into  an  ambient  which  will 
be  considered  to  be  infinite  in  extent  and  which  is  not 
required  to  have  the  same  composition  or  properties  as  the 
jet. 


29 


Ambients  may  be  described  as  stratified  or  unstratified 
as  well  as  flowing  or  quiescent.  A  stratified  environment 
is  one  which  possesses  layers  of  fluid  of  differing 
properties  within  its  extent.  Situations  of  interest 
because  of  parallels  in  nature  are  those  involving  vertical 
stratification  of  water  either  thermally  or  with  respect  to 
salinity.  In  the  laboratory,  the  environment  stratified 
with  respect  to  salinity  is  the  easiest  to  control  and  is, 
therefore,  the  one  usually  studied.  Here,  however,  the 
object  is  to  compare  LDV  measurements  with  theory. 
Accordingly,  the  ambient  of  interest  will  be  selected  as 
quiescent,  unstratified  water  in  the  interest  of 
simplicity. 

In  addition  to  classification  as  jets  or  plumes, 
discharges  may  be  described  by  their  angle  of  discharge  with 
respect  to  the  gravity  field.  A  vertical  jet  is  one 
discharged  parallel  to  the  gravity  field  (in  either  the  same 
or  the  opposite  sense)  and  a  horizontal  jet  is  one 
discharged  normal  to  it.  Obviously  a  jet  can  assume  any 
orientation,  and  the  orientation  may  change  over  the  flow 
field  if  buoyancy  is  present.  For  the  present  work,  the 
vertical  jet  has  been  selected  to  provide  for  a  stable, 
predictable  trajectory  in  the  presence  of  buoyancy  effects. 
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B.  ENTRAINMENT  MODEL  FOR  THE  QUIESCENT  AMBIENT 

Morton  et  al.  [7]  were  the  first  to  use  the  entrainment 
model  which  supposes  that  the  rate  at  which  fluid  is  brought 
into  the  jet  is  proportional  to  the  centerline  velocity  of 
the  jet  and  to  the  effective  circumferential  area  of  the 
jet.  Mathematically, 


E  oc  2frBU 

m 

Here  E  represents  the  volumetric  rate  of  entrainment  into 

the  jet,  B  is  some  characteristic  dimension  of  the  jet  yet 

totbe  defined,  and  U  is  the  mean  centerline  velocity  of  the 

m 

jet.  This  rate  of  entrainment  must  be  equal  to  the  rate  of 
change  of  volumetric  flow  within  the  jet  or, 

dQ/dS  =  E 


Here  Q  is  the  volumetric  flow  rate  and  S  is  the  streamwise 
coordinate  of  the  jet.  The  constant  of  proportionality 
defining  E  is  called  the  entrainment  constant  or  entrainment 
coefficient,  a.  Utilizing  the  notation  thus  defined,  the 
entrainment  function  becomes 

E  =  2na BU 

in 
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The  assumptions  used  in  the  solution  of  the  equations 
for  differential  modeling  of  the  plumes  are  as  follows: 

1.  The  turbulent  jet  flow  is  steady. 

2.  Since  the  jet  flow  is  fully  turbulent,  radial 
molecular  diffusion  is  neglected,  compared  to  radial 
turbulent  transport. 

3.  Streamwise  turbulent  transport  is  negligible  when 
compared  with  the  convective  streamwise  transport. 

4.  The  variation  of  the  fluid  density  throughout  the  flow 
field  is  small  compared  to  a  chosen  reference 
density.  Density  variation  is  significant  only  in 
buoyancy  terms.  (Boussinesq  approximation). 

5.  There  is  no  variation  of  other  fluid  properties 
throughout  the  flow  field. 

6.  Pressure  is  hydrostatic  throughout  the  flow  field. 

7.  The  jet  is  axisymmetric .  That  is,  throughout  the  flow 
field  there  is  no  variation  circumferentially  of  any 
important  parameter . 

Most  modelers  have  assumed  a  Gaussian  profile  for 
velocity  in  the  jet.  That  is,  the  velocity  in  the  jet  cross 
section  will  vary  as: 

U  =  U  exp(-r2/B2) 
m 

2 

Here  B  is  twice  the  variance  from  the  mean  centerline 
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velocity  of  the  local  velocity  across  the  jet.  Using  this 
definition,  B  can  be  seen  to  represent  the  radial  coordinate 
at  which  the  mean  velocity  has  decayed  to  1/e  of  the  mean 
centerline  velocity.  B  is  therefore  representative  of  the 
nominal  half  width  of  the  jet  and  is  the  quantity  used  in 
the  entrainment  relations.  A  similar  Gaussian  profile  has 
been  assumed  for  variations  of  temperature  and  concentration 
to  complete  the  models,  but  a  spreading  factor  equal  to  the 
square  of  the  turbulent  Schmidt  number  is  added  to  the 
denominator  in  the  exponent.  See  Gebhart  [6]  for  details. 

In  the  Zone  of  Established  Flow,  which  is  of  principal 
interest  here,  buoyant  forces  and  momentum  effects  are  both 
important.  The  relative  importance  of  these  quantities  can 
be  represented  by  the  value  of  the  densimetric  Froude 
number,  F.  It  is  defined  by: 

F  =  U0/tgD(,a-,0)/,0)1/2 


The  importance  of  the  initial  momentum  is  given  by  the 
discharge  velocity  UQ  and  the  contribution  of  the  buoyancy 
force  is  represented  by  the  density  difference.  It  is  a 
measure  of  the  velocity  generated  by  the  buoyant  force. 
Hereafter,  the  densimetric  Froude  number  will  be  referred  to 
as  simply  the  Froude  number.  In  some  models,  <  ie  local 
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Froude  number  is  used.  It  is  defined  as: 


PL  ' 

Mathematical  models  using  the  entrainment  function 
provide  prediction  of  velocity,  concentration  and 
temperature  profiles,  trajectory,  and  jet  width  once  the 
entrainment  constant  is  specified.  The  determination  of  the 
entrainment  constant  has  been  empirical  in  nature,  with 
several  researchers  recommending  specific  values  or  forms 
for  correlations.  While  many  of  these  have  been  formulated 
for  variable  orientation  or  for  flowing  ambients,  they  will 
not  be  addressed  here;  rather  those  relations  germane  to 
quiescent  ambients  and  vertical  jets  will  be  discussed. 

Albertson  et  al .  [83  conducted  measurements  in  order  to 
verify  mathematical  models  and  to  determine  the  value  of  the 
entrainment  constant.  They  found  that  for  non-buoyant  jets 
(F  =«>),  the  appropriate  value  of  a  was  0.057.  Other 
researchers  have  been  in  basic  agreement  with  that  figure. 
Albertson  also  found  that  the  centerline  velocity  decayed  as 
1/S  (that  is,  inversely  as  the  streamwise  coordinate)  for 
pure  momentum  jets  and  that  the  half  width  increased 
directly  as  the  streamwise  coordinate.  These  dependencies 
are  consistent  with  the  derivations  involving  mixing  length 
hypotheses  found  in  Schlichting  [53*  The  assumption  of  the 


Gaussian  velocity  profile  in  the  Zone  of  Established  Flow 
was  also  confirmed. 


For  pure  buoyant  plumes  Abraham  [9]  proposed  a  value  of 
0.085.  List  and  Imberger  [103  and  Fan  [11]  recommended  0.082 
for  the  value  of  the  entrainment  constant  for  pure  buoyant 
plumes.  Fan  also  recommended  0.057  for  momentum  jets. 

In  an  effort  to  model  jets  which  were  neither  pure 
momentum  jets  or  buoyant  plumes,  Morton  et  al.  [7]  proposed 
that  the  entrainment  constant  be  given  by: 

a  -  0.057  +  aj/F^ 

Hirst  [12]  proposed  a  value  of  0.97  for  a£  and  that  the  sine 
of  the  angle  of  inclination  of  the  jet  to  the  horizontal 
should  be  included  on  the  correction  term  also. 

Davis  et  al .  [13]  suggest  the  entrainment  coefficient 
be  given  by: 


«  =  0.057  +  0.083/F0*3 

They  also  measured  the  decay  of  centerline  jet  velocity  with 
streamwise  coordinate  as  well  as  the  growth  of  the  jet 
dimensions.  Their  apparatus  discharged  a  saline  jet 
downward  into  fresh  water,  and  they  used  hot  film  anemometry 
for  measuring  velocity.  No  firm  correlation  was  reached 
with  respect  to  jet  growth  because  of  the  scatter  of  the 
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data.  Their  velocity  decay  data  showed  a  decrease  in  the 

rate  of  decay  with  decreasing  Froude  number.  Th-is  is 

consistent  with  the  data  from  Albertson's  studies  for 

infinite  Froude  number  and  that  of  Rouse  et  al.  [14]  who 

reported  that  the  centerline  velocity  for  a  pure  buoyant 

-1/3  .  .  . 

plume  decayed  as  S  '  .  Davis  also  verified  the  Gaussian 

nature  of  the  velocity  distribution  in  the  Zone  of 
Established  Flow. 

C.  SUMMARY  OF  PREDICTED  BEHAVIOR 

Based  on  the  foregoing,  the  jet  behavior  observed 

through  Laser  Velocimetry  is  expected  to  be  predictable. 

Specifically,  the  jet  is  expected  to  grow  in  the  half  width, 

B,  at  a  rate  proportional  to  the  streamwise  coordinate,  S. 

The  centerline  velocity  is  expected  to  decay  at  a  rate 

—  1/3 

proportional  to  1/S  for  the  pure  momentum  case  and  to  S 

for  the  pure  buoyant  case.  By  fitting  the  observed  jet 

velocity  field  to  appropriate  power-law  relations,  an 

expression  for  the  volumetric  flow  rate  Q  as  a  function  of 

streamwise  coordinate  can  be  obtained.  This  expression  can 

then  be  differentiated  with  respect  to  S  to  obtain  the 

entrainment  rate.  As  B  and  U  are  known,  the  value  of  the 

m 

entrainment  coefficient  at  any  level  can  be  calculated.  A 
mean  entrainment  coefficient  can  then  be  obtained  for  each 
jet  and  compared  with  theory.  Values  of  mean  entrainment 
coefficient  are  expected  to  vary  between  0.057  and  0.082. 


IV.  APPARATUS  AND  PROCEDURES 


The  present  work  represents  the  first  in  a  series  of 
studies  investigating  the  characteristics  of  turbulent 
buoyant  jets  using  Laser  Velocimetry.  As  such  the  present 
research  concentrated  on  the  development  and  validation  of 
the  technique.  This  leads  to  a  developmental  approach  to 
the  use  of  the  LDV  and  to  the  techniques  of  measurement 
employed.  Also,  the  LDV  progressed  through  several 
evolutionary  stages  during  the  course  of  the  work  as  the 
methods  used  were  refined  and  problems  with  original 
approaches  were  identified  and  corrected. 

The  information  presented  in  this  chapter  is  the  result 
of  a  number  of  changes  to  the  original  apparatus  and 
approaches  to  the  employment  of  the  LDV.  It  reflects  the 
form  of  the  experimental  arrangement  and  the  method  of  use 
which  proved  to  be  most  effective  in  measuring  the  fluid 
velocities . 

A.  LASER  DOPPLER  VELOCIMETER  OPTICS 

The  LDV  used  in  the  present  work  was  a  TSI,  Incorporated 
model  9100-5  Laser  Doppler  Velocimeter.  It  is  a  modular 
instrument  which  is  intended  to  be  user  assembled  with 
components  included  as  required  for  the  application.  All 
components  are  mounted  on  an  anodized  aluminum  base  which  is 
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equipped  with  mounting  brackets  for  securing  the  modules  to 
the  base.  The  optical  modules  are  secured  to  each  other  by 
the  means  of  thumb  screws  and  the  entire  assembly  is  mounted 
on  the  base  on  rotating  mounts  which  permit  the  choice  of 
angle  of  velocity  component  to  be  measured.  Figure  5 
presents  a  schematic  representation  of  the  LDV  optical 
assembly. 

The  laser  unit  is  mounted  in  the  lower  part  of  the 
extruded  base.  It  is  a  Spectra  Physics  15  milliwatt 
Helium-Neon  laser  which  operates  in  the  continuous  wave  (CW) 
mode.  Power  to  the  laser  is  furnished  through  a  separate 
power  supply.  The  LDV  is  assembled  with  the  laser  unit 
pointing  away  from  the  test  section.  The  laser  beam  is 
focused  onto  a  front-surface  mirror  located  in  a  housing  at 
the  end  of  the  base  which  reflects  it  to  a  similar  mirror  at 
the  top  of  the  housing  from  which  it  is  reflected  in  the 
direction  of  the  test  section.  Both  of  the  mirrors  are 
provided  with  adjusting  knobs  which  act  through  wedge 
assemblies  to  permit  alignment  of  the  beam. 

From  the  top  front-surface  mirror,  the  beam  is  focused 
on  an  aperture  located  in  the  center  of  the  after  rotating 
mount  for  the  optics  assembly.  It  then  passes  to  a 
prismatic  beam-splitter  module  which  separates  the  single 
beam  into  two  coherent  beams  of  equal  intensity  which  travel 
in  parallel  paths  down  the  base. 
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The  two  parallel  beams  then  enter  the  frequency  shifting 
module.  It  contains  two  Bragg  Cells,  one  in  the  path  of 
each  of  the  beams.  The  Bragg  Cells  are  acoustically 
vibrated  by  electrical  signals  provided  from  a  frequency 
shift  unit  and  power  amplifier  located  remotely  from  the 
base.  One  cell  is  driven  at  40  MHz  and  the  other  at  a 
frequency  above  the  first  by  an  amount  relectable  by 
pushbuttons  on  the  frequency  shift  unit.  Each  Bragg  Cell  is 
followed  by  a  beam  steering  lens  eccentrically  mounted  in 
wedge  assemblies.  These  lenses  can  be  rotated  to  realign 
the  beams ,  to  correct  for  the  deviation  introduced  by  the 
frequency  shifting  process.  In  addition,  the  next  module, 
the  beam  steering  module  contains  another  such  assembly  in 
the  path  of  one  beam.  This  is  provided  to  assist  in  beam 
crossing  during  final  alignment  of  the  device. 

The  next  module,  the  receiving  optics  unit,  contains 
another  rotating  mount  assembly.  As  the  beams  proceed 
towards  the  test  section,  they  travel  through  passages  in 
the  outer  portion  of  the  assembly. 

Following  the  receiving  optics  assembly  a  beam  expander 
is  installed  if  the  beam  separation  to  be  used  was  50  mm. 
If  the  beam  expander  module  is  omitted,  the  beam  separation 
is  22  mm.  The  latter  was  selected  for  the  present  work 
because  it  provides  a  superior  signal-to-noise  ratio  and  a 
smaller  intersection  volume  with  the  selected  transmitting 
lens  than  does  the  50  mm  separation.  This  module  is 


followed  by  a  beam-stop  assembly,  a  third  and  final  rotating 
mount  and  the  transmitting/receiving  lens. 

The  purpose  of  the  beam-stop  assembly  is  to  block  out 
unwanted  orders  of  the  diffraction  pattern  created  by  the 
Bragg  Cells.  The  order  used  is  the  plus  first  order  beam,  in 
which  the  proper  difference  in  frequency  exists.  The 
unshifted  components  (zeroth  order)  and  those  of  higher 
order  must  be  blocked  in  order  for  a  proper  Doppler  shifted 
signal  to  be  provided. 

The  lens  used  for  this  study  has  a  focal  length  of  120 
mm.  Other  sizes  are  available  from  TSI,  extending  in  focal 
length  up  to  600  mm  for  the  model  9100-5.  In  the  present 
work,  it  was  found  that  the  increase  in  the  length  of  the 
intersection  volume  and  the  lowering  of  signal-to-noise 
ratio  which  result  from  the  use  of  the  longer  focal  length 
lens  or  greater  beam  separation  degraded  the  quality  of 
measurement.  Accordingly,  the  final  configuration  selected 
was  the  120  mm  focal  length  lens  with  the  22  mm  beam 
spacing.  Since  180  degree  backscatter  is  the  intended  mode 
for  this  device,  the  transmitting  lens  is  also  the  receiving 
lens . 

The  focusing  action  of  the  transmitting  lens  causes  the 
beams  to  cross,  forming  an  intersection  volume.  Light 
scattered  from  the  intersection  volume  by  particles  flowing 
in  the  fluid  is  gathered  by  the  same  lens  and  passed  back 
down  the  axis  of  the  LDV.  It  is  focused  by  a  lens  in  the 


receiving  optics  assembly  on  an  aperture,  then  refocused  by 
another  lens  to  the  photomultiplier  tube  by  way  -of  a 
mirror.  The  aperture  is  located  in  a  holder  assembly  and 
can  be  changed  providing  larger  or  smaller  size  holes. 

B.  THE  LASER  DOPPLER  VELOCIMETER  ELECTRONICS 

The  optics  assemblies,  when  properly  aligned,  cause  the 
image  of  the  intersection  volume  of  the  two  laser  beams  to 
be  focused  on  the  detector  surface  of  a  photomultiplier 
tube.  This  tube  is  powered  by  an  external  power  supply 
equipped  with  an  adjustable  gain  control.  It  is  also 
provided  with  a  self  protection  feature  to  prevent  "cooking" 
the  tube  should  the  light  intensities  encountered  be  greater 
than  appropriate  for  the  selected  gain.  The  photomultiplier 
tube  converts  the  variations  in  the  intensity  of  the  light 
signal  emanating  from  the  intersection  volume  into  a  voltage 
signal  varying  at  the  same  frequency.  This  signal  is  the 
raw  Doppler  signal  from  which  the  velocity  can  be 
determined.  Figure  6  shows  typical  Doppler  bursts  both  with 
and  without  the  pedestal  signal  component. 

The  Doppler  burst  is  caused  by  the  variation  in 
intensity  of  the  light  scattered  as  a  particle  passes 
through  the  laser  beam  intersection  volume.  Passage  of  a 
single  particle  through  the  ellipsoidal  volume  illuminated 
by  the  fringe  pattern  causes  light  to  scatter,  varying  in 
intensity  at  a  rate  corresponding  to  the  rate  at  which  the 
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bright  fringes  were  encountered.  The  frequency  of  the 
optical  signal  thus  created  is  the  differential  Doppler 
frequency,  a  knowledge  of  which  permits  computation  of  the 
particle  velocity.  The  envelope  of  the  Doppler  signal  is 
the  total  length  of  time  occupied  by  one  burst  as  shown  in 
Figure  6.  A  burst  is  defined  as  the  signal  generated  by  the 
passage  of  one  particle  through  the  laser  beam  intersection 
volume . 

The  photodetector  output  is  sent  to  the  Counter  for 
processing.  It  consists  of  three  distinct  sections;  the 
input  conditioner,  the  timer,  and  the  display  section.  For 
the  present  work,  only  the  first  two  sections  were  directly 
employed . 

The  input  conditioner's  function  is  to  take  the  raw 
signal  from  the  photomultiplier  tube  and  generate  either  one 
or  two  measurable  envelopes  for  the  timer  to  use,  depending 
on  the  mode  of  operation  selected.  The  N-cycle  mode  was 
used  in  the  present  work.  In  this  mode,  the  input 
conditioner  generates  a  signal  the  envelope  of  which  is  of 
the  saune  duration  as  the  Doppler  burst  consisting  of  N 
cycles.  Here  N  refers  to  the  number  of  cycles  which  has 
been  selected  on  the  front  panel  to  be  measured  in  a  burst 
from  one  particle  passing  through  the  intersection  volume. 
A  second  envelope  is  generated  of  time  equal  to  that  for  N/2 
cycles  to  pass.  The  N/2  cycle  envelope  is  used  to  validate 
the  data  during  processing.  These  envelopes  are  generated 
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through  the  use  of  Schmitt  triggers  and  threshold  detecting 
circuitry.  A  full  discussion  of  the  electronics  involved  is 
beyond  the  scope  of  this  work.  For  further  information 
refer  to  the  manufacturer’s  technical  documentation  for  the 
LDV  [15] . 

In  the  total  burst  mode,  which  was  not  used  in  this 
study,  the  input  conditioner  generates  one  envelope, 
corresponding  to  the  time  of  the  total  burst.  Circuitry  is 
included  to  allow  counting  the  number  of  cycles  in  the 
burst.  This  mode  is  useful  if  the  data  rate  is  low  and 
correction  for  velocity  biasing  is  required. 

The  input  conditioner  also  contains  filtering  circuitry 
for  the  removal  of  the  pedestal  component  of  the  signal  (the 
low  limit  filter)  and  for  the  isolation  of  noise  from  the 
signal  of  interest  (the  high  filter).  These  are  adjusted  in 
combination  to  provide  a  signal  which  is  free  of  spurious 
noise  while  still  allowing  enough  of  the  signal  to  pass  to 
adequately  measure  the  flow  fluctuations.  Additionally,  an 
amplitude  limit  control  is  provided.  Its  function  is  to 
reject  signals  of  greater  than  a  specified  amplitude.  This 
has  the  effect  of  eliminating  signals  from  larger  particles 
which  may  not  be  adequately  following  the  flow. 

The  final  control  on  the  input  conditioner  is  a  selector 
which  provides  control  of  the  number  of  cycles  used  in  the 
N-cycle  mode  and  the  minimum  cycles  required  for  a  valid 
burst  in  the  total  burst  mode.  Additional  fittings  include 
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jacks  for  the  photomultiplier  output  signal  and  filtered 
output  monitor  for  use  on  an  oscilloscope. 

The  timer  utilizes  a  high  speed  digital  clock  operating 
at  a  frequency  of  250  MHz  and  with  a  temporal  resolution  of 
1  nsec.  When  the  envelope  signal  from  the  input  conditioner 
is  gated  into  the  clock,  the  time  for  the  envelopes  to  pass 
is  measured.  In  the  N-cycle  mode,  the  time  for  N  cycles  was 
compared  to  the  time  for  N/2  cycles  via  a  direct  division. 
If  the  deviation  from  double  the  N/2  time  exceeds  the 
percentage  selected  by  switch  on  the  front  panel  the  burst 
is  discarded.  This  comparison  is  not  used  in  the  total 
burst  mode. 

On  the  front  panel  of  the  timer  section  output  jacks  for 
the  analog  output  and  direct  digital  signal  output  are 
found.  Additionally,  selector  buttons  and  lights  indicating 
active  values  are  found  for  the  exponent  to  be  applied  to 
the  time  signal  in  the  interpretation  formula.  The 
selection  of  the  exponent  value  can  be  either  manual  or 
automatic;  however,  when  autoranging  is  selected,  a  stable 
signal  frequently  can  not  be  obtained.  Accordingly,  in  the 
N-cycle  mode  in  particular,  the  LDV  should  be  operated  in 
the  manual  mode  once  the  magnitude  of  the  exponent  for  the 
current  flow  level  has  been  experimentally  determined. 
Autoranging  should  be  used  in  the  total  burst  mode  because 
of  the  wide  range  of  the  numbers  of  cycles  and  the  time  for 
the  burst  which  will  be  detected. 
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The  front  panel  of  the  timer  also  contains  a  standard 
ICC/Cannon  37  pin  connector  for  transmission  of  control  and 
readout  information  directly  to  a  digital  computer.  The 
signals  which  are  available  via  this  connection  are  as 
follows: 

1.  A  16  bit  word  consisting  of  a  12  bit  mantissa  and  a  4 
bit  exponent  representative  of  time  for  the  burst. 

2.  An  8  bit  word  for  the  representation  of  the  number  of 
cycles  in  the  burst.  This  would  be  useful  in  the 
total  burst  mode. 

3.  A  one  bit  data  ready  signal  for  computer  handshaking. 

4.  A  one  bit  data  hold  signal. 

5.  A  one  bit  single  measurement  per  burst  signal.  This 
can  be  used  for  control  of  some  LDV  functions  from  a 
remote  location. 

C.  THE  COMPUTER  AND  SUPPORTING  ELECTRONICS 

The  computer  used  to  control  data  acquisition  is  a 
Hewlett-Packard  9826  microcomputer.  It  uses  a  16-bit  word 
length  and  is  capable  of  communicating  directly  with  the  LDV 
via  an  interface  bus  arrangement.  Because  of  the  16  bit 
word  length,  the  data  transfer  arrangement  used  in  the 
present  work  is  only  appropriate  for  the  N-cycle  mode.  This 
is  because  there  is  no  arrangement  to  hold  data  from  the 
digital  word  representing  the  number  of  cycles  for  later 
input  to  the  computer.  Using  the  N-cycle  mode,  the  number 
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of  cycles  per  burst  to  be  measured  is  entered  by  the 
operator  at  the  computer  console.  This  particular  HP9826  is 
fitted  with  a  standard  RS-232C  interface  and  acquired  data 
is  transmitted  via  modem  to  a  larger  computer  for 
reduction.  In  the  interim,  the  experimental  information  is 
stored  on  flexible  disks . 

For  any  work  involving  laser  velocimetry,  it  is 
essential  that  the  operator  be  able  to  monitor  the  quality 
of  the  signal  being  processed.  For  this  purpose  an 
oscilloscope  is  irreplaceable.  In  the  present  work  a 
recording  oscilloscope  capable  of  dual  trace  operations  and 
with  an  operating  range  of  up  to  10  MHz  was  used.  In 
general  for  use  in  LDV  applications,  a  frequency  capability 
of  up  to  100  MHz  is  desirable. 

A  frequency  generator  is  required  for  setup  and  checkout 
of  the  LDV  electronics.  It  should  be  capable  of  up  to  1  MHz 
output  and  was  sufficiently  precise  that  a  stable  velocity 
reading  can  be  generated. 

D.  THE  SIGNAL  PATH 

When  a  particle  passes  through  the  intersection  volume 
of  a  differential  Doppler  LDV,  such  as  this  one,  the  light 
scattered  in  all  directions  is  modulated  at  the  differential 
Doppler  frequency.  A  lower  frequency  modulation  due  to  the 
variation  in  light  intensity  across  the  beams  is  also 
present  and  constitutes  the  pedestal.  The  light  scattered 


from  the  intersection  volume  is  focused  onto  the 
photomultiplier  tube  by  the  receiving  optics.  The 
photomultiplier  tube  converts  this  optical  signal  to  an 
electrical  signal  with  a  voltage  proportional  to  the 
intensity  of  the  light.  This  photomultiplier  tube  output  is 
passed  to  the  input  conditioner. 

On  arrival  at  the  input  conditioner,  the  signal  is  first 
checked  by  the  amplitude  limiting  circuitry  which  eliminates 
any  bursts  of  excessive  amplitude.  For  the  purposes  of  the 
current  work,  the  amplitude  limit  was  not  set  because  the 
quality  of  the  data  gathered  indicated  that  no  particles 
sufficiently  large  to  affect  accuracy  were  being  processed. 

Next  the  signal  is  filtered,  first  by  the  high  pass 
filter  which  is  set  to  eliminate  the  low  frequency  component 
of  the  signal,  the  pedestal.  The  low  pass  filter  next 
removes  components  of  frequency  higher  than  its  setting, 
which  eliminates  shot  noise  in  the  photomultiplier  tube  and 
thermal  source  interference.  The  signal  is  then  ready  for 
processing  by  the  Schmitt  trigger  system  for  generation  of 
the  envelopes. 

This  new  signal  is  passed  to  the  timer  where  the  leading 
and  trailing  edges  of  the  envelope  are  gated  by  the  timer 
and  the  signal  representing  the  time  of  the  burst  is 
generated.  From  here,  the  output  is  passed  to  the  display 
section  and  to  the  output  jacks  for  the  digital  or  analog 
output.  In  the  present  work  the  signal  was  read  off  of  the 
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37  pin  connector  in  the  front  panel  and  processed  by  the 
HP9826. 

The  signal  received  by  the  HP9826  is  not  directly 
interpretable  as  the  time,  but  yields  the  Doppler  frequency 
when  used  in  the  following  formula: 

f ,  =  (N  X  109)/(Dm  X  2n”3) 
a  m 

Here  D  is  the  value  of  the  mantissa  of  the  16  bit  computer 
m 

word,  fd  is  the  differential  Doppler  frequency,  n  is  the 
exponent  part  of  the  16  bit  computer  word,  and  N  is  the 
number  of  cycles  in  the  burst,  which  is  entered  by  the 
operator  on  program  initialization  and  set  on  the  front 
panel  of  the  counter. 

The  computer  program  under  which  the  process  is 
controlled  applies  the  foregoing  formula  and  computes  the 
velocity  represented  by  the  sample  using  the  geometry  of  the 
LDV  and  the  following  formula: 

v  =  xf/{2sin(#/2) } 

Position  and  orientation  data,  as  well  as  temperature  of  jet 
and  ambient  fluid,  are  entered  manually  during  data 
acquisition.  All  quantities  are  stored  on  disk  in  a  logical 
format  which  is  directly  transferred  to  the  VAX  mainframe 
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computer  via  the  telephone  line.  It  is  processed  there 
using  statistical  methods  to  extract  the  mean  velocity,  jet 
half  width,  and  entrainment  coefficients.  If  desired,  the 
turbulence  intensities  and  normal  component  of  velocity  can 
be  measured  and  processed  using  the  computer  programs 
developed  in  the  present  work. 

The  flow  of  the  signal  is  schematically  represented  in 
Figure  7. 

E.  THE  FLUID  LOOP 

The  jet  studied  for  the  present  work  was  generated  in  a 
tank  measuring  0.625  X  0.333  X  0.476  meters.  The  depth  of 
the  water  in  the  tank  was  0.419  meters.  The  walls  of  the 
tank  were  constructed  of  0.00635  meters  thick  glass.  Flow 
was  provided  by  a  centrifugal  pump  via  a  control  valve,  a 
heater  and  calibrated  rotameter  to  a  nozzle  0.003175  meters 
in  diameter.  The  nozzle  was  located  in  the  tank  0.067 
meters  from  the  nearest  wall,  a  distance  of  21  diameters. 

Temperature  was  monitored  by  two  thermocouples,  the 
first  embedded  in  the  nozzle  block  for  the  measurement  of 
the  temperature  of  the  effluent.  The  other  was  located  in 
the  bulk  of  the  fluid  in  the  tank  for  the  measurement  of  the 
ambient  temperature. 


Early  experiments  revealed  that  significant  amounts  of 
heat  accumulated  in  the  tank  from  the  injection  of  the 
warmer  water  in  the  jet,  with  the  result  that  the 


temperature  could  not  be  held  constant  during  the  course  of 
an  experiment.  The  situation  was  corrected  by  the 
installation  in  the  tank  of  a  cooling  coil  connected  to  a 
constant  temperature  bath.  Subsequent  to  the  installation 
of  the  coils,  the  temperature  remained  constant  in  the  tank 
over  the  course  of  the  experiments  to  within  1  degree 
Fahrenheit.  It  should  be  noted  that  the  presence  of  the  coil 
created  temperature  gradients  in  the  tank,  and  a  better 
system  of  temperature  control  should  be  found  for  use  in 
experiments  which  may  require  more  precise  information  on 
the  ambient  temperature  fields  than  was  required  in  the 
present  work. 

Early  experiments  also  indicated  that  significant 
circulatory  effects  were  generated  in  the  tank  when  flow 
rates  in  excess  of  40  percent  of  full  scale  on  the  rotameter 
were  used.  Accordingly,  flow  rates  below  that  point  were 
maintained  for  all  later  work. 

Figure  8  is  a  schematic  of  the  fluid  loop  used  in  this 
work. 

F.  POSITIONING  APPARATUS 

Positional  control  and  support  for  the  LDV  were  provided 
by  a  milling  machine  table  from  which  the  machine  head  had 
been  removed.  This  arrangement  allowed  the  LDV  to  be 
positioned  reproduceably  to  within  0.0000254  m,  (.001 
inches),  spatially.  It  was  also  sufficiently  massive  that 
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vibration  was  minor,  and  could  therefore  be  discounted  as  a 
source  of  significant  error. 

While  the  movement  of  the  milling  machine  provided  good 
control,  some  hysteresis  was  present  in  the  left-right,  and 
fore-aft  movements.  This  effect  was  minimized  by  always 
approaching  desired  position  in  the  same  direction.  Checks 
of  the  method  resulted  in  good  accuracy  being  demonstrated. 

G.  DATA  ACQUISITION  PROCEDURES 

Flow  rates  and  trace  locations  for  the  data  acquisition 
were  selected  to  provide  similar  conditions  to  previous 
studies  in  order  to  compare  results  and  validate  the 
technique.  While  an  infinite  Froude  number  initially 
appeared  to  be  the  best  choice,  it  was  unobtainable  with  the 
experimental  apparatus.  Relatively  high  Froude  numbers,  on 
the  order  200  to  500  were  easily  attainable  and  were  hence 
used . 

The  height  of  the  first  trace  was  selected  to  remain 
beyond  the  accepted  6.2  diameters  which  delineates  the  Zone 
of  Flow  Establishment  (ZFE)  from  the  Zone  of  Established 
Flow  (ZEF)  in  order  to  keep  the  expected  velocity  profiles 
as  nearly  Gaussian  as  possible.  Velocity  traces  were 
obtained  at  regular  intervals  to  ease  positioning.  The 
number  of  points  in  the  trace  was  selected  to  provide  a  wide 
base  on  which  the  curve  fit  routines  could  operate,  as  was 
the  number  of  traces  made  in  a  given  experiment. 
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For  seeding  of  the  flow  field,  several  substances  were 
tried  without  success  before  powdered  alumina  was  found  to 
provide  very  good  reflectivity,  uniform  (relatively) 
particle  size,  and  good  flow-following  characteristics. 
Also  of  concern  was  the  abrasive  characteristics  of  a 
seeding  particle  in  the  centrifugal  pump.  The  seeding 
particle  recommended  by  TSI  is  silicon  carbide,  but  its 
specific  gravity  is  greater  than  that  of  alumina  and  its 
increased  abrasive  qualities  were  believed  to  be  potentially 
prohibitive,  as  was  its  cost. 

In  preparation  for  data  acquisition,  the  flowmeter  was 
calibrated  using  a  graduated  container  and  stopwatch. 
Calibration  curves  were  constructed  and  the  data  fit  to  a 
cubic  equation  to  allow  computer  determination  of  the  exit 
centerline  velocity.  Analytic  expressions  for  the  density 
and  viscosity  of  water  were  also  utilized  to  provide 
Reynolds  number  and  Froude  number  data  during  computer 
calculations . 

The  alignment  of  this  Laser  Doppler  Velocimeter  proved 
to  be  extremely  sensitive  to  changes  in  location  of  the  base 
and  to  changes  in  the  angles  of  the  Bragg  Cells.  While  it 
was  not  necessary  to  repeat  the  entire  alignment  process  for 
each  experiment,  it  was  essential  that  beam  concentricity, 
beam  crossing  and  power  maximization  be  checked  if  signal 
quality  deteriorated.  Once  beam  crossing  had  been  achieved 
at  one  LDV  orientation,  the  device  was  rotated  90  degrees  to 
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assure  that  the  beams  still  crossed  and  that  the 


intersection  volume  had  not  changed  in  its  position  relative 
to  the  optics.  The  latter  was  easily  determined  by 
positioning  the  intersection  volume  on  a  scatter  plate, 
marking  the  spot,  and  then  rotating  the  barrel  of  the  LDV. 
If  the  laser  light  moved  off  the  spot  during  the  rotation, 
realignment  was  necessary  because  of  lack  of  concentricity. 
Loss  of  beam  crossing  was  apparent  by  a  total  loss  of 
Doppler  signal. 

Once  calibration  and  alignment  checks  were  complete,  one 
step  remained  in  preparation  for  data  acquisition.  This  was 


the  establishment  of 

the  spatial 

reference 

for 

the 

intersection  volume. 

The 

center  of 

the 

nozzle 

exit 

was 

chosen  as  the  origin 

of 

the  coordinate 

system. 

It 

was 

located  by  finding  the  coordinates  of  the  extrema  of  the 
nozzle  structure  in  all  three  directions  and  then  applying 
knowledge  of  the  nozzle  geometry  to  find  the  center.  The 
edges  of  the  nozzle  were  found  by  observing  the  filtered 
output  of  the  photomultiplier  signal  on  the  oscilloscope. 
When  the  intersection  volume  was  precisely  centered  on  the 
knife  edge  of  the  nozzle,  the  signal  observed  on  the 
oscilloscope  was  a  sinusoid  of  frequency  equal  to  the 
applied  frequency  shift.  A  maximum  amplitude  indicated  the 
center  of  the  intersection  volume  was  at  the  reflective 
surface. 
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After  initialization  of  the  data  acquisition  program  on 
the  HP9826,  the  LDV  was  positioned  in  the  test  section.  The 
frequency  shift  was  selected  so  that  the  whole  range  of  flow 
reversals  anticipated  could  be  measured.  For  the  present 
work  a  frequency  shift  of  500  kHz  provided  excellent 
results.  Next,  the  low  filter  was  adjusted  to  remove  the 
pedestal  while  observing  the  bursts  on  the  oscilloscope. 
Occasionally,  the  pedestal  was  not  visible  through  noise  in 
the  signal  initially.  In  that  case  the  high  filter  was 
adjusted  first  to  eliminate  enough  to  see  the  Doppler 
signal.  Final  filter  adjustment  was  largely  a  matter  of 
judgement  with  respect  to  how  wide  the  bandpass 
characteristics  should  be  set  to  pass  the  range  of  signals 
characteristic  of  a  turbulent  flow. 

After  the  signal  quality  had  been  optimized  by  the 
filter  adjustment,  the  data  acquisition  proceeded. 
Provision  has  been  made  in  the  data  acquisition  program  on 
the  HP9826  to  display  the  "spectrum"  of  the  velocity  data 
being  accumulated.  The  display  also  can  be  be  zeroed  at  any 
time  to  erase  accumulated  data.  The  display  was  monitored 
to  observe  the  distribution  of  the  velocities  to  monitor 
data  quality.  If  several  peaks  appeared,  noise  was  present 
(the  distribution  should  be  nearly  Gaussian  in  shape),  and 
the  filters  were  reset  to  narrow  the  passband.  In  some 
cases,  the  setting  intervals  on  the  filters  were  not 
sufficiently  fine  to  narrow  the  band  in  the  frequency  range 
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currently  in  use.  In  that  case,  a  change  in  the  value  of 
the  frequency  shift  was  appropriate,  and  the  frequency-  shift 
and  filters  acted  together  to  provide  appropriate  range. 

When  signal  quality  and  the  distribution  of  the 
accumulated  data  indicated  that  good  quality  information  was 
being  received,  the  operator  stored  the  information  at  that 
point  and  proceeded  to  the  next.  When  to  store  data  was 
determined  by  the  number  of  samples  desired  for  the 
statistical  base  and  the  quality  of  the  signals  being 
processed.  In  the  present  work,  the  statistical  base 
desired  was  a  minimum  of  100  data  points  at  each  location. 
Early  experiments  provided  good  experience  in  filter  setting 
and  selection  of  the  frequency  shift,  and  little  trouble  was 
experienced  in  later  experiments  in  obtaining  both  good 
signal  quality  and  a  Gaussian  distribution  of  the  velocity 
spectrum  during  acquisition. 

Several  experiments  could  be  conducted  without  the 
necessity  of  completely  reducing  the  data  for  any  previous 
ones.  The  time  required  to  gather  data  on  one  velocity 
component  in  a  jet  at  five  levels  and  seven  to  nine 
positions  on  each  level  was  about  one  day.  The  number  of 
points  which  were  measured  for  each  experiment  was  limited 
simply  by  the  time  over  which  conditions  remain  constant  in 
the  apparatus. 
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H.  DATA  REDUCTION 


Computer  programs  were  developed  for  data  transfer  and 
reduction  as  well  as  the  acquisition  control  process.  These 
programs  were  modularized  to  the  extent  practical  in  order 
to  provide  maximum  flexibility  in  application  to  later 
research.  Functionally,  the  programs  can  be  classified  as 
acquisition,  transfer,  intermediate  reduction,  and  final 
reduction. 

The  acquisition  program  resided  exclusively  in  the 
HP9826.  Its  function  was  to  automate  the  gathering  of  the 
experimental  data  and  to  convert  raw  time  data  into 
velocity.  Additionally,  the  acquisition  program  recorded 
the  relative  coordinates  at  which  each  sample  was  taken,  the 
orientation  of  the  LDV  for  the  data  set,  the  flow  meter 
reading,  the  temperatures  of  effluent  and  bulk  fluids,  the 
date  of  the  experiment,  and  identifying  information  for  the 
data  set.  Its  final  output  was  a  set  of  files  on  flexible 
disk  numbered  by  experiment  number  and  data  set  within  the 
experiment. 

Data  transfer  programs  were  necessary  on  both  the  HP9826 
and  the  VAX  computer.  They  acted  in  conjunction, 
communicating  over  telephone  lines  to  transfer  the  files 
from  floppy  disk  to  internal  storage  in  the  mainframe.  The 
process  was  automatic  and  progressed  without  operator 
intervention  once  program  execution  began  as  long  as  all 
files  to  be  transferred  resided  on  the  same  flexible  disk. 


If  a  particularly  long  experiment  was  spread  over  more  than 
one  disk,  the  operator  reinitialized  the  process  in  order  to 
change  to  the  new  one. 

The  purpose  of  the  intermediate  reduction  program  on  the 
VAX  was  to  consolidate  the  information  contained  in  the 
separate  data  files  associated  with  an  experiment  into  one 
file.  In  the  process,  the  velocity  data  was  processed 
statistically,  applying  Chauvenet's  criterion  to  reject  data 
values  which  deviated  excessively  from  the  Normal  Error 
distribution.  This  program  preserved  and  recorded  all 
identifying  data  such  as  position  and  temperature,  but 
presented  the  output  in  the  form  of  a  single  file  which 
reported  the  velocity  data  in  terms  of  mean  velocity  and 
turbulence  intensity  for  the  components  measured  at  each 
position.  Because  this  program  expected  to  find  files  for 
streamwise  and  radial  components  in  its  input  data  set, 
dummy  files  were  included  for  the  second  component  if 
studies  were  conducted  in  the  streamwise  direction,  only. 
This  was  most  easily  done  in  the  data  acquisition  phase  of 
the  process. 

Final  data  reduction  was  accomplished  using  a  program 
which  was  applicable  to  the  study  of  axisymmetric  jets  in  a 
quiescent  ambient  only.  It  arrived  at  values  of  B,  the  jet 
halfwidth,  and  U^,  the  jet  centerline  velocity,  by  fitting 
the  data  from  each  trace  directly  to  an  equation  for  a 
Gaussian  profile.  The  values  of  these  quantities  for  each 
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trace  were 


then  fitted  to  power  law  curves  to  arrive  at 
expressions  for  each  as  a  function  of  the  streamwise 
coordinate. 

The  assumption  of  the  Gaussian  profile  enabled  easy 
integration  to  obtain  a  simple  expression  for  the  volumetric 
flow  rate  in  the  jet,  Q,  at  each  level  as  a  function  of 
centerline  velocity  and  jet  half  width.  As  each  of  these 
was  now  a  known  function  of  the  streamwise  coordinate, 
substitution  and  differentiation  with  respect  to  S  yielded 
directly  an  expression  for  the  entrainment  coefficient.  The 
coefficients  obtained  from  the  power  fits  were  used  to 
compute  local  values  of  a  and  to  arrive  at  a  mean  value  and 
standard  deviation  for  use  in  comparisons. 

The  final  reduction  program  also  computed  the  values  of 
the  exit  Froude  Number  and  Reynolds  Number  using  properties 
calculated  from  accepted  analytic  expressions.  All  results 
were  tabulated  in  program  output  including  coefficients  of 
determination  for  curve  fits  and  values  of  standard 
deviation  for  quantities  for  which  means  were  computed. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Ten  experiments  were  conducted  during  the  course  of  the 
present  work.  Additionally,  many  sample  traces  were  made 
and  trial  measurements  were  taken  for  the  purpose  of 
selecting  the  optimum  seeding  particle,  optimizing  filter 
settings,  and  testing  of  the  alignment  of  the  LDV.  Of  the 
ten  experiments,  the  first  six  served  to  develop  the 
technique  to  be  used  and  to  debug  both  the  apparatus  and  the 
data  acquisition  computer  program. 

The  final  four  experiments  were  conducted  using  the 
procedures  refined  in  the  previous  six  and  constitute  the 
trial  set  for  comparison  with  previous  work  and  validation 
of  the  technique.  Items  of  specific  interest  in  these 
experiments  are  the  following: 

1.  The  verification  of  the  assumed  Gaussian  velocity 
profile  over  the  radial  cross  section  of  the  jets. 

2.  The  decay  of  the  jet  centerline  velocity  with 
streamwise  coordinate  within  the  Zone  of  Established 
Flow. 

3.  The  growth  of  the  jet  through  entrainment  of  the 
surrounding  fluid  as  demonstrated  by  the  variation  of 
jet  half  width,  B. 


4.  The  mean  entrainment  coefficient  for  each  jet,  as 
computed  from  the  local  values  at  each  trace .  level 


measured. 

Each  of  these  will  be  addressed  separately  in  the 
following  sections,  and  numeric  results  are  tabulated  in 
Appendix  B. 

A.  VERIFICATION  OF  THE  GAUSSIAN  PROFILE 

Velocity  traces  were  obtained  by  measuring  five  to  ten 
mean  velocities  across  the  diameter  of  the  jet  at  various 
downstream  positions.  Each  trace's  data  was  fit  to  a 
Gaussian  curve  during  the  first  step  in  the  final  data 
reduction  program.  Figure  9  is  a  plot  of  one  trace's  least 
squares  fit  curve  and  the  data  from  which  it  was 
constructed.  This  trace  is  typical  of  the  data  gathered, 
and  illustrates  both  the  goodness  of  the  curve  fit  and  the 
relative  uncertainty  as  to  the  actual  position  in  the  jet  at 
which  the  measurements  were  made. 

While  the  centerline  of  the  jet  nozzle  is  considered  to 
have  been  located  quite  accurately,  examination  of  the  data 
for  the  Gaussian  fits  and  the  curves  themselves  provides 
indication  that  the  jet  may  not  have  been  discharged  exactly 
vertically.  This  could  lead  to  some  error  in  the  maximum 
value  of  the  centerline  velocity  obtained  in  lower  traces, 
but  the  error  will  decrease  in  higher  levels  of  the  jet  due 
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to  the  spreading  of  the  jet  and  the  attendant  decrease  in 
the  rate  of  change  of  the  jet  velocity  radially. 

As  can  be  seen  from  the  sample  trace,  the  assumption  of 
a  Gaussian  profile  is  valid  to  within  the  precision  of 
measurement  possible  in  a  turbulent  flow. 

B.  THE  GROWTH  OF  JET  HALFWIDTH 

Figure  10  presents  the  behavior  of  the  jet  width  as  a 
dimensionless  quantity  plotted  versus  a  non-dimensional 
distance  down  the  jet  axis.  Experiments  8  through  10  are 
remarkably  consistent,  particularly  at  higher  values  of  the 
streamwise  coordinate.  The  values  of  the  data  for 
experiment  7  deviate  in  both  slope  and  intercept  for  reasons 
which  will  be  discussed  later  in  the  chapter. 

For  comparison,  the  predicted  value  of  the  quantity  B/D 
using  the  relation  proposed  by  Albertson  [8]  is  plotted. 
The  results  of  experiments  eight  through  ten  indicate  a 
uniformly  greater  width  than  the  model  for  jets  of  infinite 
Froude  number.  This  is  to  be  expected  in  view  of  the 
increased  rate  of  entrainment  which  is  predicted  with 
decreasing  Froude  number. 

C.  THE  DECAY  OF  CENTERLINE  VELOCITY 

The  decay  of  the  centerline  velocity  found  in  the  four 
experiments  is  plotted  in  Figure  11,  along  with  a  curve 
formulated  from  Albertson's  model  for  infinite  Froude 
number.  Once  again,  the  data  for  the  final  three 


experiments  shows  remarkable  consistency,  and  that  for 
experiment  seven  deviates  substantially. 

The  data  for  experiments  eight  through  ten  demonstrate  a 
decreasing  rate  of  decay  of  the  centerline  velocity  with 
decreasing  Froude  number  as  evidenced  by  their  slopes  on  the 
plot.  The  slopes  fall  between  the  value  predicted 
(variation  as  1/S)  by  Albertson  for  F  =  <ao  and  by  the 
variation  as  S  '  predicted  by  Rouse  [14]  for  F  =  0.  The 
decrease  in  decay  rate  can  be  explained  by  the  effect  of 
increasing  buoyancy  forces  with  decreasing  Froude  number. 
The  experimental  data  presented  here  is  considered  to  agree 
very  well  with  theory. 

Extrapolation  of  the  plots  of  data  in  Figure  11  reveals 
one  more  item  of  interest.  The  definition  of  the  Zone  of 
Flow  Establishment  fixes  the  demarcation  between  it  and  the 
Zone  of  Established  Flow  as  that  point  at  which  turbulent 
mixing  reaches  the  centerline  of  the  jet.  At  this  point, 
the  centerline  velocity  will  start  to  decay.  Applying  this 
definition,  the  Zone  of  Established  Flow  can  be  seen  to 


start  well  before  the  6.2  diameters  predicted  by  Albertson. 
The  change  is  attributable  to  the  increase  in  the 
entrainment  rate  for  a  buoyant  jet  over  the  pure  momentum 
jet  studied  by  Albertson. 


D.  THE  ENTRAINMENT  COEFFICIENT 


Figure  12  plots  the  mean  entrainment  coefficients 
determined  in  experiments  seven  through  ten  with  a  curve 
computed  from  Davis's  relation  [13].  The  values  computed  for 
experiments  eight  through  ten  continue  to  conform  to 
predicted  values  while  experiment  seven  deviates. 

Examination  of  the  values  of  entrainment  coefficient 
computed  locally  in  the  jets  shows  a  decrease  with 
streamwise  coordinate  in  all  cases  except  for  experiment 
eight,  where  it  remained  relatively  constant.  For  an 
initially  buoyant  jet,  the  decrease  in  entrainment 
coefficient  with  distance  is  expected  as  the  local  buoyancy 
difference  decreases.  This  is  reflected  in  an  increase  in 
the  local  Froude  number  and  a  decrease  in  a  which  is 
predicted  by  all  entrainment  models.  Experiment  seven 
exhibits  this  behavior  even  while  not  conforming  to  the 
value  predicted  for  the  mean  value  of  the  entrainment 
coefficient . 

A  comparison  of  the  obtained  local  entrainment 
coefficients  with  theory  was  not  possible  because  of  a  lack 
of  local  temperature  data.  However,  the  maximum  entrainment 
coefficient  predicted  by  the  model  of  Hirst  [12]  for  the 
four  experiments  is  0.061,  a  value  applicable  only  at  the 
nozzle  exit  for  experiment  ten.  At  the  earliest  point 


measured  in  experiment  ten,  a  distance  of  eight  diameters 
downstream,  the  measured  coefficient  was  0.076.  Since  the 


entrainment  coefficient  can  not  increase  with  distance,  as 
buoyancy  forces  are  decreasing  in  that  direction,  data  from 
these  experiments  can  be  seen  to  be  at  variance  with  the 
model  of  Hirst. 

E.  THE  JET  OVERALL  VELOCITY  DISTRIBUTION 

Figure  13  is  the  distribution  of  velocity  throughout  the 
jet  measured  in  experiment  eight.  The  plot  was  constructed 
using  the  coefficients  of  the  curve  fits  for  B  and  Um  and 
the  Gaussian  velocity  profile.  It  is  representative  of  the 
data  obtained  in  all  four  experiments.  On  this  plot  the 
variation  of  halfwidth  and  centerline  velocity  and  the 
Gaussian  distribution  of  velocity  with  radial  distance  can 
all  be  seen. 

F.  THE  VARIANCE  OF  EXPERIMENT  SEVEN 

Experiment  seven  was  conducted  with  the  highest  exit 
velocity  of  the  four  final  trials.  The  flowmeter  setting 
for  this  set  of  data  was  forty  percent  of  full  scale,  a 
value  which  was  estimated  from  earlier  trials  to  be  low 
enough  to  eliminate  the  effects  of  recirculation  in  the 
tank.  The  recirculation  effects  at  higher  velocity  were 
observable  due  to  the  presence  of  the  alumina  seeding 
particles  in  the  tank.  It  resulted  in  an  ambient  which, 
instead  of  being  quiescent,  was  coflowing  with  the  jet 
discharge.  Entrainment  coefficients  in  early  trials  were 
significantly  lower  than  expected,  and  the  rate  of  jet 
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growth  and  decay  of  centerline  velocity  decreased  as  a 
result. 

The  phenomenon  can  be  explained  by  considering  the  cause 
of  entrainment  of  fluid  into  the  jet:  viscous  shear  at  the 
jet  boundary.  The  viscous  shear  is  a  function  of  the  rate 
of  change  of  fluid  velocity  at  the  point  in  question,  and 
that  rate  of  change  is  lower  if  the  ambient  has  some 
component  of  velocity  in  the  same  direction  as  the  jet. 
Reduced  viscous  shear  results  in  reduced  entrainment,  which 
in  turn  results  in  a  decrease  in  the  rate  of  decay  of 
centerline  velocity  as  well  as  jet  growth. 

While  no  recirculation  was  apparent  in  the  tank  during 
data  acquisition  for  experiment  seven,  the  behavior  of  the 
jet  indicates  that  some  amount  was  present.  Specifically, 
the  lower  rate  of  decay  of  centerline  velocity  and  the  lower 
growth  rate  of  the  jet  dimensions  indicates  the  lower  rate 
of  entrainment  into  it.  This  is  directly  reflected  in  the 
much  lower  values  of  a  found. 

G.  RECOMMENDATIONS  FOR  FURTHER  STUDIES 

The  length  of  the  Zone  of  Flow  Establishment  is  a  topic 
for  which  no  predictive  relations  have  yet  been  developed. 
The  Laser  Doppler  Velocimeter  is  a  tool  which  can  be 
effectively  used  to  study  flow  field  changes  during  the 
transition  from  establishing  to  established  flow.  It 
possesses  the  ability  to  be  positioned  precisely  and  will 
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not  affect  the  flow  characteristics  because  of  its 
non-intrusive  nature.  Further,  the  intersection  volume  can 
be  made  small  enough  through  the  selection  of  the  proper 
combination  of  beam  spacing  and  focal  length  to  provide 
extremely  fine  spatial  resolution  of  measurements. 

The  body  of  literature  on  turbulent  buoyant  plumes  in 
cross  flow  environments  currently  in  existence  assumes  axial 
symmetry  of  the  jet.  In  fact  this  is  not  the  case.  The  LDV 
can  be  used  to  fully  characterize  the  flow  field  in  this 
environment.  For  this  use,  however,  a  very  fine  spatial 
grid  will  have  to  be  studied  in  order  to  yield  the  third 
component  of  velocity  through  the  use  of  the  continuity 
equation.  This  requirement  renders  an  exact  knowledge  of 
the  trajectory  of  the  jet  mandatory  in  order  to  optimize  the 
measurement  procedure.  The  trajectory  of  the  jet  could  be 
determined  by  dye  injection  or  by  investigation  of  the 
temperature  field  to  determine  in  advance  the  grid  of 
positions  at  which  the  velocity  should  be  measured. 
Temperature  data  throughout  the  jet  would  also  be  useful  in 
developing  exact  correlations  involving  the  local  Froude 
number,  which  would  be  a  more  appropriate  quantity  than  the 
exit  Froude  number  in  characterization  of  local  behavior 
within  the  jet. 

In  consideration  of  the  apparent  differences  reported  in 
this  report  with  the  model  of  Hirst,  further  studies  on 
buoyant  jets  discharged  into  a  quiescent  ambient  should  be 
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conducted  to  confirm  the  difference  and  to  propose  an 
improved  value  for  a2*  Hirst's  empirical  constant  in  the 
relation  for  entrainment  coefficient  as  a  function  of  Froude 
Number.  To  accomplish  this  with  accuracy,  a  larger  tank  and 
improved  system  for  flow  control  should  be  installed  to 
improve  reproduceability  of  the  settings  and  provide  for  a 
more  nearly  infinite  ambient  into  which  the  jet  may  be 
discharged.  A  larger  tank  would  also  obviate  the  need  for  a 
system  for  the  removal  of  heat  from  the  tank,  thus  removing 
the  artificial  temperature  gradients  introduced  to  the 
system  by  the  coil.  A  system  for  measurement  of  the 
temperature  field  should  also  be  installed  to  provide  for 
determination  of  the  local  Froude  number. 

H .  SUMMARY 

The  Laser  Doppler  Velocimeter  is  a  viable  and  accurate 
tool  for  measurement  of  the  velocity  distribution  in 
turbulent  buoyant  jets.  The  studies  conducted  during  the 
present  work  showed  good  agreement  with  the  results  obtained 
by  other  methods  in  separate  studies .  Inaccuracies  in  the 
results  can  be  attributed  to  non-availability  of  specific 
data  (e.g.  local  temperature)  or  to  imprecision  of  the 
positioning  mechanism  or  metering  devices  used,  and  not  to 
the  LDV  itself. 
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APPENDIX  A 
ANALYSIS  OF  ERRORS 

Error  analysis  for  the  present  work  is  accomplished  in 
three  levels.  The  first  of  these  is  the  determination  of 
the  precision  to  which  basic  quantities  such  as  spatial 
location  and  time  can  specified.  Where  available,  the 
tolerances  reported  in  the  manufacturer's  specifications  for 
the  device  have  been  used  as  the  source.  If  the  quantity 
can  not  be  determined  from  the  technical  documentation,  the 
precision  has  been  estimated  from  the  fineness  of  scale  or 
reference  used  in  the  operation  of  the  equipment.  The  basic 
quantities  and  their  estimated  error  are  as  follows: 

1.  Time  (Counter-Tracker),  1  nsec. 

2.  Time  (flowmeter  calibration),  0.1  sec. 

3.  Volume  (flowmeter  calibration),  5%  of  the  quantity 
measured  (1000  ml) 

—5 

4.  Spatial  location,  0.001  inches  (2.54  X  10  meters) 

5.  Beam  separation  at  the  transmitting  lens,  0.1  mm. 

6.  Focal  length  of  the  transmitting  lens,  0.05  mm. 

-12 

7.  Laser  light  wavelength,  5.0  X  10  m. 

For  the  second  level  of  analysis,  the  uncertainties  of 
the  variables  which  were  computed  from  formulas  involving 
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the  basic  quantities  were  estimated  using  the  method  of 
Kline  and  McClintock  [16].  The  uncertainties  are: 

1.  Radial  position  in  the  jet  (left  to  right),  0.110  mm. 

2.  Streamwise  position  in  the  jet,  0.110  mm. 

3.  Radial  position  in  the  jet  (in  and  out),  0.205  mm. 

4.  Mean  flow  velocity,  (flow  meter),  5%  of  the  scale 
reading. 

5.  Beam  crossing  half  angle,  0.000819  radians. 

6.  Frequency  of  the  Laser  Doppler  signal,  62.5  Hz  at  a 
signal  frequency  of  1  MHz. 

7.  Particle  velocity,  0.032  m/sec  at  a  signal  frequency 
of  1  MHz,  corresponding  to  a  velocity  of  3.50  m/sec. 
or  approximately  1%  of  the  velocity. 

For  the  third  level  of  analysis,  the  root  mean  square 

error  between  the  curve  fit  computed  values  and  the  input 

data  values  for  U  and  B,  both  at  the  trace  to  trace  level 

m 

and  throughout  the  jet  was  computed .  The  standard  deviation 
was  also  computed  for  the  value  of  the  mean  entrainment 
coefficient.  The  error  of  the  predictive  equations  in 
reproducing  the  data  was  on  the  order  of  ten  percent  of  the 
value  of  the  data  point. 

The  percentage  error  expected  in  the  values  of  the 
measured  velocity  using  the  LDV  is  seen  to  be  less  than  one 
percent.  The  quantity  is  very  sensitive  to  the  precision  of 
the  alignment  of  the  optics.  The  dominant  factor  in  the 
computation  is  the  error  in  beam  separation  distance. 


Accordingly,  it  should  be  noted  that  extreme  care  must  be 
taken  while  crossing  the  beams  to  move  only  in  a  direction 
perpendicular  to  the  plane  in  which  the  beams  lie.  Movement 
within  the  plane  of  the  beams  will  alter  the  crossing  angle 
and  subject  the  measurement  process  to  substantial  systemic 
errors . 

The  reproduceability  of  the  positioning  data  is 
considered  to  be  good,  but  the  uncertainty  of  the  actual 
position  of  the  intersection  volume  in  the  jet  itself  is 
considered  to  exceed  the  values  quoted  for  spatial 
accuracy.  The  reason  for  this  is  the  indication  of  a 
deviation  from  vertical  of  the  jet  trajectory  by  as  much  as 
one  degree.  Accordingly,  the  radial  location  reported  for 
the  jet  data  points  is  considered  to  be  truly  accurate  to 
within  no  better  than  one  millimeter.  This  is  considered  to 
be  the  principal  cause  of  error  in  the  fitting  of  the  data 
to  Gaussian  curves. 

The  probable  presence  of  circulatory  effects  in  the  tank 
can  be  considered  to  be  the  primary  cause  of  the  deviation 
of  the  data  for  experiment  7  from  predicted  behavior.  It 
should  be  noted  that  an  "inversion"  of  the  behavior  of  both 
centerline  velocity  and  jet  growth  occurred  in  the  body  of 
the  jet.  Specifically,  the  centerline  velocity  increased 
from  one  level  to  the  next,  and  the  half  width  decreased. 
Given  the  accuracy  of  the  LDV  in  determining  velocity,  it  is 
considered  that  the  trend  accurately  reflects  the  actual 
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behavior  of  the  flow  field  and  that  the  change  was  the 
result  of  forces  external  to  the  jet  itself.  Similar-  local 
inversions  of  the  expected  trends  were  observed  in  jets  of 
earlier  experiments  in  which  recirculation  within  the  tank 
was  clearly  visible. 


APPENDIX  B 
TABULATED  RESULTS 

Tabulated  below  are  values  of  interest  which  were 
obtained  from  the  curve  fitting  routines.  The  local  values 
of  the  centerline  velocity  and  jet  half-widths  are  computed 
using  the  Gaussian  profile  assumptions.  They  were  in  turn 
fitted  to  the  power  curves  yielding  the  equations  given. 
The  coefficients  were  then  used  to  compute  the  value  of  the 
entrainment  coefficient. 

Experiment  7 

Exit  Froude  Number  456  Exit  Reynolds  Number  7836 

Nozzle  Exit  Centerline  Velocity  2.54  m/sec 


S  (mm) 

U_  (m/sec) 

B(mm) 

111 

25.4 

.806 

3.97 

50.8 

.725 

5.79 

76.2 

.415 

9.55 

101.6 

.457 

8.  74 

127.0 

.333 

12.88 

Power  curve  fit 

equation 

Coefficient  of  Determination 

U  -  .116S"*549 
m 

.852 

B  *  . 0505S  *  699 

.934 

Mean  entrainment 

.  coefficient 

.049 

Standard  Deviation 

.012 

Experiment  8 


Exit  Froude  Number  365  Exit  Reynolds  Number  7004 

Nozzle  Exit  Centerline  Velocity  2.22  m/sec 


S  (mm) 

U  (m/sec) 
m 

B(mm) 

25.4 

.983 

3.05 

50.8 

.543 

6.59 

76.2 

.376 

9.45 

101.6 

.457 

12.02 

127.0 

.236 

16.70 

Power 

curve  fit  equation 

Coefficient  of 

U  -  . 
m 

0408S-*867 

.996 

B  *  . 1327S1 *  023 

.996 

Mean  entrainment  coefficient 

.073 

Standard  Deviation 

.007 

Experiment  9 

Exit  Froude  Number  317 

Exit 

Reynolds 

Nozzle 

Exit  Centerline  Velocity 

1.91 

m/  sec 

S  (mm) 

U  (m/sec) 
m 

B(mm) 

25.4 

.680 

3.92 

50.8 

.534 

5.80 

76.2 

.327 

9.24 

101.6 

.260 

12.94 

127.0 

.196 

17.82 

Determination 


Number  5954 


Power  curve  fit  equation 


Um  =  .0429S"*783 
B  =  .ins*940 


Coefficient  of  Determination 
.944 
.965 


Mean  entrainment  coefficient  .072 
Standard  Deviation  .008 
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Experiment  10 


Exit  Froude  Number  239  Exit  Reynolds  Number  5123 

Nozzle  Exit  Centerline  Velocity  1.59  m/sec 

S  (mm)  U  (m/sec)  B(mm) 

m 


25.4 

.527 

50.8 

.427 

76.2 

.275 

101.6 

.203 

127.0 

.173 

Power  curve  fit  equation 
U  =  .0409S-'726 
B  «  .0974S'897 


3.78 
6.17 
9.24 
13.90 
14.93 

Coefficient  of  Determination 
.940 
.983 


Mean  entrainment  coefficient 
Standard  Deviation 


.069 

.009 


APPENDIX  C 


SOURCE 


FIGURES 


PARTICLE 


Figure  1.  Geometry  of  Doppler  Shift  Model 


LASER 


APERTURE  AND  REAR  ROTATING  MOUNT 
BEAMSPLITTER  MODULE 

FREQUENCY  SHIFTING  MODULE  WITH  BRAGG  CELLS 
BEAM  STEERING  MODULE 
RECEIVING  OPTICS  ASSEMBLY 

6.  BEAM  EXPANDER 

7.  BEAM  STOP  ASSEMBLY 

8.  TRANSMIT/RECEIVE  LENS 

9.  PHOTOMULTIPLIER  TUBE 


Figure  5.  LDV  Optics  Arrangement 


Fig  6.  Doppler  Bursts 


Figure  7.  Signal  Flow  Path 


TEST 

TANK 


s*s  ELECTRIC  HEATER 


VT7  FLOWMETER 


Fig  9.  Jet  Velocity  Profile 


Fig  10.  Growth  of  Jet  Width 


Dimensionless  Axial  Distance,  S/D 
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Fig  11.  Velocity  Decay  of  Jet 


Entrainment  Coefficient,  a 

0.03  0.04  0.05  0.06  0.07  0.06  0.00  0.10 


Fig.  12  Entrainment  Coefficient 


LEGEND 

□  Experimentally  Determined 
—  Davis  et  al. 


300  350  400 

Froude  Number 


< 


? 


A 

4/ 

tT 

r  *'4 

0£ 

X 

1 

'u' 

UJ 

o 

o 

o 

x 

CO 

u_  >c 

o 

i— 

ir 

x 

c 

x 

• 

UJ  ZZ  UJ 

* 

EZ 

HH 

t— 

_J  X 

■ 

3 

at 

cc 

O  UJ 

©  « h- 

04 

DC 

S\ 

a 

a.  x 

'X  X 

\^‘ 

f — 

© 

o 

Z2  -  X  X 

N 

o 

© 

x 

u _ 1 

Oh  O 

a 

UJ 

CVI 

o 

O  _J 

1  “7“  ^ 

4> 

X 

V 

cc 

3  3  >- 

,  crt 

to 

CS 

co 

Ql 

coz:  »  ■ 

1 

♦> 

t— 

UJ  UJ 

(L-XZ 

£ 

« 

Cf 

x 

CO  _J 

J-  X  UJ 

cs 

X 

TJ 

» — t 

x  *-< 

X  X  •  DC  UJ 

4» 

X 

1 

o 

u_ 

UhCOt 

CD 

a. 

X 

— •  X  Ul  O 

HC 

~7~ 

■V  u> 

• 

Ixi 

ujqui  co 

■ 

o  • 

Ul  uj 

X  CO 

X  X  Ul 

04 

m  In- 

:  x  uj  in 

>  X  Cl.  I 

I -  UJ 

I  CO  ’T 

©  OJ  ■ 

■  ©  <S >  00  -* 

©  ©  ■  04 


Ul  I—  _l  CO 
X3  coo 
3  l  J  X  X  >-• 
ZZ  jl-I 
hh  I—  X 
X  CO  I —  O  'I 
•-«  UJ  _l  X 

x  x  co  x  u 


3 

1 

CO  — < 

X 

04 

04  CO  — <  -*-• 

X  ©  X 

X 

X  oo 

r 

tO 

o 

J— 

II 

II  II  II  OJ  II 

--  II  II 

(=> 

TT“ 

^  d 

toco 

CO 

© 

'A 

hZZCM  C 

II  •-  c 

CQ  CQ 

c 

X 

v  Hi 

•— 1 

X 

© 

«-> 

V  L>  ti  II 

oo  ©  vt  d  o 

X 

— i  ~ i 

0C 

♦->  — « 

X 

o 

r 

•h  O  r p  a  S 

©  CO  ©  <->  4'  ■■■« 

H-4 

CO  CO 

£• 

o 

o 

o  uC 

o 

II 

•rt 

X  C  SZ  4<  1 

O  — «  w  ih  il*  »— *  l.A 

■J 

o  o 

V 

o 

JO 

QC  UJ 

'rt 

o 

■A  Of  4*  <A  «-> 

II  II  II  II  4)  U  L. 

—  U  CO 

a: 

3  CO 

CtC  J— 

• 

«-> 

a 

1  ^  IX 

X  >  X  >  1  1  4<  -* 

© 

X 

3 

UJ  X 

HH 

c 

i 

1?  !•  1  S  IT* 

i  i  i  i  d  *>  x  z>  -< 

/ 

X  CO 

H-H 

•H 

;< 

4»  Z>  O  d 

£ZCXX**SZLUiZXII 

c 

i-i  O 

xac 

o 

d 

L  b  O  (Ml 

d  d  d  OX  O  3  o 

•H 

— 

—  «U- 

ox- 

—  — 

X 

X 

X  X  X  X  X 

ZZZZRfi£ilO£L'l 

cn 

1 

©  ©  ©  ©  ©  ©  ©  ©  © 
*-•  O.J  CO  'T  ID  '  £>  r  -  ©  ON 


Ul 

• 

• 

3 

~T~ 

El 

H- 

•3 

3 

M 

/ 

as 

3 

1 - 

c 

i — 

O 

• — 1  . 

0 

0 

X 

to 

■H 

Ul 

to 

w 

♦> 

a_ 

3 

(A 

co 

i=t 

<3 

/ 

(U 

X 

O 

3 

Ul 

X 

X 

■  w 

Tt  Z 

x 

z 

S 

1  3 

i— 

Ul 

3 

X 

T3  4- 

OS 

4. 

h- 

C  «J- 

>- 

3 

C  3 

X 

4fc  *J- 

O  O 

X 

13 

Al  O 

« 

-  c  «■»  ♦» 

O  Al 

3 

M 

Cu  Al 

d3IC 

Al  G. 

• 

a. 

u_ 

U  £ 

Ul 

O  «-•■  OD'-i 

COCO 

CO 

co 

• 

O  (A 

X 

z  *>  d  0  -h  0 

Ul 

1 — 1 

Ul  X 

(A  1 

1— 

►-4-3  1  a-  a 

CQ  CD 

~7~ 

a 

X  3 

1  3 

to  A'  1  — •  JZ  1 

3  3) 

1 — 1 

1—  X 

3  O 

2 

3  3  A'  'll  13 

co  to 

H- 

I— 

rH  'll  4. 

»H 

C^ijj  ♦>  O' 

OO  ♦* 

3 

3  0 

II  i-iZO 

L3 

(VI  0  d  d  41  4. 

L3U  >ll 

O 

>r  ui 

•A  il 

Ul 

-« 0 1-  _l  UJ  G 

■A 

or 

as  x 

ul  CQ 

X 

—  — 

1 

CQ 

•• 

(=i  co 

C  C  3  3 

OS  X  X  X  X  X 

3 

3 

■— 

3  CO  CO 

UJ33333 

'll 

CO 

3 

• 

O  OOO 

• 

1—  CO  CO  CO  03  CO 

~T~ 

4. 

(L'JUIJ 

CM 

••  x  0  0  0  0  0 

£  UJ  i J  L 3UUU - 


S  O  O  O  O  Q  O  O  Q  O  ©  O  O  '3  O  <3  ©  O  *3  Q  ©  O  ©  «3  ©  ©  ©  ©  ©  ©  ©  © 

» cr*  ©  *-«  o  j  co  *rr  in  *0  r>-  ©  ©  ©  0  j  co  rr  in » o  r-  co  ©  ©  ~«  oj  co  ^  uo  ©  n.  © 

*>  ©  rr  rr  rj-  rr  -r  t  -t  -t  rr  in  10  in  in  in  in  in  in  in  in  ©  kd  *d  vd  id  *d  *0  ‘•d 'd  *d 


•  '.*-■  .%*'  *-•  *  •«.  • 


•  *_VO  •  'O  v‘<l4<»  1 


700  IF  Count  =26  THEN 
710  GOSUB  Find-Peak 

720  Count =0 

730  END  IF 

7* 4 V\  f‘ r.u^*  a.1 


■  -  *  w  r* 

* 

K 

Li 

1 


■.  -w*  ■:•  v 


‘•*vl 


m 

fe 

-.> 

a 


tsj 

m  N 
•>: 
3 


m 

!i-*j 


I 

■j 


?: 


IV 

m 


I-  ♦> 

•>  c 

2  *h 

3  O 


LU 

o 

<U 

4-  3 

C  u-  '!>  C 

il'  *1*  +*  SI  'll  'll 

il)  r—i  <*_  [T*  £  L  Li 

Ct<  liP-H  J  <*-  <*- 

O  <*.  — •  L  Oil 

V)  1  1  1  O  ‘Il 

1  'li  il'  ■!)  1  O  If 

3  ♦*  3  3  *>  O  *1>  C 

■1)  -it  O  O  '-*  L_  JZ  C' 

1^ 

x  o 

O  1— 

CD 

TT  L  7“  rti  f  1 

U  ' 

o 

3  *3  r  -  i  o 

CO  CD  CD  CO 

CO 

O  1- 

1 — 

=1  O  D  Z  m  CD  D  Cli 

i-  x 

CO 

to  1—  CO  CC*  3*  3*  00  X 

O  o  O  o  to  CO  O  00 

»—i 

a  o 

>- 

UUUUOOUO 

LU  CL. 

Ul 

13  13  (J 

to 

i<£ 

s  s  :  s  : 

xa:  - 

• 

Kill  1  A:  i  %L  : 

1- IL 

-ICO 

cc i—  i  i  i —  ouo 

a;  a:  co 

a:  i— 

LU  »-«  1  1  mQ^UJUJ 

LU  f=l  *-* 

H-.QL 

_l  DC  1  1  3UIIE 

i —  a 

«— 

•te 

o  x 

0  3^  1  GNUIl 

X  Ll. 

•If 

LU  DC 

xox 

to 

Ql  DC 

o  o 

CO  LU 

_ 1  — 1 _ 1 _ 1 _ 1 _ 1 _ 1 _ 1 

O  DC 

a 

1- 

lit  111  111  lit  111  III  ll!  Ill 

lu  a 

i.u 

IUZ 

CQi^lQlOCQlBGQlX) 

LU  CO  LU 

i 

3:  •-* 

I  'X  'X  >X  I  I  'I  'X 

X  x  o 

3 

H 

_ 1 _ 1  — 1 _ 1 _ 1 _ 1 _ 1 _ 1 

1—  X  *x 

— •  to 

CO 

X  _J 

II  X 

Q_  CC 

o— *oorrr^*MinoN 

O  ql 

to 

X 

cl  lu 

+*  CQ 

(=i 

>-  >-  > 

LU  X  LU 

•Z  X 

1—  LU 

111  111  III  III  111  111  111  111 

M  1—  CD 

•*■•  O') 

LU  CO 

o  o 

CO  X 

—  —  —  —  — 

-  0-  u 

— »»  — 

■•oooooooo 

■  ■ 

30lC 

)  II  -r* 

■  ■ 

Of 

I  irt  cf 

to 

«-> 

»  ♦»  T- 

>■ 

C 

r  £zaz 

Ol 

3 

O 

>  O  O  X 

1 

X 

o 

CL 

)  CL  1 —  DC 

a 

DC 

1 

X 

OX 

3 

X 

o 

1 — 

Ul- 

«■> 

1- 

LU 

LU 

•ll 

LU 

•1; 

CL 

CL - 

- CO 

CL - 

- N 

i  0  O  Q  Q  Q  Q  O  Q  Q  0  Q  Q  O  CD  ©  O  O  O  0  0  Q  ©  *3 

'  lO  1 D  N  >30  OS  0  ^  CM  CO  rr  If)  1  £•  r-  HO  OS  ©  CM  CO  rf  in  '  £|  r-  CO 

rv.  fs.  is-  r>-  oo  co  oo  oo  co  oo  co  co  oo  co  ctmt>  cr»  <y*  ctmt>  •t*  «x>  a- 


©  © 
©  © 
a  *.  ©  © 
JS  ^ 


©  ©  ©  © 
•>i  o’)  *r  in 
©  ©  i©  o 


92 


V,  / 


1060! 

1O70! 

1080  Move_lef t : ! 

1090  !  MOVE  THE  SPECTRUM  DISPLAY  CURSOR 

1100  !  ONE  STEP  TO  THE  LEFT. 


1120  l_choice=I_old-l 

1130  IF  I .choice >135  OR  l_choice<0  THEN  l_choice=0 

1140  I_old=I_choice 


X 

y 

3 

*A 

00 

O 

Ul 

1—  (=1 

+-> 

V./ 

3 

X 

T**  ~7~ 

o 

a: 

1— 

-«  X 

•  rt 

■3 

3 

O  X 

0 

3 

O 

3 

X 

X 

3 

•v-H 

V 

1— 

•■w 

+ 

■ 

>X  U_ 

z 

X 

x 

o 

1—  Ul 

O 

■3 

3 

•H 

X  3 

ii 

■ 

a. 

O 

n 

X 

mo 

3  • 

-C 

X  • 

"?* 

3 

uo 

o 

null 

3 

3 

0-1 

t=i  x 

1 

00-0. 

O 

ST 

X 

U 

M 

a.  x 

o 

3  - 

/-S 

E  »— » 

X  3  O' 

r 

X  3 

/-S 

3  X 

X 

Ul  IJ 

3 

3 

X 

o 

LU  Ul 

>-N 

LU  Mrt 

h-  3 

♦> 

ZIUI 

XX 

+ 

■3  X 

+  uo  c 

3  1-  X 

w 

1— 

Ml 

r 

111  1— 

3  3  «1<  -rH 

X  1 — 

ii 

♦-> 

•H 

a. 

rH  *— «  |J  0 

-•*" 

V 

3  3 

c 

* 

o 

3  >3 

o  An  a 

OImIl 

>* 

v—  *— % 

•rt 

1 

a 

1— 

1  'll  O  1 

X  Q 

>; 

3  (=1 

0 

i 

Ul 

►-<0X3 

i—  a 

3 

X 

X 

1 

3 

x  a. 

II  ^  o  O 

Ul  O' 

X 

X  X 

V 

< 

O 

1-  Ul 

01  O  IX 

3  1 —  Ul 

1 — 

•3  3 

:* 

** 

c* 

t— 

O  X  I—  3 

Ul  X  X 

O  X 

3 

CO 

UJ  3 

•»h  o  :: 

3  3  X 

1— 

V 

0  13  3 

X  O  3 

21 

—  —  — •  — 

-  V 

V 

CO 

b  x 

£MrH^ 

3  0  0 

Hrt 

•  • 

1— 

< 

3 

x  o 

o  o  oo 

X 

Ml 

c 

oo 

111  13 

- - — - - — 

-  X 

V 

M 

r“ 

o 

— -  —  — •  — 

.  W  M  W  (J 

•  • 

o 

ii 

c 

IJ 

♦> 

o 

Ml 

♦* 

C 

, — « 

3  a 

r~ 

3 

3 

o  x  z 

V* 

0 

1 

o  3  cr 

•rt 

O 

■1> 

r-l  3  c 

X 

• 

1 

T" 

X  3  U 

X 

1 

X 

X 

*■* 

— 1 

3 

3 

•Ii 

2' 

ii 

»— 

:> 

1— 

3 

1 — 

£ 

3 

o 

3 

a 

3 

0 

X 

Z2 

X 

3 

X 

O 

O  O  O  O  O  O  Q  O  O  O  O  O  O  O  O  O  O  O  0  O  O  O  O  O  O  O  O  O  O  O  O  3  O  C3 

UO  '■£>  I"--  CO  <T>  3  ^  CM  CO  rr  UO  MO  r-  CO  O'-  3  01  CO  rf  UO  <D  0-  CO  3  3  «  0-1  CO  rf  UO  M5  (V  3 

_  01  OJ  01  04  OJ  0J  OJ  01  0J  01  00  O'?  00  CO  00  00  00  00  CO  00  T  ’J-  ’J"  T  ’J"  -4-  *T  ’S’ 


■>  .  i  J 


m 


<n 

_l 

Vt- 

* 

co 

1- 

© 

H>  ^ 

* 

I- 

X  • 

Li- 

c  ** 

X 

hi  : 

X 

4>  *i 

«♦- 

»— « 

o  -w- 

LU  1- 

X 

£  •!> 

0> 

o 

a.  lli 

X  o 

n 

-rt  vi 

!«- 

Ol 

T=~ 

1—  X 

i—i 

V-  1 

i 

xx 

a 

X 

z 

»-  x 

©  > 

‘ii  >A 

a 

i 

X  = 

II  X 

r*4  •  H 

X  CS 

3 

X  vl 

1=1 

—•  G£ 

■-  © 

LUPX 

a* 

^  ♦> 

O 

a: 

l 

w  ^  ^ 

c 

-•  c 

_l  LU 

1  l- 

©  ©  C  >  d 

i 

i  — 

_i  z 

<«—  2 

C  o 

x  x  c  n  >i‘ 

X 

A  O 

x  a: 

i 

x  cu  a> 

3-  Z>  ZlL  £3 

-  Q_ 

x 

XI- 

oc  a  —• 

i;  ii  n  ii  n 

II 

D  — 

LU 

e  x 

CL  CL  ©  O  O  CL 

.  ^ 

1 : 

I-  © 

L  >— (  LU  ©  1 —  Ll  LU  — •  CO  00  O"  m  U5 

*-•  _l 

LlI  LlJ  LlI  hi  llj  DD  LU 

w  \y 

5QC 

Q£  hi 

(ra.miiia;  ©x© 

o  o  o  o  o 

o 

u  O 

2  Ll. 

Q  CO  CO 

■4-  <♦—  'H  ^  •*-H 

■4- 

^  U.  r 

Ll. 

X  o  o 

iz  — 

—  — -  —  — 

•>  M 

LlI  lj?  i j} 

1— 1  »— 1  ►— |  ►— 1  M 

H-l  I 

I  I  c 
>  M  I 
— '  I— 

II  II  II 

/\  /"S 

r-  co  it. 


+j-  -rt  ■ 

£T  ~U  ii 
Gl  I 

•h  <X«  »1« 
JD  r— f  £ 

£  N  3 
S  N  O' 
1  O  — • 

i —  x  u_  * 
ii  ii  ii 
/\  /\  . 

©  — <  eg  i 

^ 

%_-•  w 

o  o  o 


«  CM 

»• 

11—  © 

Ol 

1  X  1 —  1 —  CL 

P"l 

>  *— «  » —  •  L U 

•H 

>  CL  O  X  © 

<*- 

Id _ 1  X  ©  X 

1 

©  lU 

0; 

X  © 

H> 

1— 

•H 

UJ  © 

i- 

CO  OL - 

>.  wmm 

i  CO  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©©©©©©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  © 

<r> _ i  ©  — •  eg  eo  -g-  in  ’-o  ix  ©  ©  ©  —<  eg  eo  in  gi  ix  co  <7*  ©  — •  eg  eo  t  in  'o  ix  co  '7'  ©  — « eg  eo 

-g-  >-•  in  in  in  to  in  uo  in  in  in  in  <o  co  <o  •©  <o  «o  co  <o  >o  *o  r-  rx  r---  r'-  rx  r--  ix  ix  rx  r-  co  co  co  co 


y-% 

•H 

73 

M 

C 

\/ 

iu. 

e 

x  - 

HH 

cf 

UJ 

•w 

4-5- 

y- 

-o 

X 

O 

mw* 

a> 

© 

J— 

<*- 

73 

2 

*» 

CD 

c 

1 

<3 

rH 

UJ 

HH 

ii 

•H> 

mw. 

•  w 

i. 

X 

z 

: 

o 

•rH 

*— • 

; 

o 

O  , 

^  fi 

^  f! 

i—  ««v 

y\ 

1 — 

a 

-« 

o  — 

r 

~7~ 

i 

✓  /-S. 

4->  ^ 

*»  + 

» — * 

X 

D  = 

riz 

n  ia 

¥> 

Cf 

V  «Vk 

tj  — 

c 

Of 

X 

21 

“  M 

1  ~ 

c 

jf* 

1 — 

N 

-4 

m  i  Ip- 

<2)  ^  -H 

:4c  ci 

X 

4-> 

•n 

«.  r 

-  i-  O 

a 

•r-4 

“3 

V  'S 

C  'll  v 

O  X 

_ 

o  ■ 

O 

<  O 

•H  rH  IT 

Q  H 

H-  = 

>* 

o 

X  — 

« u_ 

0'H  h 

OJ  m  - 

V 

J—  • 

rH 

O  II 

.2: 

XU-X 

i  '-<r  o 

— « >- 

to 

j  *-» 

■a « 

II  ^  X  II 

r-H  ijj 

d  X 

> 

•/I  ' 

o  - 

L.  •  H  f— "«  rH 

•  rH  rH 

i— i  or 

If 

3 

- 

j—  f— 

Q  Ll.  X  — < 

U_  -r* 

_J  X 

C  <*- 

)l— 

X>  J— 

tiOfl  J 

®  U_ 

UJ  X 

T> 

■H  | 

.2: 

wq. 

O  s  <h. 

r 

*7 

O  >v 

•HH 

ii  i —  • — ■  • — > 

C  1-  1 

x  oj 

J— 

VihCl.  ri 

•  QC 

►h  n  a; 

_■  1—  X 

i_3  J— 

UJ  X 

e  + 

A  1*- 

>Q-H- 

o  a.  h- 

£llZ  clwZI — 

•■•■■■•  HH 

T>  i- 

X 

X  X 

i-  H  H  lu 

CO  W  HH 

X  o 

o  ™) 

~>  X 

UJ 

O  UJ 

CO  Ltl  iJZ 

1—  X 

73  If 

^  — 

U_  X 

1 —  O  X  X  X  X  3 

1  'T* 

X 

Q  <3  O  O  O  'S  O  O  O  O  '3  O  O  O  O  O  O  O  O  G  O  O  O  O  O  O  O  O  O  ©  O  O 

tf  if)  t.o  fs.  CO  •7*>  O  «-h  OJ  CO  Tf'  |f)  CO  f'-  CO  '7*'  '-J  '  C'J  CO  liO  CD  fN-  CO  '7-1  O  ^  OJ  CO  T  UO 

00  CO  CO  CO  CO  CO  'X*  'X*  '7>  >X«  iTN  CT*  '7*»  'X*  '7«  i7'  «Ij  ®  '3  '10  ©  ©  ©  O  O  O  h  *h  »h  »h  »h 

»H  r-(  tH  iH  tH  tH  tH  rH  ■*— «  rH  tH  tH  tH  rH  C'*J  C'J  C'J  C'J  C'J  C'J  C'J  OJ  C'J  OJ  OJ  OJ  OJ  OJ  OJ  OJ 


2160  0pen_file:! 

2170  !  1.  CREATE  A  FILE  NAME. 

2180  !  2.  CREATE  THF  DATA  FILE 

2190  !  3.  OPEN  THE  FILE. 


1D-A132  31 6  MEASUREMENT  OF  VELOCITV  DISTRIBUTIONS  IN  TURBULENT  JETS 

USING  A  LASER  DOPPLER  VELOCIMETERCU)  NAVAL  POSTGRADUATE 
SCHOOL  MONTEREV  CA  M  D  WESSMAN  JUN  82 

F/G  20/4 


UNCLASSIFIED 


bS 

3 

.ft 

fUiiAf 

iUi 

LL 

Hi 

liki 

|  ft  . 

1  Hi 

Lil  •{" 

**&  it) 

i i 

L 

i  Jti 

jfc;N 

IB  1  *« 

i'Ji 

III  ll 

.11'’  u 

h 

iifl 

J 

I  ill  L 

uun 

[  H  Bn 

kPl 

.if 

I 

1  IWfflflt 

if 

jl 

(i 

H 

MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1%3-A 


lb 

i_ 

I 

X 


r 

o  • 

* 

o- 

r— 4 

1 

CO 

>■ 

_J 

i- 

M 

&  -<-* 

X 

<b  O 

V 

D-  O 

a  rH 

LU 

D'  a> 

-J 

d  :> 

ISJ 

M  • 

rsj 

Ah  2 

o 

X  o- 

X 

«->  /N 

V* 

C  CO 

«-> 

j— 

O  X 

to 

X 

£!  ST 

• 

Ml 

d  -O 

CO 

1 

LU 

•*h  :  'w-  JZ  <*- 

X  I-H 

d 

rr 

UJU  O-  <7- 

+ 

X 

•#-> 

M 

*»  1  ^  O  <-•  X 

3  Oi 

h-  > 

d 

CD 

—l 

c  <b  i-  x  lu  x  in 

o  m 

XI- 

1=1 

X 

ib  —  a>  ox  i  i 

^  L 

o  — 

<& 

Z  _l 

X 

C  N  «->  X  3  X  ij" 

I  d 

OJ  o 

r 

X  X 

j- 

•-«  N  lb  LU  CO  **  ‘b 

X  X 

-HO 

i 

Cki  »-H 

X 

v-  O  Z  X  -  i- 

'  \ 

XX 

X  m 

L3  1- 

X 

H  X  3  CO  X  OX 

—  3 

ox 

1-  ^ 

O  — ■ 

o 

|2  »  OK£  O  •> 

o  o 

I> 

CC  O 

X  X 

X  =  -*•-*  -HI 

O  -H 

a  O  *»• 

CL  H4 

1— 

X  O-  X  U  X  u  o- 

—  1 

X  X 

1  irt  ■  >!■ 

X 

a  /\  a  M  ffl  TZr- 

<b  x 

XX 

s  JiiCQ  £ 

LU  LU  • 

—> 

:  X:  (=1  -OJ 

z=-  -> 

1—  J— 

cS  O  «  d 

X  X  CO 

o-  H  o-  CD  =  X 

i  -h 

o»z 

1 —  1 —  X 

IL. 

XvUWI— Ov 

.c©  o 

cox 

*>>-<  Hi 

o 

o 

X  X  V  5- 

m  .  o 

1-0 

C  X  CO 

fl-ijw 

GO  X  i-h  3-  J—  1 — 

U3  ■— 1 

21  M 

lb  1  d  1 — 

ZD  X  J— 

X  . 

Z  X  1=  X  J—  -h  U_ 

X  03  lb 

X  X 

ei)Z 

— ■  CO 

o 

3  J—  X  X  i-h  O  >-H 

I-h  3- 

CO  X 

•w  ->  ib 

J—  M  X 

>.  w 

X  LU  X  X  O  O  X  5- 

^  X  * 

X 

3  l~  ^ 

LU  03  3 

d  i— 

Z  X  D-  O  _l  03  d 

*»  >b  o 

XI— 

4i  (iCC-H 

C OXO 

pH  M 

1—  X  -XX  — i 

•H  • 

X  X 

tt.  SOU. 

a  co 

xh.x>x>>  a 

O  Cl/3 

X 

X  OCQ0 

—  —  «  —  — 

-•  in  o 

x  o  x  x  z>  o  <n 

O  d  03 

XX 

LU  O 

Q_ 

~ZZ  H  O  K  2!  *H 

-h  a:  -h 

o 

II  LU  X 

•• 

3 

— 1  x  X  -  J—  03X3 

•b  n  i 

O  *-H 

W-CQI— O 

in 

1  OJ 

X  _J  Z  X  03  X  3  1  Z>  Z>  II 

XI—  • 

lb  3  X  -h 

c 

L  A 

U1NZI-U1IQ  L- 

II  1  o 

LU  H  o 

Z  CO  X  CO 

o 

d  > 

XIMOXXLDX  d 

'll  d  L 

r\j  03  x 

dOd  03 

■H 

Qj  A 

X  O  X  O  O  I-H  X  lb 

0*  *»  QJ 

1  O  X 

ZUUI 

♦> 

— 1  X 

LUZXZXXX  — 

C-h  N 

Hfl.N 

«n 

O  S 

:::::::  (J 

d  tb  1 

(Li 

X  «  >-h  — 

3 

X  J— 

1 —  1 —  1 —  1 —  1 —  1 —  i —  CQ 

of 

3  3 

33333333 

X 

1 

CO  X 

XXXXXXXCO 

X 

X 

♦> 

O  X 

xxxxxxxo 

X 

3 

|J  IH 

mhhhhhhU 

3 

H 

l- 

J— 

UJ 

•H 

X 

X 

Li. 

* 

X  — 

9  O  ©  OOOOQOiSSIO&IS)Q<)>OQQS)OCIOCIQI3|5QOOG)OSI09S 
-*  oj  co  -r  in  'o  r'-  03  cr>  ©  -h  oj  co  t  in  'b  'o  f--  co  'T.  ©  -h  oj  co  -i-  m  i-id  r-  co  >.t.  ©  -h  oj  co  rr  in 

M  04  04  O  J  OJ  OJ  OJ  O  J  OJ  CO  CO  CO  CO  03  CO  03  OD  03  CO  T  T  T  rf  T  -T  rt  rr  T  rr  UO  UO  UO  IT>  UO  UO 

vl  OJ  OJ  OJ  OJ  OJ  OJ  04  OJ  04  OJ  OJ  OJ  OJ  OJ  I  OJ  0  J  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ 


e . 


0* 

. 

c 

I 

I'M 


3 

CD 

C 

1 

> 

■b 

** 

3 

«■> 

a* 

a* 

c 

Ml  Ml 

«->  1 

C  1  111 

C  X 

O  o’ 

Qj  «v 

■->  «->  o 

0/  -h  r  : 

«->  d  * 

r-*  1 0  •  O' 

do  o 

N  £  X  ^ 

4^  ■  w  vs. 

N  CfOCO 

£Z  =  = 

O  1  • _ 1 

CD  O'- 

Cl— 1 —  *-« 

•H  ••  ✓“*. 

I  M  ItL 

r 

S 

2 

J—  r  (OV 

ox  «• 

-•0-0 

-  X  u. 

s  n  x  3- 

i  ID  O 

o-w  n 

O-  X  Ml 

UJD3  J 

n.  UJ  IX 

• 

CO_IUl 

1 —  CD  d  UJ 

y— 

TU.ZU1 

X  X  -►  X 

•H 

Ul  CO 

00  X 

UJ  O  d  O 

X  ►- 

X  c. 

1-  UJ  1- 

3-XOX 

CD  4- 

1 —  UJ  ZD 

UJ  XX 

II  - 

5- 

••H  0-J 

X  CD 

_J  UJ  1—  u. 

o  |—  I 

y— 

X  1  . 

I'M  *-■  O 

■Ts  UJ  C-  1 — 

II  .v: 

~r  of 

NCQO 

-CO  X  CO 

o 

/\  i2 

OQUI 

OXI —  X 

LlI  \jl 

o 

3 

M  T  tL 

XX  o 

OCCffiD 

_J 

►-< 

1-  X 

> 

3-  *-« 

X  1-  X  CD  > 

UJ 

,  Sw-  II 

»-*  CO 

d 

UJ  UIOI- 

II  XO  d 

3- 

♦> 

COXUI 

XIJM 

O  »  X  X  ■-« 

C 

OLU-I 

li 

y-  y-  co 

v  UJ  u 

o 

X  3  O 

O.CQO 

UJ  o 

CO  1-  X  Ml  Ml 

y- 

_J->  rr 

£.  > 

y-  u_xx 

/\  rilll  1^*H 

CD 

3-  3>o 

•O 

X  O  1— 

CO  O  CO  CO  — •  o 

Ul 

XX  CM 

ox 

1 

I'M 

UJ  »-«  UJ  O  1 

o 

UJ  Ul 

_i  u. 

V. 

/N^ui  - 

UJ  3-  X  _J  >v  U 

—I 

XX  o 

H-O 

d 

003-3- 

HUJHOO  d 

X 

l_»—  1 — 

LD  UJ 

CD 

vv*0  ■*> 

UXX><-  ID 

3- 

o 

xcox 

H-  1 —  XX 

UlCLOON-i 

3  -■  —  IJ 

1—  UJ 

O 

:  :  :  : 

Oi  r  I  HO 

3- 

O  w  w 

ac  cd 

v  Ml 

(=1 

X 

S£h-  X 

CD 

y-y-y-t- 

1-  1—  1—  iD  00 

_l 

II  CM  CM  UJ 

WCD 

3CLCL 

00X0 

ZZDD^D 

/-V  A  A  x 

xox 

CO  UJ  UJ  Cl  X  >~i  Cl. 

O  «-•  X  X  O  CO 

l - 

D-  •-«  .sc  .sc  —  I— 

O  UJ  UJ 

XXXX 

— I  X  X  X  N  o 

X 

V  •-*  c  C  1 

nmm  wma  >h 

h-t_IIHCJIJ 

UJ 

CO  3  3  •-« 

X 

3- 

CD  *>  ►-*  I 

•• 

1— 

X 

OX  •-<  II  /"*. 

c 

X 

o 

— *  ii  ii  in  ii  u.  u.  -H  H 

o 

UJ 

o 

||  X  i  -1  _ |MM  WM 

•H 

M 

3  ID  C  II 

♦» 

X 

O  •*«  3  X 

Ml 

o 

—  — 

JI1H  •— 

r 

3 

• 

<r 

1- 

•• 

i 

. 

o 

♦» 

u 

X  X 

c 

X  X 

Of 

o 

•-I  o 

3 

o 

1- 

C 

0> 

UJ 

o 

CO 

X - 

—  o 

O  O  O  G)  O  O  O  O  O  G>  O  S> 

O  O  O  O  O  i3  Q  O  O  O  O  O  O  O  ©  O  O  O  O  O  C 

^CN  CO 

—  CM  CO  *T  liO  CD  r- 

CO  |7>  O  O  J  CO  T  IO  ID  N  CO  (T'  Q  O  J  CO  M-  UO  ’D  O-  CO  C 

in  in  in  in  us 

'O  '.£>  •  X>  '  Q  *jD  '0  '0 

<£>  <  o  rs.  rx-  o-  is-  rv-  rs-  o-  rs-  o- 

CO  00  CO  CO  CO  CO  CO  OO  CO  0 

CM  CM  CM  CM  CM  CM  O J  CM  CM  CM  CM  CM  3  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  0 


97 


■  /  -  i  .‘j 


/% 

Ml  U 

►H 

1  U. 

V 

o- 

♦» 

- 

to  - 

f— « 

•H 

V  II 

Si 

u.  cs 

♦> 

_J 

1  Ui 

•H 

* 

A 

ct 

-Q 

/S 

CO 

— •  C  U- 

_l 

OJ 

1 

3  O  : 

UJ 

* 

▼H 

4J- 

Ml  ■» 

X 

/% 

i 

c 

<M  Ml  X 

H 

f— f 

t-t 

fltt 

(X  'U  *> 

*- 

V 

c 

\  K  -H  w 

21 

CM 

4- 

o 

ffOO  . 

o 

'w' 

CM 

a 

+  c  o  - 

+ 

+ 

X 

UJ  0-,<r' 

_J 

C 5 

«•> 

UI 

(SO  OMS) 

Ui 

c 

5- 

*v-x 

■  »-* 

z 

Mi 

to 

h- 

*-  : 

TH  >%.  H4 

z 

“•HI 

c 

M 

CM  — 

*  O:  X 

X 

♦> 

o 

o 

*  — 

Ml  C  II  •» 

X 

C 

a 

o 

t5  l- 

ty  X\  *H 

o 

C* 

in  x 

_l 

Ml  X  X 

Zj  CO  O 

Ui 

—•X 

*-•  ui 

Ui 

mi  m  *-* 

o  i?  x  o 

II 

II 

•«  X  o 

X  d>  X  —• 

z  • 

<r 

CS 

1— 

o 

*»  H-  Q_ 

O  h  O  'll 

ox 

CO 

oo  in 

T“ 

Oh  C 

Ui 

C  <S> 

1  IU.v> 

z> 

CO 

tl  f —  <rt 

UI 

tl  0/ 

1— 

O  .  II  « 

>>  II  >v 

ui  ac 

M 

2 

«■»  CM  C 

X 

X©  X  <x 

OX  CO 

U  1— 

J" 

Ml  ©  ■►» 

o 

c  — ■  O 

a. 

II  II  *>  UCl 

C  *>  CQ 

X  o 

X 

(d  N  Cw 

Q_ 

4i  ii  a,M 

4->  +>  '•*  : 

4i  •>■•  i  CE 

ac  ui 

ac 

•H  •— f  Qt 

X 

C»-*  X 

o 

<  •—  o 

3  O 

x  a. 

X 

♦>  SH 

Ui 

O  UI  h- 

o 

3  3  O  Q_  Q_  O. 

?  0  0.0. 

o  CO 

CQC  X 

Cl  QC  X 

1.0  Ml  •“*  Ui  Ui  CO 

Ill  >— i  CO  UI  u_ 

O  HI 

X  O  Ui 

«x»  fr  «l>  UJ  UJ  M 

U_  |'U  M  UJ  M 

—  — 

XU-  X 

—  —  — . 

Ui  U.  X 

— -  —  — . 

QC  Qi>  CQ  OQ  «  Ui  U-  >«  CQ  X 

■■ 

* 

CO 

pa 

♦> 

O  yN  u'  ^  (T*  ijs  OS  ON  O  O  O  O  O  Q  O  'S  O  *2hh  -h  ^  ^  OJ  '.\J 

J  CM  0J  O- J  oj  o.  J  O-l  0.1  r\i  0.1  0.1  ro  oo  o\  ro  m  i*o  m  i*o  ,v*  (-.->  i-o  co  oo  ro  oo  ro  m  ro  m  m  »v»  m 


i>x.  N  .  i^. r«  *j,  ' 


V: 

I 


►  **. 


I 


fi 

16' 


<p 

■H 

x 

i 

>* 


o 
o 
*— • 

>■ 

V 

3  S 
O  *> 

‘H 

i  o 

>-  o 


■ 

• 

«->  *“• 

E 

X 

•h  Ob 

ZD 

ZD 

O  3- 

cm 

cm 

O 

t— 

h- 

*-•  o 

o 

O  LU 

a>  x 

LU 

LUX 

3-  X 

/N 

CL 

X  o 

+ 

cf 

CO 

CO 

*  :>  ^ 

♦> 

o 

1  o  -• 

0- 

UJ  1 — . 

d  i  + 

-O 

I- 

»— 1 

X 

«->  1 

1 

»H 

'w' 

t-Uh 

r-4  >vM 

•-t 

o 

«-> 

XX 

0»  *ov 

V 

o 

IDE 

2!  »-H 

O  -i  *» 

♦> 

-I 

00  IT 

ooo 

*  OX 

x 

LU 

tH  -H 

X  10- 

H  O  f 

ip 

3- 

0b  : 

H-O 

V 

ox  ■ 

XXX 

>*  to  to 

Ob 

LU 

H  ^  - 

U  CO  1— 

*»3-X 

x 

I 

»-H 

►I  UJ  X 

i  /\  II 

■* 

h- 

Q  _I 

ujx<=> 

O  > 

cf 

ii  LU  LU  *1 

X  X 

O  *»  n  n 

♦> 

? 

— 1 3-  ca 

o  > 

—«-iv  || 

CC 

OXI- 

UJUh 

to  O  ♦>  tf 

-o  . 

cm 

QCX-IOC 

T  »H 

~  3-  O  X  -» 

Is 

a 

O  LU 

l-LUO 

1—  +  — ■  T  C. 

M 

•IL  X 

3-  O 

X  i  _J  to  iTj 

_l 

t-  cm  _i 

LU  11  n  3-  1>  1 

LU  Ui 

LU  ~ >  LU 

O  Q_  n  1 —  X  n 

3- CD 

mm 

00  3- 

II  LU  XX 

o  x 

£ 

—  — - - —  — 

HH  C£  2*  »-» 

z  _i 

3 

1- 

M 

«•> 

«•> 

O 

X 

0p 

IP 

a 

•H 

z 

irt 

X 

m 

X 

i 

cm 

X 

ZD 

3 

ZD 

1 

H 

O 

\— 

♦> 

LU 

i_ 

LU 

OP 

x 

« 

cm  — 

- IJ 

O  Q  G)  O  CO  O  Q  CO  O  O  CD  CO  CD  CD  O  iO  O  CO  CD  CD  CD  CD  CD  CD  O  O  CD  CD  CD  1 

co  r^-  co  t.  o 

— <  o  j  co  ’T  in  '  U  r  -  co  it>  o 

— <  O  J  CO  T  liO  'O  P- 

CO  «T>  ‘3D  »-«  O  J  CO  Tf  1 

OJ  OJ  OJ  CVJ  CO  CO  CO  CO  CO  CO  CO  CO  00  CO  ’T  *T  ’T  *T  rf  rf  rf  rf  tt  TTIO  ID  ll-'  UO  llO  1 

ro  ro  ro  ro  ro  ro  fn  co  co  ro  co  co  ro  ro  oo  in  co  CO  CO  CO  oo  CO  CO  CO  CO  CO  co  CO  CO  < 

LU 


X 


x - 


99 


*  *,*•»'  .  •*/»  ,*«  IV’/f  /«  -w  >  ,%  „>  v«s  •  .  *  ^  •  3  "  «  •  •  •  •  .•  -  . 

-ia*-""'‘ - - -  •■*  .-•■-^V-V-V-V.,  /  V -VA 


kM  ^..1, ,  in 


s-  •' 

*  *  1. 


=5 

3- 

UJ 

o 

(X 

1 

2 

•  on 

J- 

cf  s 

1— 

>-  UJ 

o 

♦>  A 

CCIVI 

UJ 

_i 

in 

Q£ 

a. 

CD  UJ 

CO 

Q£Q£ 

in  in 

05 

1=1  2 

C 

■W 

CO 

'  .  - 

22 

-*■« 

V 

■  u. 

+  +• 

5- 

:  2 

S 

O  OJ 

•  CO  /\ 

)')'>>>' 

1— 

(=1  2 

1 

V-  \ 

2  - _ 1  ^  LlI 

I  1  1  1 

»H 

C=l  O 

23 

4»  X  > 

UI05UII— 2 

>■ 

X  C  C  X 

o 

C5 

UJ  2 

N  1  1 

UuJfljXw 

n 

X  >*•  &  -H  *H  d  Z'  Z' 

o 

-  OH 

2  7* 

3 

1  C  X 

CXCL3_J 

f— • 

1  iTrrr  |  f 

_l 

i=l  UJ 

»—  ■-< 

L 

P—4  -r-#  d 

Q.  X 

cx^^^^xc 

UJ 

(=5  0- 

in  ac 

♦> 

7*  7“ 

CO  UJUI 

^  1 

13  min  in  in  rs-*« 

> 

•fc 

0(02 

o 

<=»--. 

111  UJ  XX 

■-  c 

T^-1"  ^  •  •  •  m  >_  Zl 

4t=  '"N 

II  —  II 

4> 

2  0  0 

mix  t-  • 

XJ  •-< 

•>■  *»  1  !++"■•' 

3=** 

=  05 

L  z\ 

a. 

■r  •*-  u 

XI-I-UJ  >- 

l 

xxxxxxxxtno 

\ 

U> 

<D  (if 

1—  I>l> 

V_ 

1  1  1  I  1  I  1  1  ON  _| 

UJ2 

OH-v 

1 

O  N  N 

2  x  y-  ac  •-* 

CJ  X 

xxxxcccc  -• 

■n  33  O 

*-•  «-> 

3 

^  1  1 

COO  JO 

Ol  1 

O  i5  r3  >3  -<■*  in 

r 

:  •  v-z 

0  0  2 

(3 

O  M  w 

X  cuo 

i—4  C 

V-  05 

:  ;i  p*  *-< 

— 

1. 

IU  1 —  1 —  Q£  05  _J 

o 

s  s 

s  X5J 

» 

to  CU3 

— 1  X)  Z>  UJ  — •  Ul 

X  CD  H-  Ll_ 

N  :>  <2! 

OQ.ClNQ> 

iXi 

CCXCCCCO 

_J  _l 

_ l  — 1 1 

OCIXwCO 

IOC 

=5  UJ 

acexcxiirciiLXcxOiOLLiuJUJUJUJ 

O  O  UJ 

X5 

*-•  s  « 

05  > 

«(=ii^«i=ii=iC=ii=i2 

CD  CO  CO  CD  GO 

U.  2 

C=i  05 

i 

-«  oi  CO  ^  (O 

o  o 

■r  >2  a:  az  a: 

u. 

2  O 

u. 

—  —  —  —  —  —  — 

IJX- 

_i  _i 

_ 1 _ 1 _ 1  — 

M 

Ul  UJ  — 

>-■  1 

c 

to 

to 

i. 

o 

(rt 

l 

3 

to 


O' 


UJ 


SOOOOQOQOOiSOOOOOOOOOO'SfflOOOOOSOOOOiSQOOO 
O  *-*  OJ  05  -T-  IT.'  '■£>  fv.  CO  0N  C3  *•«  OJ  CO  "I"  U5  'jj  f'-  CO  '7'-  O  ■— '  01  05  ’1"  to  '  jj  f'-  05  ON  ©  OJ  05  T  U5  'i5  fv. 

co  co  *o  1.0  «.o 1  o  co  '-o  'is  co  n  o-  in.  r-_  o.  o_  r-  o-  r--  r'.  co  co  05  co  co  co  co  co  05  co  on  cr.  on  ■?>  ctn  on  ctn 

CO  05  CO  05  05  05  05  CO  CO  CO  CO  CO  05  CO  05  CO  CO  CO  CO  CO  05  >55  05  CO  05  CO  05  05  05  05  05  05  05  05  CO  05  05  05  • 


100 


•'  .*»  ,*•  . v  v-  k.  »  V 


►  A  -V 

•  4  *  • 
A- 


X 


l!) 


• 

3 

• 

01 

l  • 

1— 

X 

a_  >-  x 

w 

CD  o  o 

o 

OQ 

wxo: 

o 

O 

O  111  1 — 

_i 

CD 

oo 

LU 

O 

►  Old 

> 

UJ 

:  LU  CL 

X 

UJCtCCO 

UJ 

LU 

X 

OIL 

x 

X 

LU 

►-«  LU 

t-  ■ 

h- 

X 

ooi 

1- 

XXI- 

Zh 

X 

o  ac 

— <  x 

cf  ^ 

CD 

1  X 

M 

£  ►- 

CO 

— •  >-  o 

X  o 

1  v 

A 

:  1— 

oca. 

T  *> 

♦>  ^ 

»-*  H 

LU  1 

Av  x 

C  X 

1-0  X  • 

Q-X 

lO  ^  O' 

X  cs 

XC'wiO 

O 

00  •-<  •** 

«J  O  £ 

H  JOO) 

UJX 

£  a  i 

O  LU  CL  O 

X  CD 

*»  ►-«  X 

i  i  >— 

Q->  CL 

1—  : 

OX  II  II 

i— <  3  '-/• 

LU  O 

1—  IT  X  X 

II  O  •» 

LU  LU  X 

«  _l 

-r*  &  Cj  U_ 

ll  X  X 

XXI —  X 

X_l 

© 

O  to  £  £  *-• 

O  CD  7* 

h- 1 — 

—  <n 

II 

11  X  1  1  ►-< 

•H  -H 

X  X 

Ll  o 

5< 

•— <  *HlQ 

Old  -U 

QC>  O  1- 

Ll  XI- 

xox 

O  0C  X 

■  • 

£ 

CL  LU  X 

O  CD 

Ll _ 1  05  3 

—•  OJ 

1 

O  LU 

1  O  Ll 

X  — 

Ll  X  — 

M  iJ  M 

OOOOOOOOOCD  OOOOOOQOOQO  *S>  ©  ©  Q  O  ©  ©  C5  ©  O  >3  CD 

CD  cr>  ©  -*  >"‘.J  CO  T  U5  *D  O-  '55  >7'-  ©  — «  OJ  CO  T  U5  VD  T'w  CO  >7'  ©  —t  OJ  CO  T1‘  U5  >■£»  O-  05  O'1  'S' 

O'  O  0  O  O  O  O  O  O  Q  O  H  ^  ^  -1  Cl  OJ  OJ  OJ  OJ  O  J  OJ  01  OJ  OJ  05 

co<0"frT-rTT"r'<r"r'r’r’rT''i-TT',T->TT’j-Ti-TT'TTi-''3-Tj-''3-T'-j-T'i- 


1  •  ■  •  ^  ’  i  *  »  •  < "  i 


.  V„%  ,*•  /.  .%  A  *.  .%  .V  *.  %  •,  S  n,  n  •,  *  •,  .  -  •  -  .  . 


o . 

'■XLi  s.’. 


4310  V=Height  <  I. choice) 

4320  MOVE  KtY 

4330  LABEL  Mt" 

4340  PFH  1 


w 

:  C 

/S  O 

I 

•  *H 

1 

to 

k1 

CO  V. 

•  » 

Ok 

x  to 

H* 

X  > 

«• 

o  c 

s 

V  o 

• 

D 

o 

A 

t' 

:  X 

s 

U 

•w  »H 

■k 

■ 

•— « 

P 

to  :> 

to 

3>*» 

CO 

.  £ 

it 

o  • 

O 

:  <3 

N 

fH  : 

»— « 

<S>  X 

i  H 

X 

in 

\  ^ 

N 

ii 

a_ 

N- 

UJ  : 

**  3- 

X 

n  v 

>il 

•»*>o 

X 

in 

HI 

■ 

O  1-  X 

u 

in 

>  lu 

If  ^ 

o*-ui 

OJ 

—*  -J 

•7  to 

O  X 

o 

•  vt 

<=*  ;  — 

•J  u 

to  O  O 

X 

Z 

:  X 

.  * 

I>  — 1  UJ 

•X 

ca 

»  o 

+  UJ  X 

m. 

1-  >-  U 

x 

:>:>x 

X 

** 

X  <3  X 

o 

1  V  V 

X 

: 

H-C  r-«  HI 

*;  i 

Or  s 

X 

X  Cl  U 

£  to  •— « 

4^  »v  n 

_l 

u 

X  to  X 

o  '-s 

•— «  ^  /\ 

X 

X 

--•X 

W  -W  *i 

to*. or- 

XX  "OX 

r*  ox 

P  ^  * 

UJ  • 

X  CO 

Ll!  X  1 

N  x  v 

*  CO  CO 

X  CO 

O  X 

• 

X  X  *-  X  4k 

w.  o~« 

©  V 

1-  >- 

CO 

i-  h-  o  a:  to 

k.  i  o» 

o>> 

X 

CO  CM 

X 

to  £ 

‘  wlio 

-r*  3C  X 

X  _J 

O 

•X 

Tj-  fH  H*  ^ 

w  II  II 

CO 

O  CO  LLI 

X  X 

w  1—  X 

X 

co  in  o  =  x 

X 

k  -w  X 

x  a:  a: 

UJ  CO 

XXX 

X 

II  II 

o 

•  «  i  i  in 

0 1 —  1 — 

—I  »-• 

X  X  UJ 

UJ 

xxtm-ui 

X 

■  a  a  cm 

i 

o  o 

XI — 1 

X  X  X  X  IX  X 

X 

:  a>  a> 

•— « i~  i — 

x  x  o 

X 

IV  lY  ITI  H  tV  M 

X 

t  1  QJ  CU  LU 

n  y~r 

—  —  —  — 

■•  u  o  u 

1-1 

X  LU  O  X  X 

;  .■S4.it.->- 

HHH 

X 

uxxn 

X 

o 

to  x  x 

XX  X 

X 

» 

xZ 

3=-  O.  u. 

—  —  — 

..  KH  HH  LLi  — 

o 

| 

<3 

mm 

a 

to 

•H 

o 

■o 

w_ 

X 

1 

X 

" 

X 

x 

i- 

X 

'll 

X 

"Z> 

<3 

X 

1 

X 

\- 

►— 

X 

h- 

UJ 

c— « 

UJ 

•H 

U1 

X - 

- o 

X - 

- LL. 

X  - 

i  O  (3  Q  O  5)  O  O  O  Q  O  O  O  '3  O  O  O  O  O  '3  Q  Q  O  O  O  Q  O  Q  Q  O  Q  O  '3  ® 

in  •■i)  r--  co  >7->  o  cm  co  t  in  u>  o-  oj  co  ©  »-«  cm  co  ir  in  o-  co  co  ©  *-«  cm  co  t  in  u>  r-- 

co  co  co  co  co  m-  m-  m-  m-  m-  t  ’*•  m-  t  t  in  in  in  in  in  in  in  in  in  in  1  o  '  i:i  co  <  o  co  *o  co  ’  O 


1PJ2 


4680 

<1  £  Q  Cl 


I  II ■ I «l ■  I 

•  .  *  ■  « 


i 

* 

4 

i 

j 


- 

! 


I 

I 

« 


t 


i 

i 

t 

t 


i 

I 


C“ 

CO 

X  X 

z  o 

o  c 

v  fir 

3 

>-  w 

O  >1) 

Z  x 

UJ  u. 

©  — 

O  5 

u> 

UJ 

Z 

QL  II 

•H 

M 

u. 

UJ 

<=>  ■ 

X 

i.O 

ZZ1  2=- 

O  N 

i— 

i 

>-o« 

1-  x 

S’ 

h-  U _ 1 

•S^ 

Ul 

0; 

M  CD 

» 

O 

*> 

O  Ul  LU 

UJ  > 

u. 

oocx 

1-  CJ 

y- 

w 

DC  Z 

O 

: 

UJ  o 

UJ  UJ 

z 

o- 

>oq; 

3-  © 

UJ 

>* 

o 

z  o 

© 

o 

ZUJlL 

O  LU 

Of 

UJ  x 

O  O' 

Ul 

Ul 

c  u_ 

oc 

ZZf 

HI  (J 

Ul  o 

Li.  • 

Of 

u>z 

CCt  •**  UJ 

CD 

Ul 

coin 

X  Ul  Z 

UJ  — 1 

X  ol 

>N 

<r_i© 

»3  o  • 

X  _J 

cS  Ul 

CS 

0.0 

-<K  4(0 

h-  UJ 

r—4 

Ct  00  UJ 

li  3  Z  C 

o 

1 1—  c 

a 

cm-io; 

ui co  'if 

UJ 

lOU.Cu 

U.  QU. 

■rth  OO.  01 

ID  ID 

•H  HI 

•H 

"7J  *— <  O  CD 

Z  id 

•OX  u 

T> 

IU\UJ  — '  O  © 

X  X 

1  CD  O 

| 

V-OCC  II  ul  II 

X  O' 

■-  m 

V. 

•• 

cf_j  ii  <z  -  i  - 

O  CO 

CS  Z  1 

CS 

X 

Ol  UJ  X  O  <3  3  C# 

Ol  Ul  3 

•ll 

o 

H>  OO  0l  Ol  111 

—  —  — 

—  rH  ^  Cy 

r-l 

c 

Os  C  s  '  Z  • 

O  S  Z 

o 

Oi 

'1>  o  o 

»• 

3 

COb-  31— UJ  ICO  1 

c 

CO  1-  0Q 

CO 

# 

©  ©  *ZCD*s-0«->- 

01 

©  ©  © 

© 

to 

CD  O.  dli— «©  CCD  C 

L. 

CD  a.  05 

CD 

OZ  V-Qil  OO  O 

o  z  o 

O 

iJulULfl-OtJO 

1 

LJhiL3 

—  UJ 

1 

Ol 

0£ 

(7* 

ol 

o 

© 

c 

© 

♦J- 

K 

cf 

1- 

•w  (=1 

c 

UJ 

1— 

Ul 

3  Z 

o 

o;  — 

- o 

OL - 

- ©  UJ 

JOOOOOOOQOOOOOOO'SOOOOOOQOOOOOOQQO 
%  o  O  J  CO  T  UD  '  £>  r  --  CO  O'  ©  — *  O  i  CO  -f  li'i  '£>  T'-  05  IT-  ©  — <  I*  J  CO  T  lO  CD  fv  CO  <0*  ©  — ■ 

)  r-.  f'-  p*-  Is-  f'-  r--  r^-  r--  Is-  p.  co  co  co  co  co  co  co  co  co  oo  o>  o>  O'  O'  O'  O'  O'  o>  o>  O'  ©  © 


103 


3 


330  CONTROL  So  35  Baud 

340  CONTROL  So4  5Bits-5+4  !  BITS-'CHfiR  &  #STOP  BITS 


CO 


Cl  **t 

c  or 


cl 

3  X 

LU 

73  O 

df 

Ll 

1  **4 

73 

Li- 

o 

O  — 

3 

♦->  LL_ 

s:  d 

s  CO 

d 

1  c 

+ 

1=1  y- 

<u  ••■<  e. 

LU  _J 

o 

C  €  v  ♦>  £ 

O  _l 

, 

o  u_ 

•*•«  l  i-  i-  3 

OO  c 

H  H- 

_l  O  C-  73 

21  Li- 

oo 

1—  1  1  1 

O' 

L3  O 

O  L  i-  d 

x  cl  'ii 

CD 

1—  O  Q_  Q_ 

OO  CO 

L'XX 

0 1—  d 

•— «  Li_ 

4*  LU  LU 

U  O  O  O  <=t 

1-  o 

CL  X  X 

Uh  h 

a.  iL  i— 

XI- 1- 

=  OO  o 

UJ  <_)  O 

o 

Cu  =  DOh- 

O  LU  •  J 

•n  0 

73  —  O 

LUI'rt 

<r 

:  II  II 

O  d!  s  O 

iucz 

■  ■*. 

<x  a  l 

z:  c  u  u 

Ml  LU 

co 

•c  e  0) 

•h  L  L: 

lu  ~  c  x 

«*. 

4t  3  --- 

<n  oi  sroiii 

X  d  1— 

o 

•  73  3 

■ 

•W 

•-  C  -  |— 

-<>  c 

* 

CO 

d  a. 

LU 

C 

Ot<  4i  O  OI 

— I  ••v  1—  Q 

LU 

L3  *»  X 

0 

CD 

5w |  | —  1—  |—  (=1 

•«=  ©  II 

30 

CO 

X  d  O 

H-4 

\/ 

L:  :  :  :  : 

r  -h  O 

h-4 

00» 

►-<  (=i  O 

I— 

II 

LU 

■•4  f)  Q 

1— 

i- 

CO 

0» 

_ 1 _ I _ 1 _ 1 _ I 

LOO- 

3 

■X 

OOQ 

o 

>■ 

O  LU  UJ  LU  UJ  Ld 

—  u.  CO  L3  3- 

o 

1— 

0  01  *0 

CL 

Cu 

1—  CO  CO  CO  CO  CO 

s-/ 

CL 

O) 

L  >X  'X  •-« 

o  >r  <x  >x  a:  x 

i —  co  ca  i— 

lL  rr* 

1— 

o  © 

L3 _ 1 _ 1 _ 1 _ 1 _ 1 

rociiw 

LU 

COCO  >1» 

Kf 

©  II  II 

Cl  1 —  iL  CO 

05“ 

H  V£»  ^  CD 

||  Q.  — <  *»  s 

CL  ©  If)  CO  T'-  CO 

H  T 

►H 

O  II  II 

■30 

=*  £  ii  a<  c 

O 

<li  X'  1—  X  Ll 

LU 

Cl  CO  CO  O 

30 

mm 

(fl  j  C>  Vl- 

CL  3-  3-  3—  3-  3- 

73  O  CO  O  •— < 

O 

«« 

Ih  1— 

•X 

•<■«  73  ♦»  1  II 

lV  III  ill  ill  ill  111 

O 

LU 

■ii 

II  O  U_  Ll  O 

CL 

•H 

7)  c!  d 

LU  1L  ii  V  iL  M 

■  ■ 

CL 

>  CQ  O  »-•»-«  IJ 

1— 

1  ^  r<  »H 

1  c 

•H 

Ml 

O  d  1  d  73  II  II  II  II 

/  33  c"  rT 

(If  -H 

fr 

c 

O  O  O  O  O  O 

- 

— 

—  o 

— 

— — 

••  d 

•H  ill 

Oi 

1- 

_l  CO 

CL 

1 — 

Q  O  Gi  Q  >3  O  *2)  O  O  O  Q  O  O  <25  O  O  i 

0  0  10)  10)  <01  10)  0  101  0  10  10  10  10  10  I0i  »0  101 

©  ©  © 

r>_  os  a  >  ©  -«  cm  co  t  uo  co  r--  oo  >?•  o  -*  cm  co  t  uo  co  r--  co  o  © 

cu  co  *r  uo  co  f'-  co  o- 

©  OJ 

(*)  CO  CO  'T  M-  T  T  "I-  T  "T  7"  "f  T  lO  ll‘J  UO  UO  If)  UO  UO  UO  UO  LO  CO 

•-D  VD  VO  '•X)  '-£»  Vi)  VjD  h-  Is- 

105 


V  ",  >  "  •,  •/  ■'y  «*j  S'. V.V.V  '  V  V  V  V  V •  '  »  "  .>  .  > 


CE 

r*- 

■fc 

OJ 

VO 

: 

v 

•w 

m 

Vt 

X 

x 

nr 

X 

»—• 

o 

CO  - 

II 

II 

CD 

V* 

V* 

= 

Zk 

>• 

V. 

<1^ 

a> 

OL  i- 

X 

X 

h- 

"7" 

32! 

CO 

. 

G  1- 

Ul 

Ul 

1—  1 

X  G 

X 

X 

Ui 

V 

0# 

1—  X 

1— 

G  • 

1— 

1— 

4* 

*— • 

G  l— 

UJ 

X 

•H 

OG 

oo 

C  X 

00 

_ 1 

X 

<♦- 

c  CO  = 

O 

II 

c;  m 

II 

Ui 

CD 

i 

■H  •  || 

X 

/V 

c=t 

11 

m 

c 

5>  .  ✓*. 

XX 

V* 

irt  '!• 

** 

* 

X 

*1? 

■D  h~  X  s 

X  LU 

X 

O'  X 

'X  UJ 

>•  >. 

X 

G 

X  >—  X  UO  •»• 

t-  X 

'll 

I?  4-' 

■1' 

G 

H  CU 

1 — 

O 

X*-G  o 

1- 

Z*Z 

•*-« 

IX 

CO  X 

iiZi£ 

G 

OX  .  u>  X 

UO 

V 

IT*  V- 

V 

UJ  X 

•»w 

X 

X 

1 —  3  -  O'  d  X 

in  in 

7=~ 

O 

X  00 

T  : 

SI 

G 

OX  O  G  G  ^  X 

04  UO 

'll  <*■ 

G 

»— • 

Ul  X 

© 

CO 

L9  LLl  »-*  CD  1  cl  L  (- 

AOJ 

X 

3 

G  O 

X  - 

G 

o  . 

X  vlUO  — •  G  'if 

V  II 

•r-«  Ct* 

GX 

y-  & 

G 

XI-  '-*•*•  .  1  G  G 

/■>.  A. 

o 

—  o 

i=) 

UJ  X 

z 

1— 

1 1 1  *ii  :f*  r  *— «  ^  O 

X  +*  4* 

rr 

■3 

X 

X 

UO 

OO 

!2I 

2C  ■*“*  j3  P  v  3  C**i 

Ul  >V  >v 

'X 

*n  a 

X 

X 

U0  CZT 

o 

X 

1—  II  *>XXcSC  1  —• 

X  <u  <u 

3  vt 

'00 

04  X 

X 

X 

v-  XI  'X  1  II  II 

1—  Zs<£  Zs£ 

•»H 

O  -ii 

UO 

X  X 

II  *-< 

K 

(S  *X*  •*-  +  O  cJ 

V  V 

II 

O 

II 

>  o 

CO 

X 

ii  ♦»ocsL.Lce*t 

®  XX 

3  cf 

X  X 

*fy  G 

X 

G  3<*-GXX  '1»*J-  C5 

II  GG 

<D 

■if  jU 

CP 

oc 

c 

1— 

1 

£  G  G  CO  c=  G 

•XXX 

4-> 

X  _l 

4'  O 

a 

G  X  •  i —  1 —  G  1 

l  4/ 

3 

0. 

3 

3-  O 

'00  IT* 

o 

♦> 

G  O  -■-<  G  G  X  1  *-*• 

i  <-«  U_  U- 

a 

g 

X 

V  ft* 

J— 

CS 

CSO  o  X  X  X  G  *-•• 

;  Cl  ^  *— •  Lu 

z. 

«1>  V- 

Ul  X 

X  1—  X 

G 

*>  3X1 —  1 —  O)  v 

3 

M 

o 

X  o 

o 

X  0- 

G  G 

1— 

1 

rS  U_  LU  ID  ID  O  O 

(=l 

o 

+>  <4- 

CO 

t—  00 

X  X  O 

G 

M 

fiw-LflQO  G  'X 

. 

« 

1 —  1— 

o 

U- 

XU- 

Lu 

X  G  O 

X 

X  X 

I— • 

ti|M- 

- - - — 

»•  1— 1  — 

■•  —  —  - 

■•  H-t  O  U 

X 

—  M  — * 

1— •  H-f 

GOOOOOOOG3 

O  O  O  O  O  Q  O  O  O  O  ©  O  O  0  O  O  O  ®  G  CD  G<  G  'ZD  'ZD  'ZD  CD  CD  G  *<  04  CO  -r  UO  'X*  f'-  CO 

CO  rr  UO  CD  f'-  CO  <X>  CD  — *  04  CO  -T  UO  '■£>  f'-  'jO  CT>  G  —<  04  CO  "T  UO  VO  {'-  CO  CT*  G  'ZD  G  G  G  G  G  G  G 

o-  (-.  o-  o-  f  -  0-  co  co  oo  oo  oz>  co  co  >zo  co  g  <r»  tr»  or*  oo  cr.  o'*  cr.  cr.  cr*  cr>  »-*  »-*  *-*  ■>■* 


ZU'-- 1~*.  ■■ 


r 


li 

£ 

3 

3 

I 

cf 

cf 

« 


I 

A 

+ 

■ 

LU 

1*C 

dc 

CO 

LU 

H- < 

:  J- 

« 

/V 

2 

•  D- 
T>  a 

O 

a<  o c 

1— 

LU 

o 

C5  lL 

1 - 1 

4>  X 

CO 

C 

111  A 

1 — ( 

a*  i 

: 

■  ww 

mm. 

:  -W- 

¥t 

•  ON,  0,1 

4' 

X  — 

cc 

•H 

a:  <*- 

<*- 

LU  1 

1 

mm*  (2. 

0. 

:  X 

X 

4t  -- 

Mw 

CD 

L  V  V 

H« 

o 

CO 

-  LU  LU 

o 

•-  _l  _J 

L3 

0*  »— «  1—4 

■ 

VU.  Li- 

CO 

“7* 

S  1=  TJ 

o 

LU 

•  ml  Cu  : 

HH 

X 

*•  «»l  «-> 

i— 

H- 

CL  L-  Cf  L. 

CO 

Cl  'JJ  Cl 

*— € 

1-  L- 

DC 

£ 

2>  1 —  O  1 — 

LU 

•—4 

Cl.  3  ID 

H 

:  3 

HlLOQ. 

o 

II  II 

LU 

XL  I —  1 —  I — 

X 

4*  li 

*7* 

ODOD 

X. 

3  £ 

+ 

►—4 

CD  i-D  CD  * 

■li 

X 

X3  3 

cf 

f~ 

T 

O 

X 

_J  -u 

♦> 

ID 

UZIIO 

•H 

o 

cf 

CD 

X  LU  LULU  h- 

to 

O 

♦> 

■3 

Lli 

h-  X  X  X 

T5 

X 

_l 

Lu  cf 

1 

CQ 

1 —  1 —  1 —  'li 

O 

■ 

o 

X 

>-»  l=l 

»— 4 

3 

Oft  r— 4 

£ 

DC 

1 

D 

CO 

ir>  -4-  <jd  m-. 

| 

LU 

i- 

H4 

li- 

Cf 

qc  AinifiioiL 

0^ 

1 — 

4' 

H-4 

4-> 

L3 

O  V  II  II  II  ® 

■ 

c 

22 

X 

n 

•~y 

TV  "T"  ~7"  ~7“  "7* 

1— 

»-• 

i— « 

LU 

ou 

*— i 

ll:  q:  lii  ui  m 

DC 

_l 

DC 

•H 

1— 

X 

i 

_i 

lii  frT  iV-  iV  i'K  l_H  I'l* 

LL. 

L  111 

LU  H-f 

a 

LU  LU  LU  LU  »-«  7? 

o 

Cl  3 

LU 

"7“ 

U.  OD  O 

o 

II  1— 

o 

II  O 

X 

— 

X 

U_  li_  li_  li_  li_  0D  £ 

1— 

L  O 

f— 

L  £ 

1 — 

c 

X 

O  ' IHMW  X  1 

CL  13 

CL  1 

-W 

0; 

1— 

f— 

Cu 

LU 

■  ■ 

O' 

DC 

c 

Zj 

3' 

r 

L3 

rH 

O 

*  4 

CL 

M 

Q_ 

■  ■  •  H 

T?“ 

cf 

CD 

DC 

•it  —1 

f- 

1— 

4->  -J 

X 

DC 

V- 

3 

3 

1*- 

X 

•H 

O 

LU 

o  o 

o 

o 

'3 

D  O 

•3 

H- 

L-  f— 

1 

1  f— 

O 

L  O 

Urn. 

L  O 

4> 

IJ - 

>•  -—•  — 

•  LU  13  — 

<—  — •  — 

..  a. 

—  —  *  — 

•«  o.  I_D  — 

—  *  — *  — 

•  1 — 

3  O  O  O  CD  O  O  O  O  O  O  O  O  O  '3  O  O  <3  '3  CD  >3  3  CD  O  CD  O  O  '23  O  O  CD  O  3 

Js  Q  «— «  0-J  (V)  rj-  UD  *C'  0-  CO  O'*  3  mm  O  J  i'm  tj-  UD  *X’  F*w  CO  O'  <3  *— »  OJ  CO  TJ*  UD  'C1  f  --  CO  O’*  3 

3  ^  ^  ^  ^  ^  OJ  W  O  J  O  J  O  J  OJ  OJ  OJ  OJ  OJ  CO  CO  CO  CO  CO  CO  CO  CO  00  CO  *T  T 

H  m  I  ^4  ^>1  ^—4  «^-4  ^-4  <H  ^-4  ^4  ^4  *■  4  «*-4  <p4  ^4  ^-4  ^4  ^4  »H  ^4  ^  -^4  v4  ^4  ^4  ^-4  *H  »H  »H 


107 


^ »  •  «  '  ■  ' ' 

<A  Of  tiT sZ  V  '.-■.  V.  j 


<■.  -’■.'  -\  *  I  ■ 


1420  OUTPUT  P r 5 "  1.  Baud  Rate  ="»Baud 

1430  OUTPUT  Pr> "  2.  Bits/Char  =  "5  Bits 


4l 

X  •>. 

to  3  CVI  rr 

—•  a  ii  ii 


a i  £x  s- 

♦J- 

m 

om  id 

Oi 

• 

to 

+* 

n  <*-  o  3  x 

X 

Ml 

mm  mm  J  /w 

•H  3 

£ 

¥► 

1 

O 

r  cf  r  \  O 

"O  i/I  — •  O 

O 

3 

c 

-sr 

II  X  II  X 

X  3  O  £ 

TJ 

*H 

ii  a:  - 

O  cf  — <  3  O 

<4- 

—I 

V*  v  +> 

0 

o  >  co 

■HfijCOPU 

aw 

p-  to  to 

♦j- 

-  i-  ii 

*T  n  n  M 

r 

cm  to 

cu  ~  Cj 

* 

to  -  1 

a. 

x»-«  ^  CO 

-  C-  O-  O-  C-  II  II  II 

Cf 

CO 

£i  x> 

3 

013  3  11  1 

:  OOOO  X  X  CO  X 

ZJ 

IxJ 

-rr<>-  to 

hvCE  CO 

O-  )—  1 —  I—  f—  to  to  II  to 

>- 

i-  cf 

■r- 

a  II  £  CO  QC 

to  =  s  s  :  •-«  — 1  ift-t 

■O 

11 

to  _l 

O  to 

3  -«  O  —  l- 

C  d  il  ♦>  t 

n 

a-4 

O  *»  ■>• 

1—4 

PljOlj 

OHHI-h  3  3  •<-«  3 

to 

*= 

X  to  r 

*^4 

D3DDPPCDO 

Ul  to  O- 

3 

X  3  3  3  3 

o- 

-  4C 

O  1 

•  • 

O  Z  !Z  21  !Z  121  HT  21  X 

0 

Cf 

r  Cf 

O 

oo  rr 

•a*  w  whh «LlILULlIUJ 

+> 

O-  «J-  ■<-1 

to  3 

X  3333 

cf 

#  <3  to  I 

■  4  vw 

333331—1 —  1 —  1 — 

to 

Z) 

ZJ  to  *>■ 

•H  I 

Ixl  Ixl  Ixl  Ixl 

iH 

*>  to 

<#. 

s 

to  3  3  3  3  —«  00  CO  in 

H 

:> 

C  (P  Cf  Ml 

to 

O'  )—  l~  ►—  1 —  II  II  II  II 

<4- 

a 

to  C  ■*■>  l 

C$  £ 

C  ■*.  i-  i_  l_ 

IJZ 

i-H  O  C 

♦J-  C5 

rj  O-J  CO  Tt"  to  to  to  to 

cf 

UJ  Ul 

>H  L  "U  *H 

c:  i_ 

i-  *-  1-  W-  *-  W-  X  II  II  II  II  *>«»*>  *> 

a 

3  m  3 

•-  C  Cn 

0 

a.  a.  a.  q.  a.  q.  o  x  x  x  x  3333 

C  PC 

K  *h  h* 

Q|  *H  4J-  Of 

:  to 

5  OOOOOOOO 

to  1/1 

O  O'  Ml  CO 

a  m  r— 4 

1 —  1 —  1 —  1 —  i —  ►— 

•H  *H  *H  .■”  t.  ^  £ 

Cl  ,H 

O  : 

X  to  0  11 

a.  -H 

OOOOOiOH-XXXX  O  O  O  O 

0  31 

||  —1  IxJ  CO _ 1  *>  s 

Ol  U_ 

0.0.0.0.0.0.3)333:3:0000 

+*  Irt  = 

r  2  =  Of  aH 

s 

1 —  1 —  1 —  1 —  1 —  t —  Cl. 

to 

to  <—  11 

Ml  = 

1— 

3D33D3 

ZU.U.lLlUxlLU.U.3 

— 

—  —  — 

■•  II 

Ifl:  -Mr  f —  1 —  1 —  1  II 

0  1 — 

OOOOOO 

M  M  M  M  »H  1— 4  ►— 1  1— 4  1— 4 

O 

O 

1  3  0  O  O  Q 

3  0  1 

£ 

SI-UO-iLLL^  3 

1—  3 

1 

mm 

3 

r>  _i 

333  0  3  Ixl  O  3 

, 

to 

to 

•0 

e  3 

WMMp  JIJJOw 

C 

1—4 

<3 

0  3  u. 

_l 

•H 

•H 

H4  1—4 

Ul 

i 

_l 

<•- 

CS 

u. 

1 

0  — 

M 

OClOQSIOOOOOOOOOOSIQOOOOOOSOOOOOO'S'SSSIOOS) 

rr  in  'D  n  co  a-,  o  — «  oi  co  rr  in  ‘x>  o-  co  cr*  o  — *  oj  oo  rr  in  10  r-  co  oi  o  — <  oj  co  *r  in 1 o  r-  co  cr>  o 
rr  rr  rr  -r  rr  t  in  in  in  in  in  in  in  in  in  in  1  o  'O  'Xi  id  io  'O  'o  >o  co  io  r-  o-  o-  o-  r^-  o-  r-  o-  o-  o-  co 


.>  .N  /«  _%•_%  . 


1810  I 

1820  END  IF 

1830  ! 

1840  !  For  LDV  data  after  the  first  data  set 

1850  ! 

I860  IF  Ldu$=" 1 “  THEN 


to 

at 

i 

cS 

*> 

n 

a 

V 

4* 

_l 

© 

> 

ff 

1 

© 

1- 

© 

n 

« 

© 

1 

© 

: 

t- 

<3 

♦» 

*> 

to  L. 

+  C 

lit :  to 

*»  to 

1©©  JQ 

to  £ 

*>  ©  ©  s_ 

lit  -rH 

lit  lO  liD  3 

1  L 

cS  1  1  C 

cS  to 

_ 1  =  II  1 

♦>  a. 

/*.  to 

Cf  X 

iJi-ec 

o  © 

to  ©  *>■  to 

il  ii 

at  cS  © 

*>  4* 

1  1-  TO 

to  to 

CS  ©  1  © 

at  •-< 

♦J-  ©  »-«  © 

1  -i 

CS  t—  V  © 

rS  <*. 

«  ©  o'  O 

♦»  I 

o  ©  u 

cs  a 

Li_ 

«x- 

►H 

i  il  to 
II  CS  GO 


lu  —  a 


a. 

CD 

II  H- 

2£  ID 

00.0.0. 
to  LU  LU  I- 
_£I  LU  LU  © 

oootoo- 


Q  CD  O  O  '©  O  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  ©  © 

r-  co  co  ©  -h  f.j  co  Tt  ud  •£»  r--  co  co  ©  — ■  oi  co  t  ud  id  r>-  co  co  ©  oi 

00  CO  00  CO  iX>  CO  iX>  OO  iX>  <X>  iX*  *Xi  'X>  ©©©©©©©©©©'-<’—•  »-« 

^  ^  ^  ^  ^  ^  ^  -H  O  J  OI  OI  W  W  IM  W  Oi  ©  CM  CM  CM  CVI 


2130  Check=Ra<I) 

2140  1=1+1 

2150  END  WHILE 

2160  ! 

2170  ASSIGN  PFile 


h- 

X 

x 

UJ  OJ 

• 

X  + 

1 — 1 

x  oc 

t—  *J- 

M 

Ul  UJ 

•H 

a*. 

j—  y- 

/\  i.O 

1—4 

03 

xx  o 

1 — 1 

XV 

x  x 

00  X 

•f* 

00 

xx 

-«  II 

X 

a— 4 

az  o 

cn  X 

£ 

\-X 

xo 

<3  *> 

3 

** 

o 

♦>  IJ4 

4-xX 

X 

i — 

<3  C 

n  ... 

X 

UJ  CO 

3  0» 

*— «  z 

o 

xo 

xv 

|XN_J 

-  X 

••w 

ox 

» — *  : 

M  » 

Z 

<3 

v  « X 

W  * 

X 

Q£U J 

*5- 

=  UJ 

: 

•w 

lilX 

<3 

a  ~x 

X 

Jfc 

CQ  1— 

3x-\ 

X  **  J- 

X  £  UJ 

s 

X 

1  4* 

X  X 

♦c  3  X 

3  O 

— •  X 

£  xx 

U»  X  •— ' 

u 

ZK 

v  £ 

xx  30 

C  ^  CO 

~7~ 

<3  3 

:  21  xX 

0>  X  3 

4—4 

< 

X  UJ 

x  x 

tH  V  V 

X  3 

CO 

X 

NX  V 

S  to 

X  o 

3 

} 

o  *-■ 

z 

II  o-* 

O  II  CO 

- 

Zl- 

— 1  UJ 

»a.  (it 

1—  o 

o 

UJ 

x_j 

3  II  O 

•4-1  | - 

CO 

>  --« 

co  <r 

3-  II 

3  X 

— 1  L  3 

«  n 

II  X 

J'HVlL 

II  Cu  X  •— < 

1— 

«  a. 

— ■  —  —  — . 

J*  ^ 

XX  10  ' 

1-4  £  1— 

3 

.  £ 

X  O' 

V  o  u. 

3  3  1— 

X 

•  3 

■■ 

£  C 

D.hQ 

O'  X  O  X 

1— 

3 

3  <U 

U.  X 

O  UJ 

3 

:  c f 

Of 

X_J  — 

i— i  UJ  — 

X  X  — 

O 

I  *>  X 

•  «  £ 

'<=»  3 


X 

£ 

X 

X 

1 

X 

3 

3 

3 

H 

C 

»- 

<=> 

UJ 

<v 

UJ 

X  - 

-  CO 

X - 

- UJ 

OOOOOOOOOOOQOiSOSJQOOOOOOQOOOO 

co  <r>  ®  oj  co  ’T  in  co  r-  co  cr->  ®  -<  oj  co  t  ho  *.o  r--  co  cr.  ®  oj  co  uo 

— <  — i  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  co  iV>  I>0  00  co  CO  CO  00  I'O  ro  ’T  TJ-  -f  rf  Tf 

OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ  OJ 


110 


f orma t ( 2x >  ' da t a ( ' # i 5 >  ' )  =  '/lfl5.5) 


AH  data  has  been  written  onto  disk 


i'S'wwv 

».\ 

)."  ' 


E 

n 

i 

i 


i 

I 

$ 


•.  • 
r-‘ 


tt 

I 

S$ 


m 

'.*. 

:.S 

1 

1' 

M 

O 


> 

ft) 

•w 

e 

•  • 

X 

<0 

_l  <0 

ft) 

44 

41 

Q 

*•> 

E 

X 

E  <0 

ID 

s 

4-/ 

X 

O  TJ 

C 

>» 

. 

u 

.  x 

> 

4^ 

<0 

•  <*-  C 

4) 

>- 

<0 

ft) 

X 

c 

<0 

4-S 

•  • 

> 

X 

ft)  C 

TJ 

CO 

Q 

o 

m  o 

C 

-J 

H- 

to  ■— 

o 

<0 

o 

ft)  X 

0rnm 

41 

-X 

t-  ft) 

CO 

X 

0>m 

• 

X 

(0 

*-> 

4) 

4) 

o 

c 

W 

<0  *** 

•— 

S3 

•  o 

TJ 

E 

c  — 

X 

to  u 

3 

CO 

o  44 

V 

o 

TJ 

C 

4) 

4-  <0 

44 

—  <0 

E 

44  44 

2 

c 

O  - 

O 

«* 

T3 

<0  c 

o 

c  44 

C. 

C 

4) 

*4  ft) 

C 

>  ft) 

>♦- 

4) 

E 

c  •- 

4) 

c 

ft)  c 

6 

c. 

ft)  x 

.* 

o 

a:  «) 

Jf 

•9m 

O 

—  o 

ID 

> 

CO 

X 

*- 

x 

w 

*4 

TJ  3 

•*- 

ft) 

o  — 

<0 

C  <0 

T3 

a 

4? 

to 

X 

E 

to  X 

X 

C. 

r—  44 

*  4-# 

X 

U 

4l 

ft) 

« 

«J  c 

O 

o 

to 

X 

u  o 

X 

«*- 

ft)  ft) 

ft) 

to 

M 

c 

•*»  44 

X 

4 

*4  *4“ 

4)  • 

-4- 

■  44  <0 

44 

«— > 

X  X 

c-  C 

•»•  L 

o 

ft)  o 

ID  O 

CT 

o  o 

O  Q 
—  x 
ft)  O 
>  O 

c 

C 


>•-  >  X 


o 

c 


c 

ft) 

to 

V 

x 


«J 

TJ 


I  I 
I  I 
to 

10  *-> 
*->  V 


ru  o 
■»  to 

%  TJ 

—  c 

«  <o 


UJ 

S 

> 

o 


«x 

oc 

o 

o 

K 

Q. 


nj 


UJ 

CO 

o 

a 

at 


3 

o 


<0 

<0 

a  x 

4> 

ft)  10 

H- 

ft) 

TJ 

X 

u  4) 

X 

to 

<0  - 

E  • 

>“ 

X  S3 

<0 

44  C 

4>  • 

Jt 

> 

>- 

E 

to  44 

ft)  o 

X  CO 

ft)  10  • 

Q 

>" 

«  3 

44  (0 

to  ■— 

4^ 

X  —  C 

— ) 

X 

X  C 

ID 

to  TJ 

44 

4> 

*4X0 

X 

X 

X 

TJ 

*0  *4* 

CO 

3 

X 

X  *4 

4-* 

C 

44  tO 

3 

ft)  C  *4 

<0 

T5 

ft  ft) 

T!  ft) 

(0  o 

O  ID 

*4  O  <0 

X 

X 

to 

4> 

TJ  > 

o  a 

A0 

3  TJ 

ft) 

E 

O  ft) 

4>  ID 

a  <o  — 

TJ 

X  <0 

3 

u 

^  T3 

E  *4  — 

(0 

ft)  44 

CO 

e 

o  *  ft> 

ft) 

X  <0 

CO 

•  4 

10 

L  V 

<_)  TJ  > 

GC 

3  TJ 

to  S3 

to 

s  o 

X 

ft) 

X 

E 

3 

C 


C 

4) 

E 

X 

ft) 

a 

X 

ft) 


<0  .O  O  TJ  < 


UUUUUUUUUUUUUUUUUUUUUU<J<JUOUU<JUUUUUUC 


113 


•  ^  .*•  .■■ 
>  *.»■>.»  Im 


/•  ’■  ,N 


data  set  number 
mont h  . 
day . 


yeoi'i 

x(mm)  . 
y  (mm)  . 
z (mm )  . 

Nozzle  temperature  (C) 


> 

O 

*■> 

rvi 

■#— 

• 

m 

< 

c 

•• 

«* 

<  -o 

4-> 

ft) 

C 

— 

■** 

m  <o 

a> 

O'  u 

u 

'  c 

c 

n 

4) 

TJ  — 

*— 

ID 

3 

t)  ft) 

£ 

C  •— 

L 

•#- 

3 

H-  >*- 

t— 

• 

o 

z  •  c 

e  re 

to  rv»  ry  « 

4»  *  «  e 

C.  *  «  w 

V  «  «  « 

Ifl  3  > 

3 

I  C  C 

I  use 

•  4)  4)  <0 

a  «  E  £  It 

yl  01  w  w  £ 

N  N  N  N  w 

£  fftHVl  V 

u  u  «  «  > 

w  w  «  «  > 

3  >  3  >  3 

C  C  C  C  C 

<0  <0  <Q  10  <0 

CD  CD  <b  4>  4> 

e  e  6  e  e 


<♦-  c  -c  —  —  e  c  o  a  o  c. 


T> 

4-* 

•*- 

4> 

4-> 

4-> 

4> 

s 

a 

<0 

e— 

X 

o 

4) 

X 

•• 

w 

O 

**» 

X 

VI 

4) 

c 

ft) 

^■1 

ru 

c 

*0 

(0 

c 

ft>  ft) 

ft> 

V 

e  e 

e 

to 

<0 

B 

w 

<0  (C 

ID 

10 

z 

3 

c  c 

C 

O 

♦-  >*- 

>♦- 

CO 

"O 

t_ 

■O 

JC 

</> 

3 

•  • 

• 

H- 

3 

CO 

<0  £3 

u 

JUUUUUUUUUUUUUUUUUUUUUUUUU 


-  -  „  *  «.  n  „ 

“  . 

-_L  - 

-T- 

TT 

r.% 

K-J 

jM 

* 

»  » 

L*'*1 

Ki 

• 

■ 

c. 

<0 

•V. 

-O 

$ 

> 

3 

c. 

10 

jO 

«• 

i 

=T 

ru 

> 

E 

c_ 

a  I 

>  w 

% 

o 

to  o 

£. 

H- 

r  » 

4-» 

10 

X 

* 

a  c 

.O 

4) 

c 

•  ] 

••• 

ru 

X 

<0 

•  t 

-  X) 

3 

4^ 

3  *. 

% 

c. 

(0  /-% 

y> 

4> 

wm 

•w  o 

E 

o 

JO 

a  o 

3 

E  ■ 

/*4 

—  IT 

(0 

3 

• 

.■w 

«h  w» 

% 

c_ 

C 

• 

n  w 

ru 

10 

c. 

• 

E  3 

E 

■W  • 

41 

r* 

4>  O 

3 

X 

C  t. 

13 

c 

►  % 

*-  Q 

Cfl 

<0 

41  41 

E 

4/ 

• 

•*-  X 

41 

E  X 

3 

n 

j|i 

4J 

w-4 

C_ 

—  E 

C 

E 

• 

%»  * 

E 

C-  3 

3 

r* 

a  /-k 

3 

4) 

41  C 

4-4 

C 

c_ 

X  o 

CO 

X> 

O 

c 

4J 

V  o 

%> 

X  4-4 

u 

+4 

£ 

—  rO 

> 

• 

O 

41  4) 

E 

4) 

E 

;; ' 

«*  w 

<u 

-y 

4-» 

CO 

•»* 

CO 

3 

X  to 

Tl 

co 

4) 

c_ 

C 

<0  -w 

4)  (A 

CO 

-C  <0 

4) 

<0 

F  10 

E  > 

X 

4> 

V  4-4 

O 

4-4 

4-4 

■ 

«n  X 

10  * 

V— 

<0 

X 

(0 

4) 

m 

4)  > 

c  > 

E 

in  X 

41 

X 

CO 

v  W  ^ 

-X  4J 

O 

*♦» 

*«■* 

,  •?* 

<0  * 

*■  X 

<- 

4J 

41 

O' 

<0 

4m*  CO 

( n 

*  X 

JZ 

c 

4-4 

,v 

</>  o 

4)  3 

O 

X  W 

4-4 

10 

E  KV 

6  % 

10 

% 

X 

c 

X 

4-*  W 

<0  C 

4-/ 

c- 

>* 

X  CO 

4-4 

c  c 

X 

•  4)  (0 

C  10 

(0 

•*» 

>- 

r 

3 

4-4 

ID 

X« 

13  W  4m> 

-X  4J 

X 

(0 

>- 

Q  Q 

a  e 

CO 

E 

3* 

*  to 

E 

o 

X 

4)  e 

C  x 

4)  *- 

<0 

** 

*■>  X 

—  > 

4) 

X 

c.  >- 

^  H-4  (L 

/-n  ao  a?  <-» 

_J 

<U 

4)  3 

% 

X 

4) 

X 

4)  V 

O  • 

-  co  rvj 

- 

CO 

« 

*  C 

*•4 

X 

X 

-C  x 

O  •- 

G 

^  o 

«k> 

•-  C 

t-  (0 

*- 

z 

3  ! 

— •  X  ^ 

X  4-S 

X 

oj 

o 

<u  a; 

X 

%r\j  ♦ 

rvj  *  +■ 

ru 

<« 

* 

c  •*«* 

*-  E 

<D 

sO  W  4 

>o 

W  *  X 

W 

\> 

4)  to 

O  3 

4> 

• 

W  to 

w 

4-4  Lp  W 

*- 

ir> 

V* 

0>  c 

<0 

a 

10 

4>  <0  w 

4) 

(0  W  4) 

<0 

w 

ii 

4j  4) 

c.  #— 

E  X 

4-» 

EX*' 

E 

X 

-a 

4-«  e 

<0  (0 

-  L  It 

t.  <0  — 

C- 

to 

3*1 

c  — 

X  <0 

• 

(-0  4) 

c 

O  4)  <- 

o 

a 

—  T > 

u  c- 

3  >*-  C 

I  --  C-  3  <*- 

C 

’  •* 

•:■} 

O 

ru 

o 

o 

o 

u  u  u  u  u  u 

u  u 

u  u 

u  u 

c 

c 

c 

c 

to 


IS 


115 


Form  the  vertical  file  namet 


cont i nue 


-  •  *V  V 


■  a 
?• 

a> 

JC 

a  v 

0 

E 

w 

V 

i 

*a 

4J 

w 

c 

® 

o 

<D 

<i» 

c 

<b  T> 

a* 

C 

a* 

• 

4-> 

W 

y  ; 
• 

a 

•o 

to  C 

c 

o 

<D 

<D 

I*  't 

10 

® 

—  <0 

II 

3 

o 

X> 

k'-j 

w 

V 

aa 

K> 

V 

•-* 

ru 

> 

<0 

•  ‘ 
V 

<0 

b 

—  0> 

*-• 

•— 

1 

*U 

f\J 

10  "D 

E 

<r~ 

>  — 

4 

e  a 

V 

«■  <0 

* 

a- 

r<\  * 

• 

« 

ro 

P 

1 

0)  "a 

3 

C 

>0 

a 

* 

0> 

* 

II  ♦ 

a 

II 

u  <0 

o 

4a  <a-  a 

w 

a 

■k 

o 

c 

«x 

“a 

o 

u» 

-a 

z 

4>  a 

<U 

II 

• 

• 

,*  •— 

o 

:*. 

4a 

c  *a 

10  —  X 

*J 

II  "O 

o 

■O 

U 

4a 

)'•, 

c 

II 

<0 

•-  —  rvj 

•*- 

.*  <T> 

V 

V 

-* 

<0 

0)  II 

[*- 

o 

3 

II  4a 

• 

®  w 

<- 

c 

0) 

3 

II 

3 

u 

0) 

-C 

c 

oi  «*- 

w 

(A 

to 

T3 

II  o  4a 

2 

3 

^  L 

c 

C 

0) 

(. 

u  •*- 

CT  k- 

r- 

•*- 

4> 

E 

3 

<0 

N— ' 

O 

-*- 

-K 

•■» 

aC 

•a 

<0 

4) 

to  — 

4->  E 

c 

O' 

4a 

u 

4^ 

u 

■O 

o 

** 

u  — 

V  b 

*1 

C 

«a 

0) 

c 

w 

TJ 

*3 

• , 

c 

•— 

a  <o 

to  o 

o 

o 

o 

II 

JC 

o 

H- 

C 

•* 

►  J 

V 

o 

c 

—  o 

•<-  a- 

0 

"O 

o 

u 

u 

t) 

u  u  o  u  u 


v  *-•  ;v  y-»  k  »  v*  *  *  k  •  »  •  -  -  *.  v 


,s  *  •  ’  *■  * 


c  — 

V  •>> 

-C  O'  —  10 


O'  I  «3 
«  O'  "O 

x  o'  >  > 
io  i  «n  *-* 
E  *-  in 
>11  Q 

V  J-> 

"D  <-*  <0 

II  v 


0 

* 

0) 

t 

w 

<0 

D 

*-* 

<0 

E 

> 

cn 

« 

> 

O 

*-» 

*— 

rvj 

V 

<0 

• 

<C 

ST 

a 

<0 

Q 

*0 

T3 

"C 

a 

<T5 

—* 

u 

>*• 

V 

> 

(/> 

•— 

> 

> 

U  H- 

o 

a: 

> 

a> 

o 

4»» 

^4 

• 

■c 

4-o 

T) 

a> 

-D 

II 

a 

V-/ 

"O 

C 

o 

t 

D 

+ 

<0 

c 

o 

*> 

O’  *- 

.o 

C 

X 

• 

4) 

Oi 

o 

«*• 

II 

<0 

w 

E 

T> 

II 

C 

II 

> 

> 

*0 

c 

o 

0) 

4> 

C 

tl 

o 

T) 

*0 

4) 

u  u  u  u  u  — 


Compute  the  turbulence  intensities 


<D  Ct)  <D 
■O  T5  TD 
D  >  > 


i 

1 

1 

c 

c 

c 

0) 

<0 

•D 

ID 

E 

4 > 

41 

4) 

4) 

E 

E 

E 

" 

c 

4* 

3 

> 

> 

o 

•4 

w 

%■* 

w 

»*• 

*— 

« 

« 

■It 

4*  • 

/•* 

CD  s 

II 

»-s 

c  /-s 

- 

— 

«*■ 

*#• 

L  O 

<0 

w 

w 

w 

o  o 

• 

E 

<0 

ID 

ID 

<4  in 

4) 

w 

fci 

* 

C  w 

4 

ID 

ID 

ID 

*•-  4 

■c 

-C 

T5 

D 

C 

3 

> 

3 

4*  O 
C  C 

C 

1 

1 

1 

/-N 

O  x 

4) 

/-N  C 

C  -N 

c 

E  O 

XT 

C 

c 

C 

co  cd  co 

CD  CO 

<D 

*4  S 

4 

ID 

ID 

ID 

EOF 

O  E 

<D 

c 

4> 

41 

4) 

o  E  O 

E  0> 

E 

O  4 

B 

E 

E 

4^4 

>  4* 

> 

G  >■ 

co 

3 

> 

3 

•#»  •*  — 

«  -e- 

« 

x  CD 

E 

w 

V-# 

C  C  C 

c  c 

C 

o  C- 

O 

w  <D  W 

CD  w 

CD 

u 

4 

■*■ 

4#  O  4 

O  4* 

O 

O  CD 

co 

CD  E  CD 

E  O 

E 

j: 

c 

E  — 

f\J 

rO 

O  D  O 

>  O 

D 

4  C 

41  E 

E 

E 

»-  W  ^ 

W  r— 

s«* 

CD 

• 

~  3 

3 

3 

"4  v.  4 

\  4 

N 

O 

—  10 

(0 

(0 

v  C.  \ 

c  \ 

C 

c.  O 

*_ 

C 

(D  f\J 

<D  i*n 

CD 

O  4 

• 

»  II 

II 

II 

ed  e 

X)  E 

X) 

4  C 

o 

o 

o 

D  r\J  D 

f\J  D 

> 

CO 

• 

• 

• 

V) 

II  — 

OJ 

m 

co  d  co 

>  CO 

D 

o 

o 

o 

E 

—  e 

E 

F  4; 

4> 

3 

3 

3  3 

ii  ii  ii 

II  II 

II 

• 

II 

II 

II 

4 

O  10 

(0 

(0  c 

u 

o 

•*“ 

L  L  L 

c.  c 

C 

** 

rvj  ki 

p» 

o* 

4-» 

IQ  <0  <0 

CD  CD 

CD 

6 

e 

e 

W 

c 

£  3X 

-Q 

D 

3 

3 

4 

o 

o 

rvj  rvj  r\i  rvj  > 

> 

in 

CO 

CO 

TS 

u 

3  3  > 

>  3 

D 

O 

o 

u  u  u  u  u  o  uu  cruu  uuuuuou 


120 


* 

4 

I 


ifdset  .eq.  setmintl)  then 
iout  =  0 

el  se 

iout  =  iout  -  1 


Write  out  the  means  and  standard 
deviations  onto  disk. 


o  o 


■ 

4-4 

o  « 

4  * 

•> 

nj 

<M 

E 

4-* 

O  Q 

ID 

*  4/ 

X 

C 

in 

«  4/ 

-X 

•  #- 

% 

4-» 

a 

o 

4>  1 

X  1 

a5 

X 

cn 

4/ 

Ml 

4> 

—  o 

•4-  O 

• 

•*» 

o 

o 

% 

i  - 

1  — 

4-4 

a 

^■s 

« 

4—* 

4*4 

4*» 

fO 

X 

O 

»4  ® 

CC 

30 

O  4*S  « 

X 

X  X 

cn 

-X 

V 

o 

o  4-» 

o  o  a  =r 

3* 

•n> 

* 

-X 

4> 

—  4) 

«•  X 

W 

f\J 

E 

\ 

\  in  + 

4- 

ro 

S  \  4) 

<4-  «♦» 

t> 

-X 

4-» 

<D 

4-^ 

4-4  •*■ 

-X 

a  a  — 

i-i 

E 

* 

4» 

C 

1) 

dj 

Oj  l^l 

\ 

X  X 

r\j  ro 

<D 

•M 

in 

.X 

in 

in  i  4~ 

4-4 

4*4 

\ 

4)  V  1 

0 

o  o 

91 

C 

-X 

-X 

4) 

4- 

4)  (M 

•4- 

X 

X  X 

w* 

^  4-4  in 

in 

in 

-X 

w  v-4  Q 

4) 

4<  4) 

1*  *  •; 

4-» 

*->  <u  •<- 

•^* 

\ 

W  4-4  X 

*P* 

V 

4t 

« 

* 

c 

Cl/)'-' 

V- 4 

w 

\ 

C  C  4) 

w 

N— 4  W 

!%" 

c 

L. 

c 

L. 

•4* 

•<-  — ■  C 

c_ 

C 

*— * 

C- 

L  C- 

!*,  * 

•4- 

0) 

<y 

4) 

ID 

ID 

ID 

X 

ID 

to  to 

'*  4 

4-4 

4-4 

tl 

II  II  X 

X 

JZ 

II  II  II 

X 

JZ  JZ 

'*.4 

3 

u 

O 

U 

o 

o 

U 

II 

u 

u  u 

o 

<0 

<0 

ID 

nj  ro 

—  f\l  1-0 

Ml  a 

c 

c 

l_ 

c 

4-4 

*-’<->  it 

II 

II 

4-4 

a  a  a 

II 

II  II 

fe?  ox 

■O 

ID 

<0 

ID 

dj 

4)  4) 

4/ 

XXX 

iv  *-  c 

3 

-C 

X 

in 

in  in  «-*  f\j 

rO 

in 

4)  4)  4) 

r\j  r4| 

;  -j  #>  4/ 

</) 

u 

u 

u 

•*•  jt 

-X 

-X 

-X 

-X 

-X  -X 

M 

[v 

B 

1  -  O  U  U  U  U  O  O  O  o  uu 


123 


I 


kexp  =  kl//k2//k3 


d'//kexp//kset//' .dot 


program  entrain 


4> 

b  4>  x 

I)  £  w 

o  > 


m^ 

a 

c.  *— 

44 

> 

4) 

E  N 

o 

a> 

O  C 

<A 

*■» 

> 

•00  *• 

C 

> 

o  o 

■D  ffl 

b  0-0 

% 

c. 

44 

(A 

4>  E 

4> 

a  o 

(A 

a 

3 

(A 

b  b 

<0 

C 

3  O 

3 

o  -o 

C  • 

O  O 

AJ 

> 

O 

a> 

O  <n 

J£ 

~  c 

-O 

•— 

^  w 

c. 

c 

c 

•*“  w 

tn 

tn 

c 

c. 

O  4) 

4> 

X 

m  a> 

• 

*> 

<0 

4) 

c  u 

c. 

44 

to 

tn  •— * 

TJ 

*o 

tn  b 

f\i  <o 

a 

<A 

44 

4) 

4) 

to 

c 

c 

w  b 

X 

a 

to 

<A 

JD 

b  V 

b 

E 

4) 

4> 

>  *•> 

t) 

44 

3 

O 

a  x 

4) 

O 

c  tn 

w 

U 

0*  >• 

E 

to 

<3 

X  *l 

■o 

c. 

o  c 

4) 

% 

X 

O 

<A 

a> 

"4- 

•*»  o 

E 

4> 

O 

c 

3 

e 

«l 

* 

i->  -0- 

<G 

0mm 

o  o 

> 

E 

> 

(0 

X 

, — 

c 

*0 

to  tn 

M 

ru  o 

3 

X 

44 

a 

o 

<0 

—  c 

■O 

M 

w  ru 

<0 

tmm 

44 

c  w 

u 

a> 

3  4> 

o 

3 

T3  X 

(A 

(A 

a>  c 

c  •• 

X  E 

4) 

c 

*■  4> 

% 

% 

tn 

C. 

44 

c  o 

« 

E 

to  — 

«  Cl 

O 

<0 

a> 

c 

4)  -•» 

•00 

'0)  t 

•w  T3 

M 

4> 

e  <c 

E 

E 

u 

a> 

CT  ~ 

X  o 

M 

X 

O  b. 

44 

3 

3 

4! 

s 

w  < — 

o 

44 

ryj  *-< 

•— 

* 

Ti 

o 

u 

C  X 

(0 

o 

X 

c.  — 

c 

W  X 

X 

<0. 

c 

0) 

V 

•«— 

« 

44  4) 

X 

o  «* 

o 

M  % 

O’ 

44 

E 

c. 

to  c 

u 

X 

X 

a 

4> 

44 

*i  0-0 

4; 

c 

44 

to 

*4- 

w  w 

a>  4* 

c 

4) 

X 

0-0  O 

b 

3 

«k 

E 

4-* 

4> 

x  tn 

o 

#— 

<0 

b  • 

4-* 

44 

o  o 

«. 

O 

to 

fO 

<0 

o 

o  — 

•00 

4> 

s 

U 

o  ru 

3 

o 

44 

c 

*c 

u 

■o 

b 

*-*  o 

E 

b  w 

o 

4> 

ru  w 

o> 

% 

E 

a 

o 

*0 

> 

O  <0 

w 

a 

w  tn 

4) 

4) 

c 

3 

o 

> 

a> 

*0  > 

4) 

<0  (A 

X 

(A 

>  b 

b 

E 

4> 

w 

C. 

44 

X 

*0 

E 

w  41 

• 

4^ 

4> 

k  w 

o 

<0 

4 

a 

w 

tn  ■— 

E 

<0  — 

~  c 

u 

c. 

U» 

ru 

c 

C 

c 

u 

4>  U 

V 

*c  •- 

3  b 

4> 

O  >*  X 

W 

-X 

O 

44 

(A 

c 

c 

>  O 

W 

«  .  . 

b*  ^ 

w» 

c  4; 

4_ 

a> 

b  «» 

*00 

4) 

~  c 

<A 

44 

ru  c 

-X 

3 

X 

a> 

X 

3  4> 

X 

X  C 

3  4) 

4> 

w  % 

X 

—• 

<0 

(A 

> 

•w 

O  > 

<0 

V- 

O  = 

(. 

X  4-V 

* 

44 

% 

o  ru  o 
c  o  to  —• 

O  <"U  *  w 

w  a>  (o 
C  8 

4)  —  —  — 
E  L  IC  # 
—  0  4>  41 
"C  >  b  b 
O 


<  <c  <c 
b  TJ  *> 

0)  *. 

«-  b  C 
u  a>  tj  c 
to  a  o 
c  <u  —  e 
<e  **  «o  •£ 
x  c  a>  o 
o  —  b  o 


uuuuuuuuuuuuuuuuuuuuuu 


i di aqn 


>„’.N 

.**  s 
►  .  ■; 

a 


v, 


4) 

c 


« 

4) 

3 

ID 

> 

4) 

c 

3 

** 

ID 

C 

41 

a 

E 

4) 


4) 

X 


>-<  w  «n 

v  v  L. 

E  E  0 

'*■  •*"  0 

w  L  (.  >. 
w  w  Q  Q  D 

3  >  3  >  C 


«  «  «  «  « 


■c  X3  T  -O  -O 

ID  <0  ID  ID  ID 

4)  4)  4>  4) 


« 

o 

4-> 

CP 

c 

•0 

c 

c 

3 

C- 

C  ID 
ID  **  *■> 

■*•  -*- 
4)  C  <0 
*-«*-*-' 

3  E  e 
Q  3  3 
E  W  0 
O  © 

U  II  II  *— 

CO  o 


P  E 
3  3 
0  0 


O  — 

a 

■o 

c 

4) 


C 

3 

41 

0 

O 


c 

4) 

E 

(. 

4) 

a 

X 

4) 

4) 

X 


CA 

c 

4mf 

o 

•f* 

C 

3 

0 

. 

4> 

»“4 

c 

CO 

3 

c 

ID 

•0* 

C 

V 

0 

o 

a 

4b* 

o 

E 

c 

t 

o 

4) 

4) 

o 

4m* 

•  E 

o 

O 

04 

4) 

o 

O 

V 

C 

L 

00 

04 

04 

ID 

3 

'v 

•0 

4) 

«A 

04 

•0 

E 

ID 

•P" 

*0 

w* 

4) 

•0 

•  0 

(A 

4) 

E 

© 

© 

© 

N-4 

C 

-C 

o 

o 

© 

© 

O 

4> 

0) 

*-» 

•— 

o 

o 

© 

o 

o 

E 

E 

(A 

•— 

04 

0-0 

•-* 

00 

0* 

•0 

•0* 

> 

4) 

» 

'v 

N. 

N. 

c 

c_ 

4) 

V 

/■N 

00* 

a 

Q 

(- 

3 

4) 

•0 

•0* 

•0- 

•0- 

•0- 

3 

> 

II 

a 

o 

•0 

•0 

•0 

•0 

•0 

•0 

n 

•I 

04 

E 

•0  »0 

<0 

% 

S-0 

y-4 

N-4 

w 

w 

0"V 

*S 

•0 

O 

\  V 

0— 

—4 

X 

> 

M 

3 

> 

•0 

•  0> 

•  0 

u 

ID  C 

CL 

•00 

•0 

4m*  +0 

II 

II 

it 

it 

II 

ii 

w 

if) 

E  E 

4) 

4) 

c- 

3  3 

•0 

00* 

F 

F 

V 

cn  (A 

•0 

•0- 

•0- 

•0- 

•0 

•0 

(A 

.0- 

•0 

»#0 

•0 

c 

l. 

> 

•  1  II 

rvj 

w 

p— ' 

w 

w 

V-4 

a 

a 

0; 

04 

X 

> 

M 

3 

> 

3 

> 

L. 

o  c 

E  E 


O 

■o 


o  u 


u  u  u  u 


u  u  u 


I 

!s 


129 


f*  •  •  "  «  “  •**  •*  ^  *“  p.-*  ■»  4  V*,  ^  *  •  •  v  *  *  "  *  « 


•  l  •  - 

»JUi  1. 


6/S*(?£-uu,<)  = 
fe/S*(2£-*u,»)  = 
0001/P 


[ 


c 

•—  II  <0  c 

4-  E  6 
C  "C  **  4- 
O 
o 


4 

«-4 

© 

t r 

r*- 

>»■ 

o 

O 

rP 

X 

rg 

• 

* 

4) 

• 

I 

•* 

* 

-f* 

rp 

<0 

-0 

4-> 

rP 

•* 

£ 

a> 

4> 

* 

•K 

■LS 

© 

— , 

* 

c 

* 

• 

C 

E 

P- 

44 

e 

t 

V 

c 

« 

4> 

r-N 

« 

P* 

rg 

rP 

c. 

in 

1 

«■* 

* 

«f 

« 

a> 

© 

« 

X 

0) 

i\j 

rP 

£ 

E 

© 

• 

4- 

3 

CC 

© 

4* 

< 

Z 

ru 

rP 

Pi 

fP 

© 

• 

« 

© 

4-S 

« 

• 

4* 

«* 

© 

4-N 

-o 

1 

rg 

<0 

© 

X 

. _ 

rg 

E 

© 

e 

o 

« 

4-J 

• 

<0 

c 

* 

c 

« 

4* 

o 

> 

c 

E 

IP 

PJ 

c. 

4) 

p 

1 

«C 

£K 

4-* 

Q> 

<* 

4-S 

« 

IP 

© 

6 

o 

T> 

Pi 

• 

IP 

H- 

L. 

C 

© 

4> 

O' 

1 

* 9 

GT 

© 

© 

© 

X 

© 

IP 

• 

p- 

4 

E 

o> 

p* 

© 

1 

Pi 

10  K> 

4> 

o 

© 

(0 

© 

o 

• 

•4" 

o 

• 

E 

© 

c- 

© 

44 

o 

1 

<+•' 

• 

w 

* 

c 

9 

C 

* 

© 

l 

w 

« 

flj 

«♦* 

E 

GT 

© 

e  ■v 

C 

« 

c. 

Q 

C 

1 

© 

*♦- 

rvj 

a; 

*o 

H- 

O 

£ 

0) 

«■* 

4C 

X 

« 

V 

c 

cr 

O' 

4T 

P* 

4C 

•-< 

rP 

a 

« 

« 

w- 

/"N  tP 

■c 

© 

© 

IP 

4> 

IP 

©  IP 

* 

<o  p* 

O 

• 

© 

a> 

sr 

cr 

fp 

GC 

O' 

•*-  © 

IP 

O' 

• 

44 

00 

K> 

•«» 

• 

« 

O' 

in 

>oc 

» 

w 

© 

3 

a* 

1 

iT 

4- 

<0 

• 

\  IP 

• 

O 

4* 

a 

o 

4)  ro 

© 

© 

E 

4- 

© 

o  ru  v 

\ 

t. 

r- 

E 

• 

o 

• 

© 

O 

4-< 

1P 

L. 

c 

c 

a  it 

O 

1 

• 

-f- 

© 

o 

« 

0* 

• 

«  4 

6 

E 

0) 

o 

(_> 

IP 

© 

*«* 

© 

© 

p* 

rP 

"D  4> 

3 

3 

V 

• 

«■« 

•« 

o 

« 

• 

w 

«  a 

X 

o 

«-* 

1 

© 

\ 

C  • 

II 

II 

<T5 

It 

• 

CT 

o 

5T 

w» 

4> 

• 

a 

E 

E 

« 

• 

* 

II 

i 

3  to 

X 

X  4- 

Q 

•* 

«■* 

-Q 

C* 

n 

C. 

<0 

eo 

<0 

II 

C 

c 

E  r- 

II  >— 

E 

<D 

F 

II 

X 

CO 

E 

II 

E 

«o  a 

H 

c 

3 

to 

3  T3 

a 

*-> 

o 

o 

• 

E 

to  <*- 

C 

c. 

0— 

> 

«* 

l 

•* 

c. 

1 

a 

3 

L.  .m. 

a> 

a*  ■*- 

<0 

u 

o 

u 

u  u  u  g 


o  o  o  o  o  l 


Initialize  file  containing  the  error  information 


y  o  u  o 


rs 

<0 

m 

* 

— 

4> 

0> 

2 

JC  • 

c 

c 

X** 

*■<  4> 

• 

1— > 

ST 

*— 

II 

K 

4-  •— 

<n 

<r 

•> 

O  —  *4- 

3 

m 

•fc'  u.  o 

4^ 

a: 

C  c. 

<0 

o 

c 

--  a. 

4~* 

LL. 

0) 

• 

z 

-O 

>  >•  c 

* 

•— 1 

e 

4->  (0  <0 

*» 

■  L  •— 

c 

H- 

z 

o  t.  m 

o 

l—l 

o  <o  in 

c 

u. 

4^ 

—  3 

c 

c 

4>  —  <0 

a> 

2 

a> 

>  10  CD 

• 

-I 

6 

c 

n 

i—i 

Q  — 

»•  O  fl 

« 

to 

x  c 

o  — 

•— 

co 

4>  0) 

in  o 

*— 

13 

Cl 

w  c  *•> 

h- 

<t 

X 

41  OJ 

*-*  CD 

^  uJ 

3  f  a 

ST 

3  - 

•-I  • 

—  c 

II 

O  - 

O  * 

It  TJ  — 

*< 

rvi  x 

r\j  x 

>  -O 

-  nj 

*  rvj 

41  c. 

c 

sr  — ' 

cr  w 

T  C  O 

3 

w 

<0  0  0 

w 

u 

4>  u 

4)  <0 

a>  <o 

a:  io  io 

C 

*->  F 

4-  e 

« 

i-  fc. 

•4-  L 

O 

t-  O 

C  O 

• 

o 

*  ■*- 

X  ^ 

K* 

o 

, 

o 

o 

<\l 

rvj  u 

C 

C 

c 

c 

c 

*j  x  3  O’ 

C  4> 

3  II  II  •  —  — . 

—  O 

— 1  + 

II  4)  -X  -X  c 

II  3II4-4-WW3.*  — 

4)  C  4)  4>  O 

O©*-1®-—  U  U  U  II  II 

<0  C  C  <0  <0  w 

L  II  3  II  C  3  II  II  -*  — 

*j  O  O  O  ****'•- 

.x  u  —  u  u  -x  •—  x  > 


x  t  race ( k  )  =  x  ( 1  ) 


^ 

f 


U 

4) 

/■N 

**• 

41 

41 

4^ 

• 

O 

u 

4) 

■O 

<0 

4) 

c 

c 

^  (- 

C 

<0 

4)  ^ 

45 

^  > 

O  > 

> 

45  * 

4)  * 

O 

O  4) 

C  4> 

4) 

•#- 

<0  O 

4-»  U 

C 

u 

C  <0 

•f  4) 

4) 

o 

C 

w  L 

X 

*— 

rvj  4> 

•  > 

Si 

w-  X 

X  X 

> 

<\l  * 

*  % 

4^ 

-X  ^ 

<*s 

u 

V 

%  p— 

41  — 

•— » 

c 

/-\  W 

a  w 

-r- 

4l  M 

45  N 

n 

«— . 

O  *• 

C  % 

c 

c. 

4) 

V  *-s 

<0 

«> 

C  45 

—  4> 

4-  U 

w  o 

10 

c 

—  4? 

^  to 

V 

a> 

w  C 

C 

•X  c 

**• 

o 

*-•  4^ 

45 

*  4* 

** 

X  — 

x 

4)  — 

*#* 

o 

^  W 

O  w 

o 

4)  OJ 

4i  ru 

o 

U  -* 

C  X 

<0 

4?  - 

4-»  % 

•» 

V 

C  ^ 

p— 

>  ^ 

> 

3 

4^  a> 

w 

*•  45 

• 

r— 

X  O 

N 

4)  O 

« 

(4 

<0 

/-n  *  45 

^  U  4) 

c 

45 

> 

—  «>  «- 
w  vj >  ** 
N  <0  *«* 

*.  w 


<0  C 

e  — 


/•N 

II 

4k/ 

II  *- 

4- 

XJ 

* 

45 

*_ 

X 

-X 

*  -X 

4) 

C  C 

c 

O 

^  , 

s-* 

/•N 

% 

% 

/"N 

*-< 

^  X  % 

U 

45  45 

45 

o 

45 

■•s 

D 

45 

-X 

X 

c 

45  w  _x 

45 

O  X 

••• 

C 

%«/ 

u 

w 

v— ' 

II 

3 

U  W 

C 

o 

*-» 

4^ 

^4  . 

II 

<0 

4-1 

V/ 

O 

4)  4-» 

4-/ 

45  a 

X 

(ft 

X 

c 

•*• 

c 

o 

L  ^  L 

X 

% 

%*/ 

r-s 

4-# 

H- 

O 

c 

4-/  h-  O 

4-/  O 

E 

c 

a 

X 

•«- 

c 

D 

M  IP 

—  DC 

II 

IP 

+* 

c 

45 

X 

v- ' 

V- / 

41 

c 

w  «-t 

w  45  C 

*4 

4/ 

O 

4-» 

II 

45 

45 

4) 

CT 

4) 

41  CT.  4) 

4) 

eft 

o 

II 

o 

<J 

4; 

H-  O 

U 

o  o 

0 

X 

tn 

<p 

— » 

/“N 

4) 

4) 

*— 

^— ■ 

</> 

•—  *J 

45 

45 

w  ■*- 

w 

c 

c 

C 

C 

p— 

#— 

O 

L 

C 

II 

at  -o 

to 

45 

o 

W 

4-» 

4-# 

4) 

4) 

p— 

it  O 

4) 

4k/  45  45 

*/ 

o 

O  — 

c 

X 

rvi 

^4  . 

> 

M 

u 

o 

o 

V>  O' 

(ft 

N  O  O 

— 

X 

O  "C 

2  J 

H- 

o 

X 

uoooooooooo  o 


k2 ( j ) =  sart ( l /k2  C  ]  )  ) 


cont i nue 


<n  <*- 
v  oi 
—  o 
n  o 

I0t 

Cl 

—  c 

in  4i 
■C  >  6 
*  C 

M  — 
«•  tft  <0 
Cl  Cl  t 

C  C 

«*_  r»  nt 


uruoooociocio 


-X 

(A  41 

V 

CO 

C  -C 

ru 

4>  4-» 

-O 

E 

% 

D 

X) 

H- 

m 

C 

O  "C  O 

s 

4> 

4) 

4-* 

• 

• 

c 

c 

m  a> 

2 

% 

K 

% 

<0 

<0  3 

4) 

*“N 

a-  ru 

ru 

ru 

4m* 

«— 

C 

<0 

-O 

4) 

in  id 

• 

I— 

%  % 

* 

-D 

c>  > 

II 

3 

•  » 

10 

0. 

a  o 

<0 

JO 

c. 

*■»  c 

D 

I— 

UJ  UJ 

% 

% 

c 

•«“  ID 

4m* 

3 

CD  *- 

4-4 

ru 

#— 

k-'  Cl 

<0 

c 

S  «_> 

J* 

.x 

D 

C  6 

4m* 

3  3 

% 

% 

x> 

(D 

CO 

z 

Z  O 

4> 

4) 

to 

3  Cl 

4k 

t— 1 

Z 

u 

U 

+* 

a  jc 

• 

<. 

1-  o 

<0 

<0 

M 

c 

QC 

z  u 

c 

c. 

41 

01 

c_ 

H- 

UJ 

V 

4-* 

£ 

£  01 

4m* 

Z 

S  «-  1- 

M 

N 

4m* 

«•>  4-> 

c 

Ul 

w  >Z 

% 

% 

3 

4) 

a:  -  lu 

4> 

4) 

4- 

O'  Q 

• 

S 

Ul  <0  S 

U 

U 

3 

C  6 

II 

«r 

aarw 

<0 

<0 

O 

—  O 

4> 

QC 

X 

X  *  £Z 

c. 

c_ 

*■>  O 

*— 

o 

01  UJ  O  aJ 

4m* 

A# 

41 

3 

•*- 

o 

-  e  a. 

4> 

Q  O 

cr 

k  X 

w 

w 

6  in 

4> 

a. 

#— k 

X-'UI 

c 

o  — 

ru 

O 

» 

4— 

rvj  ru  - 

4-* 

4m* 

s 

U  <1 

II 

O— 

4-4 

*  «-»  k 

4«» 

X 

—4 

x  —  x 

* 

ru 

*  N  *  CU 

* 

2 

c  ru 

w 

ru 

w  ru  w 

o 

O 

3 

w 

w 

w 

Q 

a 

tv* 

4m* 

4n*  4<k> 

41 

<0 

41 

<0  4>  <0 

— 

C 

4-/ 

E 

4m* 

E  4-*  E 

•— 

01 

c 

»-• 

C.  •*-  C. 

<0 

<0 

a 

C- 

O 

c 

O  t-  o 

o 

u 

o 

2 

2 

*•»  2  *•“ 

o 

~  AJ 

^4 

0—4  *4 

u 

o  u 

c 

c 

c 

c 

c 

c 

u  u 

o 

•4  ^ 

i 


Stop 


subroutine  f  name2(  i  exp/ knarne) 


f  v'-*.TS?r«L’ri 

if 


m  . 


[/■'  , 


ft! 

•*» 

E 

Q 

« 

X  1 

C 

4) 

-V 

—  o 

% 

o 

O 

1  »« 

X 

«  a  *■» 

/■s 

4> 

K> 

o  ^  >  ec 

3C 

CO 

a 

,X 

o  o  a  5T 

O’ 

X 

v-* 

—  «  X 

a 

w 

r\j  e 

\  V  4>  4- 

4* 

-X 

4» 

-X  G.  TO 

a  o  — 

-X 

s 

e 

^  x  C 

XX  •— 

f\J  f«"l 

\ 

s 

<0 

^  41  x 

4)  4)  1  0 

a 

o 

V 

- 

C 

-X  -X 

•*■  ■>■  X 

X 

X 

rvj 

H- 

—• 

w  w  0  4J 

4> 

4> 

-X 

0) 

u  u  X  ■•■ 

s 

o 

4) 

*  *  <* 

C  C  «i 

>w 

W* 

s 

u 

C 

C  (.  c. 

•*•  •#•  c 

L 

L. 

«■* 

• 

41  4l  4l 

<0 

TO 

TO 

-X 

u 

v  w 

II  II  II  X 

X 

X 

II 

3 

o  o  o 

u 

o 

a 

II 

O 

10  <0  TO 

(VI  ro 

4) 

c 

L  L  C 

a  q  a  ii 

II 

ii 

a 

E 

J3 

<0  <0  <c 

XXX 

X 

<0 

3 

£££ 

4i  4)  4>  ^ 

r\j 

hO 

4) 

c 

» 

o  o  o 

**•  -X 

-X 

JX 

-X 

je 

\:C 

$  1 

ir.‘ 


o  o  u  o 


oooooooot 


do  540  i=l,nitems 

if(x(i)  .It.  xmin)  xi»1n  *  x(U 
340  continue 


c 

E 

» 

II  C 
c 

C  3 
•.  w  TJ 
eve 
•o  c  v 


c 


o 

-o 

% 

4>> 

c 

E 

« 

o 

o 

c 

c 

* 

V 

•*>  •• 

T3 

c 

•O  6 

% 

% 

*-  c. 

« 

IT 

fVi 

o  o 

E 

< 

o  *■ 

«l 

3 

u 

r 

<0  « 

c 

« 

< 

x 

V 

a- 

% 

*  «-> 

E  -* 

X 

o 

3  » 

% 

X  X. 

*>  E 

N 

w  O 

*  V 

% 

N 

<\J 

K>  ~ 

C 

>  « 

« 

F  — 

E 

> 

•* 

o 

3  C 

O 

<0  c 

u 

(O  w» 

4* 

**  3 

X 

^  M 

10  u 

« 

V 

rvj  «. 

c 

T> 

X 

E  « 

W 

u 

<1 

w 

3  « 

4-* 

V  v 

«rt  E 

X  X 

II 

c 

-  V 

*4» 

«-  o 

*.  • 

•«  *J 

a 

3  * 

W  E  •<« 

o 

V 

X 

**  rvj 

E  3  C 

Q 

•-  C 

V-/ 

V  «f 

0) 

U»  V 

N 

ct  * 

w  »  x 

€ 

—  OJ 

c 

c  «*  c 

•#* 

•• 

%  o 

V  <  Iff 

a  c 

V  —  V 
<0  E 

C  V  •- 

H»'  c  ■o 


uuuuuuu  uuuuuoooou  ouc 


138 


10 

£ 

V 

V 

•*■* 

/-s 

C 

/— <» 

e  v 

o  aj 

w 

c  s  « 

X 

DO* 

w 

T3  CiO 

*->  C*  ^ 

V  D  E 

ru  o  rvj 

— '  "C  3 

♦  —  * 

*-  —  V  10 

-It  <0  * 

E  E  ~  ^ 

3  3—1 

*  »  F? 

«  «  3  E 

W  •*•  •#-  w 

-*  9  !ft  3 

WWW  X 

E  E  «  10 

W 

w  X  NJ  w 

3  3  m  « 

0s  w  w  g 

«  ID  EM 

CT 

o  a  a  o 

1  1  3  <0 

O 

—  o  o  — 

rvj  nj  m  ■*■ 

#—> 

40  —  « 

e  e  i  i*o 

W  <0 

w  10  ID  W 

3  3  9  E 

e 

10  10  E  3 

i>  + 

♦  +  ♦  ♦ 

*  *  3  10 

** 

10  K1  10  * 

o  o  o  o  o 
•  •  •  •  • 
o  o  o  o  o 

II  II  II  II  II 

»-  rvi  ^  3  i/i 
e  e  e  e  e 

3  3  3  3  3 

W  M  *1  *>  •» 


—  aj  sr  in 

C  E  £  F  E  £ 

*■3  3  3  3  3 

— •  <o  to  « o  to  to 

ii 

ii  ii  ii  ii  ii  a; 

3 

o  ^  in  c 

ineeeee  — 

O'  3  3  3  3  3  *■> 

to  to  to  to  to  C 

O  O 

•o  o 


£  e  *  ~ 

#  3  «- 

**  to  e  a 

w  Gj  w 
C  w  C. 

0.  —  O 
II  x  C  •—  <“* 


D 

w  (0 

<0 

c 

E 

■— * 

E 

C- 

O  II 

II 

D 

D 

c 

II 

D  — 

AJ 

4> 

T3  (0 

<0  «- 

C 

Cl 

C- 

o 

in 

O'  Cl  o 


o  o 

T> 

44 

C 

e 

<0 

a  o 
c  ~ 

s 

- 

*•“  44 

«■* 

TJ  3 

X 

At 

t~  n 

S 

rvj 

o  •»• 

« 

N 

.#» 

o  •- 

(0 

E 

w 

* 

o  44 

X 

+* 

3 

X 

(0  10 

c 

in 

•* 

At 

•4“ 

* 

o 

% 

* 

« 

S  O 

•** 

nj 

/-N 

»#* 

M 

4"\ 

04 

u 

X 

Cf) 

C  w 

«» 

X  *- 

* 

N 

E 

O  x 

X 

w  ® 

E 

O 

-c 

% 

M  E 

1 

** 

3 

w 

4-»  4- 

(ft 

S  L 

Cl 

« 

•«> 

/ 

e 

O 

o 

% 

C 

/-%  ru 

o 

m  z 

-X 

o 

IT 

O  X 

w 

w 

in 

M 

•  e 

•*> 

©  C 

0 

Ol 

E 

% 

O  3 

c 

TJ  <0 

•  • 

X 

j c 

3 

/■* 

•4- 

E 

o 

to 

0) 

• 

4-» 

0>  10 

c 

* 

E 

in 

4-  || 

x  to 

o 

II 

c 

<u 

O 

E 

a 

K  3  »■ 

c.  • 

X 

4-» 

a> 

•  ru 

3 

<0 

3  * 

n  e 

4-* 

X 

<0 

44  (S 

a> 

X 

4^  rvj 

E  3 

c 

~  E 

o* 

X 

* 

0)  At 

V  M 

o 

O 

e 

•-  3 

U.  <0 

w 

N 

Cfc  S 

4>  4 

V) 

X 

©  o 

• 

• 

4 

w  0) 

Cl 

—  ru 

•  • 

o 

o 

o 

N 

c 

C  -X 

C 

c 

o  o 

• 

II 

•• 

4 

o 

II 

o 

II 

•4 

Ui 

U  44 

o 

II  II 

«•» 

3 

to 

<U  -X 

Q 

0) 

ru 

II 

(ft 

o 

o 

o 

o> 

c 

ru  s 

X 

4-* 

in 

c 

a. 

t)  — 

e-* 

0) 

X  X 

N 

N 

c 

O' 

A 

cc 

44  ® 

<0 

E 

E  E 

E 

E 

3 

3 

C  0) 

0) 

3  3 

3 

3 

o 

o 

m 

a. 

*•  L 

c 

■O 

in  <n 

(ft 

a 

•o 

u  uuuuuuuuuuu  uuu  o  o 


140 


V’/v-Vi.v^;i 


sumx4  =  sumx^  ♦  x ( i ) * x ( i ) * x ( i  )  * x C 1  ) 
sumz  =  sumz  ♦,a1oq(z(i)) 
sumzx2  =  sumzx2  t  a  I oq ( z ( i ) ) *x ( i ) *x ( f ) 
points  =  ooints  ♦  t.O 
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